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while  the  ordinates  for  the  others  are  in  terms  of  the  square  root  of  image 
irradiance  (Vl).  Legend: 

_ unaberrated, _ • _ 0.5  waves, _ 1 .0  wave 

4  5  The  amount  of  edge  ringing  is  shown  for  vario'js'amounls  of  the  79 

aberrations  of  (a)  spherical,  (b)  y-axis  coma,  and  (c)  defocus.  Edge  ringing 
is  defined  as  the  ratio  of  the  peak  closest  to  the  edge  to  unity.  The  plots  in 
the  left  column  are  calculated  from  the  edge  image  irradiance  while  the 
others  are  from  the  square  root  of  the  edge  image  irradiance.  The  solid 
curves  are  for  the  unapodised  cases  while  the  dashed  curves  are  for  the 
apodised  cases. 

4  6  The  amount  of  edge  shift  is  shown  for  various  amounts  of  the  aberrations  81 
of  (a)  spherical,  (b)  y*axi$  coma,  and  (c)  defocus.  The  solid  curves  are  for  the 
unapodised  cases  while  the  dashed  curves  are  for  the  apodised  cases.  The 
edge  shift  is  defined  as  the  distance  of  the  half  peak  irradiance  point  from 
the  geometric  edge.  The  plots  in  the  left  column  are  calculated  from  the 
edge  image  irradiance  while  the  others  are  from  the  square  root  of  the 
edge  image  irradiance. 

4  7  The  acutance  of  the  edge  is  shown  for  various  amounts  of  the  aberrations  83 
of  (a)  spherical,  (b)  y-axis  coma,  and  (c)  defocus.  The  solid  curves  are  for  the 
unapodised  cases  while  the  dashed  curves  are  for  the  apodised  cases.  The 
plots  m  the  left  column  are  calculated  from  the  edge  image  irradiance 
while  the  others  are  from  the  square  root  of  the  edge  image  irradiance. 

4.8  The  coherent  image  of  a  slit  through  an  optical  system  having  0.0, 0.5,  85 

and  1.0  waves  of  defocus.  The  top  plot  is  for  an  unapodised  optical  system 
while  the  bottom  plot  is  for  a  system  having  a  Gaussian  apodiser.  Legend; 

_ unaberrated, _ 0.5  waves, _ 1 .0  waves. 

4  9  The  square  root  of  the  irradiance  of  the  coherent  image  of  a  slit  through  86 
an  optical  system  having  0.0,  0.5.  and  1  0  waves  of  defocus.  The  top  plot  is 
for  an  unapodised  optical  system  while  the  bottom  plot  is  for  a  system 
having  a  Gaussian  apodiser  Legend 
_ unaberrated, _ 0  5  waves, _  _ 1  0  waves 
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4  10  The  coherent  image  of  a  slit  through  an  optical  system  having  0.0, 0.5, 
and  1 .0  waves  of  y-axis  coma.  The  top  plot  is  for  an  unapodised  optical 
system  while  the  bottom  plot  is  for  a  system  having  a  Gaussian  apodiser. 
Legend; 

_ unaberrated, _ 0.5  waves, _ _ _ 1 .0  waves. 

4. 1 1  The  square  root  of  the  irradiance  of  the  coherent  image  of  a  slit  through 
an  optical  system  having  0.0, 0.5,  and  1 .0  waves  of  y-axis  coma.  The  top 
plot  is  for  an  unapodised  optical  system  while  the  bottom  plot  is  for  a 
system  having  a  Gaussian  apodiser.  Legend: 

_ unaberrated, _ 0.5  waves, _ 1.0  waves. 

4  12  The  coherent  image  of  a  slit  through  an  optical  system  having  0.0, 0.5, 

and  1.0  waves  of  spherical.  The  top  plot  is  for  an  unapodised  optical  system 
while  the  bottom  plot  is  for  a  system  having  a  Gaussian  apodiser.  Legend: 
_ unaberrated, _ 0  5  waves, _ 1.0  waves 

4  13  The  square  root  of  the  irradiance  of  the  coherent  image  of  a  slit  through 
an  optical  system  having  0.0, 0.5,  and  1.0  waves  of  spherical.  The  top  plot  is 
for  an  unapodised  optical  system  while  the  bottom  plot  is  for  a  system 
having  a  Gaussian  apodiser.  Legend: 

_ unaberrated, _ 0.5  waves, _ 1.0  waves. 

4  14  The  width  of  the  irradiance  distribution  of  the  image  of  a  coherently 
illuminated  slit  is  shown  for  various  amounts  of  the  aberrations  of  (a) 
spherical,  (b)  y-axis  coma,  and  (c) 'defocus.  The  solid  curves  are  for  the 
unapodised  cases  while  the  dashed  curves  are  for  the  apodised  cases. 

Only  slit  images  that  were  in  terms  of  the  square  root  of  irradiance 
were  considered. 


4  1 5  The  irradiance  distribution  in  the  image  of  two  point  sources  separated 
m  the  object  plane  by  (a)  2b  =  20,  (b)  2b  =  30),  and  (c)  2b  =  40.  The  optical 
system  is  unaberrated  and  unapodised.  The  vertical  axes  are  scaled  in  the 
same  relative  units  of  irradiance,  while  the  horizontal  axes  are  scaled  in  the 
same  relative  units  of  dimensionless  distance  coordinates. 

4  1 6  The  irradiance  distribution  in  the  image  of  two  point  sources  separated  in 
the  object  plane  by  (a)  2b  =  20,  (b)  2b  =  30),  and  (c)  2b  =  40.  The  optical 
system  is  unaberrated  and  has  a  Gaussian  apodiser  The  vertical  axes  are 
scaled  m  the  same  relative  units  of  irradiance,  while  the  horizontal  axes 
are  scaled  in  the  same  relative  units  of  dimensionless  distance  coordinates 

4  1 7  The  Sparrow  limit  of  resolution  for  various  values  of  (a)  spherical,  (b)  y-axis 
coma,  and  (c)  defocus.  Data  points  are  not  included  for  higher  amounts  of 
y-axis  coma  and  defocus  because  the  point  image  is  not  sharply  peaked  at 
these  higher  values. 

4. 18  The  mensuration  errorfor  a  coherent  optical  system  having  no  aberrations 
The  solid  curve  is  for  the  unapodised  case  while  the  dashed  curve  is  for  the 
apodised  cases.  The  mensuration  error  is  the  measured  separation  of  the 
two  points  minus  the  expected  separation  of  the  Gaussian  image  points. 
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4.19  The  mensuration  error  for  a  coherent  optical  system  having  0.5  waves  of  100 

spherical  aberration.  The  solid  curve  is  for  the  unapodised  case  while  the 
dashed  curve  is  for  the  apodised  case.  The  mensuration  error  is  the 
measured  separation  of  the  two  points  minus  the  expected  separation  of 
the  Gaussian  image  points. 

5.1  The  experimental  configuration  used  to  measure  the  irradiance  impulse  -  104 

response  of  the  optical  system  formed  by  L1.L2.  the  iris,  and  the  apodiser. 

5  2  The  experimental  configuration  used  to  measure  the  image  irradiance  109 

distribution  of  edge  and  slit  objects  thorugh  the  optical  system  formed  by 
L1.L2,  the  iris,  and  the  apodiser. 

5.3  A  comparison  of  theoretical  and  experimental  results,  (a)  A  single  frame  of  117 
data,  (b)  The  average  of  eight  frames  of  data,  (c)  A  comparison  of  the  data 
in  (b)  to  an  Airy  pattern,  (d)  A  comparison  of  the  data  in  (b)  to  a  central 
slice  through  the  impulse  response  predicted  using  the  measured  aberrations 
The  vertical  axis  for  each  plot  is  relative  irradiance  (I)  plotted  on  a  logarithmic 
scale.  The  horizontal  axis  for  each  plot  is  normalized  distance  v.  The  scales  are 
the  same  on  ail  of  the  plots. 

5  4  The  photograph  in  the  top  of  this  figure  is  of  the  irradiance  impulse  119 

response  of  an  unapodised  and  essentially  unaberrated  optical  system. 

The  aberrations  of  this  system  are  listed  in  column  3  of  Table  6.2.  The  plot 
on  the  bottom  of  this  figure  is  the  calculated  modulus  (vertical  axis)  of  the 
impulse  response  through  the  same  optical  system.  The  horizontal  axes  have 
the  same  units  of  normalized  distance. 

5  5  Theoretical  (solid  curve)  and  experimental  (broken  curves)  plots  of  central  120 
slices  of  the  irradiance  impulse  response  for  an  apodised  (G  =  3)  and 
unaberrated  optical  system. 

5  6  Thephotographinthetopof  this  figure  is  of  the  irradiance  impulse  121 

response  of  an  apodised  (G  =  3)  and  essentially  unaberrated  optical  system. 

The  aberrations  of  this  system  are  listed  in  column  3  of  Table  6.2.  The 
apodiser  is  described  by  (4.6)  and  shown  in  Fig.  3.9  The  plot  on  the  bottom 
of  this  figure  is  the  calculated  modulus  (vertical  axes)  of  the  impulse 
response  through  the  same  optical  system  The  horizontal  axes  have  the 
same  units  of  normalized  distance. 

5.7  An  interferogram  of  the  example  optical  system  generated  by  the  point-  1 23 

diffraction  interferometer.  The  white  dots  are  the  points  which  were 

digitized  and  entered  into  the  interferogram  analysis  program  WI5P. 

5.8  The  photograph  is  of  the  irradiance  impulses  response  of  an  unapodised  1 27 

optical  system  with  the  aberrations  of  astigmatism  =  0. IX,  coma  =  O.IX, 

and  spherical  =  -0.8X.  The  3-D  plot  is  the  theoretically  prediaed  modulus 
(vertical  axis)  based  on  these  aberrations.  The  horizontal  axes  are  in  the 
same  units  of  normalized  distance 
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5.9  Thephotographitoftheirradianceimpulsesresponseofanapodised  128 

(G  =  3)  optical  system  with  the  aberrations  of  astigmatism  s  0.1  X,  coma 

s  O.lX,  and  spherical  s  -0.8X.  The  3-D  plot  is  the  theoretically  predicted 
modulus  (vertical  axis)  based  on  these  aberrations.  The  horizontal  axes  are 
in  the  same  units  of  normalized  distance. 

5. 1 0  Theoretical  (solid  lines)  and  experimental  (dashed  lines)  results  are  1 30 

compared  for  two  systems.  In  the  first  system  (a)  astigmatism  «  0. 1 X, 

coma  s  O.lX,  and  spherical  =  -0.2X;  and  in  the  second  system  (b) 
astigmatism  s  O.lX, coma  s  O.lX,  spherical  s  -O.GX.  The  plots  on  the  left 
are  for  the  unapodised  case  while  the  others  are  for  the  apodised  (G  s  3)  case. 

5  1 1  Theoretical  (solid  lines)  and  experimental  (dashed  lines)  results  are  131 

compared  for  two  systems.  In  the  first  system  (a)  astigmatism  =  O.OX, 
coma  a  0.  f  X,  and  spherical  a  -1.2X;  and  in  the  second  system  (b) 
astigmatism  a  0.2X,coma  a  0.3X,  spherical  a  -1.9X.  The  plots  on  the  left 
are  for  the  unapodised  case  while  the  others  are  for  the  apodised  (G  a  3)  case. 

5  1 2  The  photograph  is  of  the  irradiance  impulses  response  of  an  unapodised  1 32 

optical  system  witht  he  aberrations  of  astigmatism  a  0.2X,coma  a  0.2X, 
and  spherical  a  -O.lX.  The  3-0  plot  is  the  theoretically  predicted  modulus 
(vertical  axis)  based  on  these  aberrations.  The  horizontal  axes  are  in  the 
same  units  of  normalized  distance. 

5  13  The  photograph  is  of  the  irradiance  impulses  response  of  an  apodised  1 33 

(G  a  3)  optical  system  witht  he  aberrations  of  astigmatism  a  0.2X,  coma 
a  0.2X,  and  spherical  a  -O.lX.  The  3-D  plot  is  the  theoretically  prediaed 
modulus  (vertical  axis)  based  on  these  aberrations.  The  horizontal  axes  are 
in  the  same  units  of  normalized  distance 

5. 1 4  Experimental  data  through  the  center  of  the  unapodised  (top)  and  1 35 

apodised  (bottom)  irradiance  impulse  response  of  a  system  with 
astigmatism  a  0.2X.coma  =  0.2X.  and  spherical  =  -O.lX.  Theapodiser 
had  a  value  of  G  =  3. 

5  1 5  Experimental  (dots)  and  theoretical  (solid  lines)  data  are  compared  for  136 

the  experimental  conditions  described  m  Fig.  5  13. 

5  16  Experimental  data  from  the  analog  output  of  the  l-SCAN  linear  detector  138 

oriented  perpendicular  to  the  edge  image  through  an  unapodised  (top  plot) 
and  apodised  (bottom  plot,  G  =  3)  optical  system  with  measured  third 
order  aberrations:  astigmatism  =  O.lX.  coma  =  O.lX,  and  spherical  =  -0  8X 

5  17  Experimental  (dots)  and  theoretical  (solid  lines)  data  are  compared  for  139 

the  experimental  conditions  described  m  Fig.  5  16 

5  18  Experimental  data  from  the  analog  output  of  the  l-SCAN  linear  detector  141 

oriented  perpendicular  to  the  slit  image  through  an  unapodised  (top  plot) 
and  apodised  (bottom  plot,  G  =  3)  optical  system  with  measured  third  order 
aberrations:  astigmatism  =  O.lX,  coma  =  O.lX,  and  spherical  =  -0.8X. 
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5. 1 9  Experimental  (dots)  and  theoretical  (solid  lines)  data  are  compared  for 
the  experimental  conditions  described  in  Fig.  S.18. 

5.20  Experimental  data  from  the  analog  output  of  the  l-SCAN  linear  detector 
oriented  perpendicular  to  the  slit  image  through  an  unapodised  (top  plot) 
and  apodised  (bottom  plot,  G  s  3)  optical  system  with  measured  third  order 
aberrations:  astigmatism  s  O.lX.coma  =  0.1  X,  and  spherical  s  -O.SX. 

5.21  Experimental  (dots)  and  theoretical  (solid  lines)  data  are  compared  for 
the  experimental  conditions  described  in  Fig.  5.20. 

5.22  Experimental  data  from  the  analog  output  of  the  l-SCAN  linear  detector 
oriented  perpendicular  to  the  slit  image  through  an  unapodised  (top  plot) 
and  apodised  (bottom  plot,  G  =  3)  optical  system  with  measured  third  order 
aberrations:  astigmatism  =  0.1  X,  coma  =  0. 1 X,  and  spherical  s  -O.SX. 

5.23  Experimental  (black  dots)  and  theoretical  (solid  lines)  data  are  compared  for 
the  experimental  conditions  described  in  Fig,  5  22. 

5.24  The  irradiance  distribution  in  the  image  of  a  two-point  object  through 
an  essentially  unaberrated  system  which  is  unapodised  in  the  top  plot, 
and  apodised  (G  =  3)  in  the  bottom  plot  in  each  plot,  the  solid  curve 
represents  a  theoretical  prediction  and  the  broken  curve  represents 
experimental  data. 

5.25  The  mensuration  error  is  plotted  relative  to  the  geometrically  expected 
image  point  separation  for  a  system  which  is  essentially  unaberrated 
and  is  unapodised  (solid  curve)  or  has  a  Gaussian  apodiser  (broken  curve). 

5  26  The  mensuration  error  is  plotted  relative  to  the  geometrically  expected 
image  point  separation  for  an  unapodised  (solid  curve)  and  Gaussian 
apodised  (broken  curve).  The  measured  aberrations  were  astigmatism 
=  O.lX.coma  =  0.1  X,  and  spherical  =  -O.SX. 

A2  1  The  mask  used  to  make  the  Gaussian  filters  used  m  this  experiment. 

A2.2  The  characterization  of  filter  #6.  The  onotographs  are  of  the  unapodiseo 
(upper  left)  and  apodised  impulse  responses  of  the  test  system.  Below 
each  photograph  are  plots  of  the  corresponding  theoretical(solid  curve) 
and  experimental  olots  of  irradiance  along  a  line  through  the  center  of 
the  impulse  responses.  Experimental  data  are  from  two  orthogonal  slices. 
The  ordinate  of  each  plot  is  in  terms  of  the  logarithm  of  relative  irradiance(I) 
and  the  abscissa  is  in  terms  of  the  normalized  distance  u. 


A2.3  The  charaaerization  of  filter  #  52.  The  photographs  are  of  the  unapodised 
(upper  left)  and  apodised  impulse  resoonses  of  ^he  test  system.  Below 
each  photograph  are  plots  of  the  corresponding  theoreticai(soliP  curve) 
and  experimental  plots  of  irradiance  along  a  line  through  the  center  of 
the  impulse  responses.  Experimental  data  are  from  two  orthogonal  slices. 
The  ordinate  of  each  plot  is  in  te-'ms  of  ‘he  logarithm  of  relative  irradiancefH 
and  the  aoscissa  is  in  terms  of  the  normanzea  oistance  u. 


xvii/xvili 


LIST  OF  TABLES 


TABLE 

TITLE 

PAGE 

Various  representations  of  the  first  and  third  order 
wavefront  aberrations^^. 

30 

5.1 

The  equipment  used  in  the  experiment. 

105 

5  2 

The  measured  third-order  Seidel  aberrations  after  several 
alignments  of  the  same  optical  system.  The  units  of  the 
tabulated  values  are  wavelengths  (X  =  0.6328  microns) 

115 

CHAPTER  1  INTRODUCTION 


This  dissertation  describes  a  systematic  study  of  coherent  optical  imaging  systems 
which  had  opticaf  aberrations  and  in  which  apodisation  was  employed  to  improve  their 
performance.  Apodisation  has  proven  to  be  an  extreme  r.aortant  method  of 
improving  the  performance  of  incoherent  optical  systems,  including  those  with 
aberrations.  However,  until  this  thesis,  the  benefit  of  applying  apodisation  to  aberrated 
coherent  optical  systems  had  not  been  established.  Since  all  optical  systems  have  at  least 
some  residual  aberrations,  this  study  was  well  motivated. 

In  this  thesis,  the  utility  of  apodisation  tn  improving  the  performance  of  aberrated 
coherent  optical  imaging  systems  was  investigated  both  theoretically  and 
experimentally  The  emphasis  was  on  describing  the  optical  system  performance  when 
(1)  aberrations  were  present  and  (2)  when  apodisation  was  added  to  an  already 
aberrated  system.  Based  on  the  previous  research  concerning  unaberrated  systems,  it 
was  hoped  that  apodisation  would  improve  the  performance  of  aberrated  systems  as 
well  It  will  be  shown  m  the  following  chapters  that  this  was  indeed  the  case  Among 
other  effects,  apodisation  used  in  an  aberrated  coherent  imaging  systems  was  found  to 
be  effective  m  (1)  reducing  the  ringing  in  the  image  of  an  edge  or  a  slit  and  (2) 
decreasing  the  error  in  the  measured  positions  of  the  peaks  in  the  image  of  two  closely- 
spaced  points.  This  behavior  of  an  apodised  coherent  system  had  previously  been 
reported  when  the  system  was  free  from  aberrations.  The  analogous  behavior  when  the 
system  was  aberrated  is  a  new  result. 

Apodisation  is  a  deliberate  modification  of  the  amplitude  transmittance  of  an 
optical  system  As  an  illustrative  example,  consider  a  centered,  circularly-symmetric 
aberration-free,  unapodised.  coherent  optical  system.  The  amplitude  transmittance 
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A(x,y)  of  the  exit  pupil  of  such  a  system  can  be  described  by  the  function 

where  r2  =  x  and  y  are  coordinates  of  the  exit  pupil,  a  is  the  radius  of  the  exit 

pupil,  and  circ  represents  a  function  which  has  a  -  .  3  Afxjr)  *  1  for  rSa  and  A(xj)  =  0 
for  r>a.  The  modulus  of  the  amplitude  distribution  in  the  image  of  a  point  source 
through  this  simple  system  (the  modulus  of  the  amplitude  impulse  response)  is  displayed 
in  Fig  1.1a.  When  squared,  it  yields  the  Airy  diffraction  pattern  or  the  irradiance  impulse 
response  which  is  shown  in  Fig.  1.1c.  If  now  for  example,  Alx,^)  is  modified  by  a 
Gaussian  apodising  function,  the  amplitude  transmittance  becomes 


where  p  is  the  characteristic  width  of  the  Gaussian.  The  Gaussian  is  of  course  only  one  of 
an  infinite  number  of  possible  apodising  functions.  The  modulus  of  the  amplitude 
impulse  response  for  this  system  with  p  =  .58  is  shown  in  Fig.  1.1b.  There  are  two 
immediately  noticeable  effects  when  comparing  Figs. 1.1a  and  1.1b.  First,  the  peaic 
irradiance  has  been  significantly  decreased  by  the  apodisation;  and  second,  the  outer 
rings  evident  in  Fig.  1.1  a  have  been  totally  suppressed  in  Fig.  1.1b.  It  is  apparent  that  the 
term  apodisation  is  appropriate  since  it  is  derived  from  the  Greek:  a  •  meaning  without, 
and  pod  -  meaning  foot.  The  apodisation  has  indeed  removed  the  “feet"  of  the 
amplitude  impulse  response.  This  change  has  important  implications  m  the  analysis  of 
coherent  optical  systems. 

While  incoherent  optical  systems  are  not  withm  the  scope  of  this  study,  it  is  useful 
to  compare  them  to  coherent  optical  systems.  Such  a  comparison  emphasizes  the 
important  characteristics  of  coherent  systems.  Image  forming  systems  have  traditionally 
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Fig  1  1  The  modulus  of  the  amplitude  impulse  response  through  an  optical  system  with 
(a)  a  uniformly  transmitting  apeaure  and  (b)  a  Gaussian  apodiser  in  the  apenure- 
The  irradiance  impulse  response  with  the  uniform  aperture  is  shown  in  (c).  The 
vertical  scales  are  in  relative  units,  with  (a)  and  (b)  having  the  same  scale  The 
other  axes  of  each  plot  have  units  of  relative  distance  and  the  scaling  is  the  same  . 


been  incoherent  because  the  object  to  be  imaged  was  usually  cither,  self-luminous  or 
incoherently  illuminated.  Incoherent  systems  have  the  convenient  property  that  they  are 
linear  in  irradiance,  which  is,  of  course,  the  detected  quantity.  On  the  other  hand,  when 
an  objea  is  illuminated  coherently  the  system  is  linear  in  complex  field  amplitude  and 
hence  nonlinear  in  irradiance.  Partially  coherent  systems  are  also  not  linear  in  irradiance 
but  are  linear  in  the  mutual  intensity  function,  or  for  the  non-quasi-monochromatic  case, 
m  the  mutual  coherence  function. 

Another  important  difference  between  coherent  and  incoherent  systems  is  the 
form  of  their  impulse  response  functions.  This  difference  is  also  illustrated  in  the  plots  of 
Fig. 1.1.  The  incoherent  impulse  response  (Fig. 1.1c)  is  in  terms  of  irradiance  while  the 
coherent  impulse  response  (Fig. 1.1a)  is  in  terms  of  amplitude.  These  figures  graphically 
illustrate  the  fundamental  difference;  i.e.  the  incoherent  impulse  response  is  by 
definition  real  and  positive,  while  this  is  not  necessarily  true  for  the  coherent  impulse 
response.  Thus,  the  convolution  associated  with  the  imaging  process  is  performed  with  a 
fundamentally  different  impulse  response.  This  difference  leads  to  the  deleterious 
effects  which  often  arise  in  coherent  imaging. 

One  such  effect  occurs  in  the  image  of  an  edge.  The  form  of  the  image  of  an  edge 
through  an  optical  system  is  dependent  on  whether  the  illumination  is  coherent  or 
incoherent.  This  is  shown  in  Fig.  1 .2,  where  the  irradiance  of  the  coherent  and  incoherent 
mages  of  an  edge  are  plotted  relative  to  the  geometrical  image  of  that  edge. 

These  curves  were  generated  by  calculating 
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for  the  coherent  image  of  an  edge  and 
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1  2  The  coherent  ar^  incoherent  images  of  an  edge  ar%  plotted  relative  to  the 

geometrical  image  of  the  edge  { _ coherent, _ _ _ incoherent, 

_ geometrical) 


for  the  incoherent  image  of  the  same  edge.  The  quantity  v  is  the  canonical  distance 
coordinate 


where  a  is  the  exit  pupil  radius,  A  is  the  mean  wavelength  of  the  illumination,  d  is  the 
distance  from  the  exit  pupil  to  the  image  plane,  and  y  is  the  image  plane  coordinate 
perpendicular  to  the  edge  image.  The  coherent  image  exhibits  a  pronounced  ringing 
which  is  a  direct  result  of  the  negative  regions  of  the  amplitude  impulse  response.  Also, 
the  position  of  the  coherently  imaged  edge  (the  half  peak-irradiance  point)  is  shifted  to 
the  right  relative  to  the  incoherent  image.  Such  defeas  in  imaging  are  not  limited  to 
edge  images,  but  also  occur  in  images  of  other  simple  objects  such  as  two  points  or  a  slit, 
as  well  as  in  more  complicated  images. 

The  shift  of  the  edge  in  the  coherent  case,  relative  to  the  incoherent  case  is  due  to 
the  nonlinearity  of  the  coherent  imaging  process  The  nonlinearity  arises  when  the  field 
amplitude  in  the  image  plane  is  multiplied  by  its  complex  conjugate  to  give  the  image 
irradiance  distribution.  Thompson'  has  suggested  that  the  analysis  of  coherent  images 
should  be  m  terms  of  the  square  root  of  the  irradiance.  Such  a  comparison  for  the  case 
just  considered  is  shown  in  Fig.  1.3,  where  the  amount  of  edge  ringing  on  the  illuminated 
side  of  the  edge  image  is  reduced  and  there  is  no  shift  m  the  position  of  the  edge.  There 
IS,  however,  an  effective  increase  in  the  amount  of  ringing  on  the  dark  side  of  the 
imaged  edge. 

This  study  was  concerned  exclusively  with  coherent  optical  systems.  Incoherent 
and  partially  coherent  optical  systems  were  not  investigated,  although  occasionally  the 


RELATIVE  IRRADIANCE 


normalized  distance  V 

Fig  1.3  The  irradiance  of  the  incoherent  image  of  an  edge  and  the  square  root  of  the 
irradiance  of  the  coherent  image  of  an  edge  are  plotted  relative  to  the 

geometrical  image  of  the  edge  ( _ coherent, _ . _ . _ incoherent, 

_ _ geometrical).  The 'eft  oramate  refers  to  the  incoherent  image  whiie 

the  right  one  refers  to  the  coherent  image 
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calculated  results  were  compared  to  the  equivalent  incoherent  case.  The  partially 
coherent  case  contains,  as  limits  of  coherence,  both  the  incoherent  and  coherent  cases. 
The  effects  studied  here,  e.g.  edge  ringing,  are  the  most  pronounced  when  the 
illumination  is  fully  coherent.  Thus,  the  limitation  to  coherent  systems  was  a  worst  case 
analysis.  For  that  reason  only  coherent  optical  systems  were  considered. 

The  use  of  apodisation  to  improve  the  performance  of  optical  systems  has  been 
the  subject  of  a  large  number  of  studies.  Most  of  these  have  considered  only  incoherent 
illumination.  The  original  application  of  apodisation  was  to  increase  the  resolution  of 
spectroscopic  instruments.  As  expected,  in  recent  years  the  number  of  studies  of 
aoodisation  in  coherent  systems  has  increased.  In  these  studies,  however,  the  optical 
system  was  always  assumed  to  be  free  of  aberrations.  The  researchers  who  have 
considered  the  use  of  apodisation  in  improving  the  performance  of  an  aberrated  system 
have  limited  themselves  to  incoherent  systems.  •  There  have  only  been  a  few  studies 
concerning  the  effects  of  aberration  on  the  performance  of  coherent  optical  systems. 
There  apparently  have  been  no  studies  analyzing  the  effect  of  apodisation  on  the 
performance  of  aberrated  coherent  optical  systems. 

The  research  described  in  this  dissertation  partially  fills  that  void.  That  is,  the 
oerformance  of  coherent  imaging  systems  was  theoretically  and  experimentally  analyzed 
under  the  conditions:  (1)  no  aberrations  or  apodisation,  (2)  third-order  aberrations  but 
no  aoodisation,  and  (3)  both  abe.rrations  and  aoodisation. 

In  Chapter  2,  the  basic  equations  needed  to  analyze  coherent  imaging  and  beam 
propagation  systems  are  developed  and  previous  research  relating  to  this  problem  is 
reviewed.  Chaoter  3  is  devoted  entirely  to  an  analysis  of  the  amplitude  impulse  resoonse, 
which  is  central  to  the  study  of  all  optical  systems.  In  Chapter  4,  theoretical  results  under 
differing  conditions  of  aberrations  and  aoodisation  m  coherent  imaging  systems  are 


CHAPTER  2  OPTICAL  IMAGING  SYSTEMS  •  ANALYSIS 


Th«  concepts  to  be  studied  in  this  disertetion  ere  emeneble  to  analysis  using  simple 
optical  system  models.  In  the  first  section  of  this  chapter,  these  simple  models  are 
described  and  the  formulae  for  their  analysis  are  developed  from  the  Rayleigh- 
Sommerfeld  theory.  Next,  the  diffraction  theory  of  aberratiom  is  outlined  with 
particular  emphasis  on  the  use  of  Zernike  polynomials.  Then,  the  previous  research 
relating  to  this  problem  is  described  Finally,  the  rationale  for  using  a  Guassian  apodiser 
rather  than  some  other  is  given. 

2.1  THEORETICAL  BASIS  FOR  THE  ANALYSIS  OF  COHERENT  IMAGING  SYSTEMS 

2.1.1  Beam  Propagation  Formulae 

The  optical  systems  considered  m  this  thesis  will  be  modeled  and  analyied  using 
standard  physical  optics  techniques  One  of  the  basic  tools  in  this  field  is  the  beam 
propagation  equation  This  equation  takes  two  forms.  Fresnel  or  Fraunhofer, 


Exh  Pupil  Plane  Observation  Plane 

Fig.2. 1  The  diffraction  geometry  used  m  the  derivation  of  the  beam  propagation 
formulae. 


depending  on  the  distance  of  propagation.  It  approximately  describes  the  optical  field  in 
an  observation  plane  (see  Fig.2.1).  given  that  the  optical  field  is  precisely  known  in  the 
plane  of  the  exit  pupil.  The  two  planes  are  separated  by  the  axial  distance  z;  and  a 
general  point  Pq  in  the  exit  pupil  plane  is  separated  from  a  general  point  P,  in  the 
observation  plane  by  the  distance  To,. 

The  optical  field  at  the  observation  point  P,  can  be  described  by  the  Rayleigh- 
Sommerfeld  diffraction  formula*  (an  application  of  the  Huygens-Fresnel  principle^) 


2a 


where  the  integration  is  over  the  exit  pupil  E  and  it  has  been  assumed  that  ro,  P  X.  The 
cosine  term  is  the  obliquity  factor, 

where  n  is  a  unit  vector  normal  to  the  exit  pupil  plane  pointing  away  from  P,  and  is  the 
vector  joining  the  two  points  Pq  and  P,  and  pointing  towards  P,. 

Following  the  treatment  in  the  text  by  Goodman^,  some  initial  approximations  can 
be  made  to  make  (2.1)  more  amenable  to  calculation.  First,  the  limits  of  integration  are 
made  infinite  with  the  understanding  that,  m  accordance  with  the  Kirchoff  boundary 
conditions,  the  field  m  the  plane  of  the  exit  pupil  Uo(io,>o)  be  zero  outside  the  exit 
pupil  E.  Futhermore,  it  is  assumed  that  the  distance  z  between  the  exit  pupil  and  the 
observation  plane  is  much  greater  than  the  maximum  linear  dimension  m  the  exit  pupil 
Additionally,  it  is  assumed  that  in  the  plane  of  observation  only  a  finite  region  about  the 
optical  axis  is  of  interest,  and  that  the  distance  z  is  much  greater  than  the  maximum 
linear  dimension  of  this  region.  With  these  assumptions  the  obliquity  factor  is  readily 

•Throughout  this  thesis,  quantities  like  U„(P,)  and  UjPo)  which  are  complex-valued  are 
printed  m  bold-^ace  type. 


12 


approximated  by 

“*( '‘•''o,  )  “  “  1>  <2.2) 

where  the  accuracy  is  within  5  percent^  if  the  angle  (n.Fo,)  does  not  exceed  18° 

Similarly,  the  quantity  Tq,  in  the  denominator  of  (2.1)  will  not  differ  significantly 

from  z,  allowing  (2.1)  to  be  rewritten  as* 

“•('’<)  ’  KI 1 1 V*)'  “  “  * 

The  quantity  rg,  in  the  exponent  cannot  be  replaced  with  2  because  any  error  between  r^, 
and  2  would  be  multiplied  by  the  very  large  number  2n/\  with  the  product  easily 
exceeding  2n  radians. 

However,  a  further  simplification  can  be  made  by  applying  the  fresnel 
approximation.  The  distance  rg,  is  given  exaaly  by 


r  = 
ai 


Using  the  binomial  expansion  of  the  square  root,  (2.4)  can  be  written  approximately  as 


r  =2 

01 


(2  5) 


Substituting  (2  5)  mto  (2  3),  performing  the  squaring  ooerations  of  (2  5),  and  rearranging 
terms  results  m  the  Fresnel  propagation  formula 


U  2  ,y  =  - e 

‘V  i\z 


•  ,  '  n  /■  2  2\ 


.liH/,  ,  ^  )  (2.6) 

Xz  V  0  0/  Kz  \  °  ‘  ■  ‘T .  , 


dx  dy  . 
o  •  0 


U  (z  )e  ^ 

^  o\  o'^o; 

This  IS  a  basic  formula  wnicn  relates  the  optical  field  amplitude  m  the  exit  pupn  piane 
U,j(zg,Vo)  to  the  optical  field  amplitude  in  the  observation  plane  U,(z„^,). 


*in  this  thesis  a  double  integral  with  infinite  limits  on  both  integrals  is  represented  as  m 


The  quantity  Uo(zo.>o)  contains  all  of  the  information  about  the  exit  pupil, 
including  aberrations  and  apodisation  That  is 

*hx  f  \ 

(2.71 


where  Uin(ZoO'o)  tf’*  wavefront  incident  on  the  exit  pupil,  Afroj-,,)  represents  the 
amplitude  transmittance  of  the  exit  pupil.  ♦(zoO'o)  represents  the  phase  aberrations 
associated  with  the  exit  pupil,  and  B(zo.J'o)  represents  the  finite  extent  of  the  exit  pupil 
Alx„.Vo)  is  the  term  which  describes  the  apodiser  In  general.  A[z^„y^,)  could  be  complex; 
but  m  this  study  only  real  apodisations  will  be  considered 

It  IS  instructive  to  note  that,  aside  from  multiplicative  amplitude  and  phase  faaors 
in  (2  6)  that  are  independent  of  (Xo.>o).  tbe  function  U,(z,,y,)  may  be  found  from  a  Fourier 
transform  of 


U 


(28) 


where  the  transform  must  be  evaluated  at  the  spatial  frequencies  fi  =  i|.X:  and  A  =y,..l; 
to  assure  the  correct  spatial  scaling  in  the  observation  plane 

In  the  near-field  analysis  of  an  optical  system,  (2  6)  is  used  directly  The  field 
amplitude  in  the  exit  pupil,  Uo(xo.To)  3S  stated  previously,  limited  in  spatial  extent  by 
that  exit  pupil  anc  contains  the  incident  field  amplitude,  the  attentuation  by  tne 
apodiser  and  any  aberrations  which  may  be  part  of  the  illuminating  wavefront.  The  field 
amplitude  at  each  point  in  the  observation  plane  can  then  be  determined  by  performing 
the  indicated  integration 

For  a  far-fieid  analysis,  (2.6)  simplifies  further  If,  m  addition  to  the  assumptions 


made  in  the  Fresnel  propagation  case,  the  stronger  (Fraunhofer)  assumption 


is  adopted,  then  the  quadratic  phase  factor  inside  the  integral  is  approximately  unity.  In 
this  case,  the  observed  field  amplitude  can  be  found  directly  as  a  Fourier  transform  of  the 
field  amplitude  of  the  exit  pupil.  Thus,  in  the  region  of  Fraunhofer  diffraction, 


U c  U  '  ‘ (  (*  u  f  '1  i 

j  =  IIT' 


(2.10) 


Aside  from  the  multiplicative  factors  preceding  the  integral,  this  equation  states  that 
is  the  Fourier  transform  of  the  exit  pupil  distribution  Uo(xo.>o).  evaluated  at  the 
spatial  frequencies  ft  =  xjkz  and  fy=y,/Xz. 

The  analysis  of  the  focussed  beam  propagation  system  is  likewise  simplified  by  the 
presence  of  the  focussing  lens.  A  perfect  lens  only  changes  the  curvature  of  the  incident 
wavefront,  i.e.,  the  field  just  after  the  lens,  Uo’(io..yo).  'S  given  by 

where  U,„(io.,yo)  is  the  field  just  before  the  lens  and  f  is  the  focal  length  of  the  lens, 
assuming  that  the  lens  is  in  the  exit  pupil. 

If  the  observation  plane  is  a  distance  .'"away  from  the  exit  pupil,  then  the  quaorat.c 
phase  factors  inside  the  integral  of  (2.6)  cancel  in  this  case,  the  Integral  is  agam  an  exact 
'ourier  transform  of  the  exit  puoil  function 

‘  X  '  "  /  2  2\  i2n/  \  . 
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2.1.2  Imaging  Systems 


The  analysis  of  imaging  systems  will  be  based  on  the  techmaues  of  ^ouner 
transforms  These  methods  were  developed  primarily  by  Duffieux^  partly  m 


collaboration  with  Lansraux*.  The  following  development  follows  closely  that  of  Born 
and  Wolf5,  Section  9.5.1. 


Xo  Xj 


Fig. 2. 2  The  geometry  of  a  centered  optical  system  which  gives  rise  to  an  aberrated 

wavefront  W.  The  aberration  funaion  ♦  is  in  terms  of  the  distance  QQ*  between 
W  and  the  reference  sphere  S 


Consider  the  centered  optical  system  sketched  m  Fig. 2  2  Points  m  the  image  plane 
and  m  the  plane  of  the  exit  pupil  are  specified  by  Cartesian  coordinates  and  {{,.r|). 
respectively.  The  yo.yi,  and  q  axes  while  not  shown  in  Fig. 2  2  are  perpendicular  to  the 
optical  axis  and  to  the  lo.ii,  and  ^  axes,  respeaively.  Points  m  the  object  plane  are 
specified  by  scale  normalized  coordinates  (Xo,yo)  ^uch  that,  if  Pq  at  {Xq.Yo)  is  a  typical 
objea  point  and  is  the  lateral  magnification,  then  i,  =  MX.j  and  =  MY^  so  that  the 
object  point  and  its  image  have  the  same  coordinate  values 

The  imaging  properties  of  the  optical  system  are  specified  by  the  transmission 


function  K  (lo.v^;  defined  as  'he  comolex  amplitude,  per  unit  area  of  'he  object 
plane,  at  the  point  (z„y,)  in  the  image  plane,  due  to  a  disturbance  of  unit  amplitude  and 
zero  phase  at  the  object  point  ixo.yo).  The  imaging  process  is  then  described  by 


(2.12) 


where  Uo(xo.Vo)  represents  the  complex  disturbance  in  the  object  plane  and  U,(x,o-,)  'S 
’»'e  resultant  optical  field  amplitude  at  the  point  in  the  image  plane. 

The  transmission  function  K  x.jr,)  depends  on  the  physical  properties  of  the 
optical  system.  This  dependence  can  be  derived  by  determining  the  image  of  a  point 
object  of  unit  amplitude  and  zero  phase  at  the  point  Xg  =  ZQ,  yo=yo-  case,  the 

object  is  described  by 

where  5  IS  the  Dirac  delta  function.  Then  (2.12)  reduces  to 

that  IS,  the  transmission  function  K  is  equal  to  the  disturbance  in  the  image  plane  due  to 
the  point  source  (2. 1 3)  in  the  object  plane. 

The  Gaussian  reference  sphere  S  is  centered  on  the  image  point  x,’  =  Xo'.  y,'  =  .v„' 
(shown  as  P,  in  Fig  2  2)  Let  R  be  the  radius  of  the  reference  sphere  and  let 

)  =  I  "  r-'o 

oe  the  disturbance  at  a  typicai  point  :£„.t)  on  this  sphere  due  to  the  point  source  (2.13). 
Apart  from  an  additive  factor  n/2,  the  phase  of  G  is  then  the  aberration  function  <l>  and 
the  amplitude  of  G  is  a  measure  of  the  non-uniformity  in  the  amplitude  of  the  image¬ 
forming  wave.  The  factor  t/A  on  the  right  side  of  (2.15)  was  introduced  to  simplify  later 
formulae 

Using  :ne  Sayieign-Sommerreio  oirfraction  formu.a  i2.1),  the  optical  field  m  tne 


image  plane  is  (small  angles  of  diffraction  assumed) 
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(2.16) 


where  s  is  the  distance  from  the  point  {{..n)  on  the  reference  sphere  to  the  point  (i.o'c)  'o 
the  image  plane;  and  the  integral  extends  over  that  portion  of  the  reference  sphere 
which  approximately  fills  the  exit  pupil. 

The  distance  5  is  given  exactly  by 


s  = 


and  approximately  by 
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(2.17) 


(2.18) 


where  the  binomial  expansion  has  been  applied  to  (2.17),  all  but  the  first  three  terms 
have  been  dropped,  and  the  terms  which  are  quadradic  in  object  or  image  plane 
coordinate  have  been  discarded  as  being  much  smaller  than  R. 

^hen,  the  combination  of  (2. 14)-{2  18)  results  m 
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dWn . 


(2.19) 


where  G  iS  taicen  to  oe  equal  to  zero  at  oomts  i.:\i  outsioe  the  opening  m  the  exit  puoil. 
This  IS  the  required  relation  between  the  transmission  function  K  and  the  pupil  function 
C  of  the  system  Since  K  may  be  regaroeo  as  the  disturbance  m  the  image  of  a  pomt 
tou'ce  t  can  -o/v  oe  oentifieo  also  as  .^-'o'l’.-oe  mpuise  'esponse  of  the' optical 


system  From  (2  19)  it  can  be  seen  that  this  amplitude  impulse  response  is  proportional  to 


the  Fourier  transform  of  the  pupil  function  G,  where  the  transform  must  be  evaluated  at 
the  spatial  frequencies /i  =  (z,-Xo)/\/?  and  4=0'i->oV^- 

In  a  well -corrected  circularly  symmetric  system,  K  approximates,  apart  from  a 
constant  factor,  the  amplitude  of  the  Airy  diffraction  pattern.  Since  this  study  involves 
the  performance  of  optical  systems  which  are  not  very  well-corrected,  the  form  of  K  will 
depart,  sometimes  significantly,  from  the  Airy  pattern.  However,  considered  as  a 
function  of  (zo.yo).  the  transmission  function  varies  slowly  as  this  point  explores  the 
object  surface.  More  precisely,  the  working  field  may  be  divided  into  regions,  where 
each  region  is  large  compared  to  the  finest  detail  that  the  system  can  resolve;  and  with 
the  property  that  in  each  such  region  P,  K  is  to  a  good  approximation,  a  funnion  of  the 
displacement  vector  from  the  image  point,  but  not  of  the  image  point  itself.  In  such 
cases, 


(2,20) 


A  region  P  with  this  property  is  said  to  be  an  isoplanatic  region  of  the  system.  This  study 
will  only  consider  objects  that  are  so  small  as  to  fall  within  such  an  isoplanatic  region  in 
this  case.  (2. 12)  and  (2.19)  may  be  replaced  by 


i  -  .7  1  - 
I V  I  J 


J  \  J 

U  ;  i  .  V  IKi  I  —  z  .V  -V  idz  av 

o\  0^0/  \  I  0-^1  -o/  u-o 


(2.21) 


and 
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dWn  . 


(2  22 


where  the  function  G  is  now  independent  of  the  object  point. 
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If  Uo(xo.>o).  afitl  K  are  represented  as  the  Fourier  integrals 

^  °  °Ufdg  . 


U  X  O'  1  = 

0\  0^0  ‘ 


•’t2n| 


u  (/;«)? 


(2.23a) 


dfdg  . 


'.en  Dy  tne  Fourier  inversion  formula 


u  j  f,g  1  = 


U  (  X  o  )«?' 

o\  0  I 


2n 


dx  d\ 

0  ’  0 


(2.23b) 


(2,23c) 


(2.24a) 
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K(  I  O'  je  did\ 


(2  24c) 


Equation  (2  21)  represents  the  imaging  process  as  a  convolution  Substituting 
;2.24a'C)  into  the  -ouner  inverse  of  (2.21)  and  aoDlying  '.he  convolution  theorem^  i-esuits 
in 


u 


f.8)  =  ^if.g  j^(j'-g  ) 


(2  25) 


This  equation  implies  that  if  the  optical  fields  m  the  object  ani^image  planes  are 
each  considered  as  a  superposition  of  spatial  frequencies  f  and  g,  then  each  component 


of  spatial  frequency  in  the  image  plane  depends  only  on  the  corresponding  component 
in  the  object  plane,  and  the  ratio  of  components  is  k.  A  comparison  of  (2.23c)  and  (2  22) 
reveals  that 

so  that  the  frequency  response  function  is  equal  to  the  value  of  the  pupil  funaion 
G  at  the  point 

i,=  \Rfandi\  =  \Rg. 

The  calculation  of  the  image  of  any  appropriately  small  object  can  now  be 
accomphshecJ  by  Founer  transform  techniques.  That  is, 

where  F{  }  and  f-'(  }  indicate  the  direct  and  inverse  Fourier  transforms,  respectively,  of 
the  bracketed  terms.  The  inverse  Fourier  transform  must  be  evaluated  at  the  spatial 
frequencies  f=l/kR  and  g  =  r\/XR  and  the  direct  Founer  transform  must  be  evaluated  at 
the  soatial  frequencies  f -  (x^-x^V^R  and  g  =  (^\-yo)l\R  This  powerful  tool,  (2.28),  will  be 
-sec  extensively  n  the  analysis  of  aberrated  imaging  systems  m  this  thesis 

2.1.3  Measures  of  Performance 

2. 1.3.1  Two-Point  Resolution 

The  performance  of  an  optical  system  is  measured  using  criteria  that  are 
deoencent  upon  the  intended  use  of  the  system  in  the  telescope,  for  example,  the 
resolution  of  two  closely-soaced  ooint  objects  -s  imoortant.  Lord  Ravleigh®  develooed 
the  first  resolution  criterion,  which  now  bears  his  name  iJsmg  Airy's^  result  concerning 
the  mage  of  a  oomt  source  through  a  circular  aoemure.  Lord  Ravle'Oh  stated  ’hat  two 
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point  sources  are  just  resolved  if  the  maximum  of  one  irradiance  pattern  coincides  with 
the  first  minimum  of  the  other.  As  BarakatS  has  pointed  out,  the  Rayleigh  criterion  is 
based  on  the  tacit  assumption  that  the  two  point  sources  are  mutually  incoherent.  It  is 
necessary  of  course  that  they  also  be  of  equal  irradiance.  Born  and  Wolf5  extended 
Rayleigh's  criterion  by  stating  that  resolution  is  obtained  when  the  irradiance,  midway 
between  the  geometrical  images  is  19%  less  than  the  maximum  irradiance  in  the  image 
of  the  two  points. 

An  alternate  criterion,  suggested  by  Sparrow^,  states  that  the  two  point  sources 
are  ]ust  resolved  if  the  second  derivative  of  their  image  irradiance  vanishes  at  a  point 
midway  between  the  Gaussian  image  points.  The  Sparrow  criterion  is  more  readily 
amenable  to  quantitative  calculations  and  has  often  been  used  m  assessing  systems.  It  is 
also  well  suited  to  analyzing  resolution  in  coherent  systems.  The  Sparrow  criterion  will 
be  used  m  this  study. 

The  spatial  coherence  of  the  illumination  from,  two  point  sources  can  dramatically 
affect  the  resolution  of  their  images.  This  is  illustrated  in  Fig. 2. 3,  where  Grimes  and 
Thompson '0  have  calculated  the  irradiance  distribution  m  the  image  of  two  point  objects 
for  various  values  of  y.  the  complex  degree  of  coherence  between  the  two  point  sources 
The  geometrical  separation  of  the  points. which  is  just  greater  than  the  Rayleigh  limit,  is 
moicaied  Dy  the  two  vertical  dashed  lines  "he  two  points  are  well  resolved  m  the 
mcoherent  limit,  the  resolution  decreasing  with  increasing  v  until  when  y  =  0  6  the  curve 
IS  essentially  flat  midway  between  the  Gaussian  image  points  (the  Sparrow  resolution 
limit)  For  greater  values  of  y,  the  two  points  are  unresolved.  Rojak"  has  calculated  the 
limit  of  resolution  using  the  Sparrow  criterion  for  values  of  y  from  the  coherent  limit  to 
the  incoherent  limit.  '• 
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2. 1.3.2  Mensuration  Error 


It  is  apparent  from  Fig. 2. 3  that  the  separation  of  the  two  points  when  measured  m 
the  image  plane  depends  on  the  degree  of  coherence  The  difference  between  the 
expected  separation  based  on  geometrical  considerations  and  the  measured  separation 
is  termed  the  mensuration  error.  This  mensuration  error  for  a  fixed  degree  of  coherence 
depends  also  on  the  separation  of  the  point  sources  in  the  object  plane'O.  In  the  arsalysis 
of  the  image  of  two  point  sources,  the  only  measurable  quantity  is  the  separation  of  the 
two  peaks  in  the  resultant  irradiance  distribution,  a  separation  that  would  normally  be 
considered  to  be  the  real  separation  of  the  object  points  Hence,  the  mensuration  error  % 
a  very  important  measure  of  performance  in  some  optical  imaging  systems 

2.1 .3.3  Edge  Ringing,  Shift,  and  Acutance 

It  is  important  in  other  systems  to  accurately  image  an  edge  object.  Two  major 
effects  encountered  are  edge-ringing  and  edge-shifting.  These  effects  have  been 
investigated  by,  among  others,  Hopkins’^,  SteeP^,  Canals-Frau  and  Rosseau’*,  Kinrly’S_ 
5kinner’6,i7_  Considine’8_  and  Thompson'^.  These  effects  are  illustrated  In  Fig,  l  2  where 
the  theoretical  irradiance  distribution  of  the  image  of  an  edge  object  by  a  one- 
dimensionai  ideal  lens  is  shown  for  both  coherent  arc  incoherent  illumination  It  is  also 
aooarent  from  this  ‘igure  that  the  slooe  of  the  irraciance  near  the  edge  is  diffei-ent  'or 
the  two  cases.  The  slope  of  the  edge  irradiance  is  termed  the  acutance  and  is  greater  for 
a  coherent  image  than  for  an  incoherent  image  of  the  same  edge  with  the  same  optical 
system.  When  quantified,  these  effects  can  also  be  measures  of  performance  of  the 
imaging  system. 

2. 1.3. 4  Slit  Ringing  and  Width 
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The  use  of  bar  targets  for  measuring  resolution  of  optical  systems  is  perhaps  the 
most  common  ana  rapid  method  of  evaluating  an  optical  system.  The  OTF  of  an 


incoherent  system  can  quickly  be  estimated,  using  an  object  consisting  of  many  sets  of 
bar  targets,  each  having  a  different  fundamental  spatial  frequency.  Although  the  OTF  is 
not  applicable  to  coherent  systems,  it  is  interesting  to  compare  the  image  of  a  slit  (or  a 
bar)  in  the  coherent  and  incoherent  limits.  When  coherently  illuminated,  the  edges  of 
the  slit  will  exhibit  edge-ringing  and  edge-shifting.  As  an  illustration  of  this  problem,  the 
coherent  image  of  a  slit  is  displayed  in  Fig.2.4.  Under  coherent  illumination  the  slit 
appears  to  be  narrower  than  under  incoherent  illumination.  The  opposite  situation 
would  occur  for  a  coherently  imaged  bar;  the  bar  would  appear  to  be  wider.  The  size  of 
"the  bar  (or  slit)  can  be  used  as  a  measure  of  the  performance  of  the  system.  The 
normalized  distance  v  used  m  this  figure  was  defined  in  Chapter  1. 

A  more  meaningful  depiction  of  the  coherent  image  of  an  edge  might  be  a  plot  of 
the  square  root  of  the  irradiance.  This  follows  the  suggestion  made  by  Thompson’  Fig. 
2.5  is  such  a  plot.  The  advantage  of  this  approach  is  the  removal  of  the  nonlinearity  in 
the  coherent  imaging  process.  One  result  is  that  there  is  no  edge  shift;  another  is  that 
the  edge  ringing  is  reduced 

2. 1.3. 5  Streh!  Ratio,  Encircled  Energy,  Second  Moment,  and  Near-field 
Fluctuations 

In  the  case  of  laser  systems,  the  important  parameters  of  the  focussed  beam  are 
the  Strehl  ratio,  encircled  energy,  and  second  moment  of  the  irradiance  distribution  The 
Strehl  ratio  is  formed  by  the  ratio  of  the  peak  irradiance  m  a  particular  plane  of  the  focal 
region  of  an  aberrated  beam  to  the  peak  irradiance  in  the  focal  region  of  an  unaberrated 
oearr  The  encircled  energy  is  the  total  energy  within  a  given  circle  in  a  plane  normal  to 
the  beam  The  second  moment  refers  to  the  variance  of  the  irradiance  distribution  of  the 
beam  Each  of  these  can  be  used  as  a  measurement  of  performance.  In  the  near  field  of 
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Fig. 2.4  The  coherent  image  of  a  slit  (solid  curve)  through  an  unaberrated  and 
unapodised  optical  system  having  a  circular  exit  pupil.  The  geometrical  image  of 
the  slit  is  shown  as  a  dashed  curve. 
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Fig  2.5  The  solid  line  is  a  plot  of  the  square  root  of  the  irradiance  distribution  of  the 
image  of  a  slit  through  a  circularly  symmetric,  unaberrated,  unaoodised  coherent 
optical  system.  The  dashed  line  is  the  geometrical  image  of  that  same  edge. 


the  beam  propagation  system,  the  irradiance  distribution  can  have  significant 
fluctuations  that  arise  from  diffraction  caused  by  the  apertures  in  the  system  These 
fluctuations  can  degrade  the  far  field  performance  of  the  system  by  distorting  reflecting 
surfaces  as  pointed  out  by  Avizonis,  et  al^o.  The  fluctuations  could  also  damage 
components  of  the  final  amplifier  stages  in  systems  such  as  a  high  peak  power  pulsed 
laser.  Thus,  the  amount  of  the  near-field  fluctuations  can  be  an  important  limitation  of 
performance  in  these  systems. 

2.2  THE  DIFFRACTION  THEORY  OF  ABERRATIONS 

Both  the  geometrical-  and  diffraction-based  theories  of  aberrations  are  covered 
comprehensively  in  the  text  by  Born  and  Wolf^  Those  parts  of  their  treatment  which  are 
pertinent  to  this  study  will  be  briefly  reviewed. 

Consider  again  the  centered  optical  system  shown  in  Fig. 2. 2  with  a  monchromatic 
point  source  Pq  in  the  object  plane.  The  Gaussian  image  point  P,  m  the  image  plane  is  the 
center  of  the  Gaussian  reference  sphere  S  which  intersects  the  exit  pupil  at  the  on-axis 
point  C.  A  real  wavefront  W  arising  from  the  point  source  Pq  will  m  general  not  comciae 
with  the  reference  sphere  but  will  deviate  from  it,  perhaps  as  shown  m  Fig  2.2  The  real 
wavefront  :s  drawn  intersecting  the  reference  sphere  at  the  axial  point  C. 

At  an  aroitrary  ooint  Q  on  Al  a  'ine  pe-'cenoicuiar  'o  W  i-eoresents  an  ootical  ray,  or 
wave  normal.  In  general,  this  wave  normal  will  intersect  the  image  plane  at  a  point  P,* 
wnicn  s  aifferent  from  The  separation  oi  •nese  two  oomts  P,  P,*  s  the  wave  normal 
aberration  while  the  deformation  of  W  relative  to  S  along  the  ray  QQ*  is  the  wavefront 
aoerration.  The  deformation  of  W  m  the  region  of  the  exit  pupil  will  be  described  by  the 
aoerration  runction  'P.  Assuming  tnai  tne  retractive  :nae.x  or  tne  image  space  is  unity,  'P 
(taken  as  positive  in  Fig. 2. 2)  represents  the  distance  QQ*  along  the  wave  normal 
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In  this  study,  the  aberration  function  will  be  expressed  in  terms  of  both  the  Zernike 
polynomials  and  the  Siedel  expansion.  The  Zernike  polynomials  are  a  complete  set  of 
polynomials  which  are  orthogonal  over  a  unit  circle.  They  were  introduced  by  Zernike^' 
in  his  famous  paper  on  the  phase  contrast  microscope  and  have  been  used  extensively  by 
many  researchers  since  then. 

The  circle  polynomials  of  Zernike  are  polynomials  in  the  real  variables  x 

and  y.  In  terms  of  the  polar  coordinates,  x  =  ^  sin©  and  y  =  ;;  cos0,  they  are  expressed  as 


V ^  I  ^in0,<^co.sO 


(2  29) 


where  m  >  0orm>0andn  >0  are  integers,  nijm],  and  n-lm|  is  even.  These  polynomials 
satisfy  the  orthonormal  relation 


,  ,  V'"*(i.>)v'"(i,.v)dxd>=  —  5  6  , 

,2^  2  ^  j  n  \  ■'  j  n\  ■  )  mm  nn  ’ 


(2.30) 


where  5j|  is  the  Kronecker  delta  function  and  the  asterisk  denotes  the  complex  conjugate. 

The  radial  functions  R„'"(p)  are  polynominals  in  c  containing  powers  in  Cp,  ("-2 . (hi 

Instead  of  the  complex  polynomials  V',,-",  one  may  use  the  real  polynomials 
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The  wavefront  can  then  be  expressed  as 

k  n 

'tHp.O)  =  ^  ^  cos  m9  +  D  sinmQ),  (2.32) 

n  nm  n  rn 

n  =  <3  m  ®  o 

where  Cr,rr,  and  Dr,rn  ^re  the  coefficients  denoting  the  amount  of  a  particular  aberration 


preser^t 


In  this  study  the  polynomial  R„'"  (p)  cos  m4>  is  sufficient  to  describe  the  aberrations 
because  the  system  is  assumed  to  be  symmetrical  about  the  meridional  plane  (4>  =  0). 
Table  2,1  lists  the  primary  wavefront  aberrations  as  well  as  lower  order  wavefront 
deformations.  The  first  column  of  the  table  lists  the  third  order  and  lower  wavefront 
aberrations.  The  second  and  third  columns  list  the  corresponding  n  and  m  values  used  in 
the  equations  (2.29)  through  (2.32).  The  fourth  column  lists  the  specific  functional  form 
of  each  of  these  aberrations  in  terms  of  the  Zernike  polynomials  of  (2  32)  The  fifth 
column  lists  these  same  aberration  terms  where  a  conversion  from  the  variables  p  and  0 
to  the  variables  x  and  y  have  been  accomplished.  In  this  case  the  wavefront  is  expressed 
as 

^  1 

4>U.y)=  X  S  •  (2.32) 

i*o  7=0 

where  fly  is  the  coefficient  controlling  the  amount  of  a  particular  aberration  which  is 
present.  In  the  last  column,  the  aberration  s  are  expressed  as  terms  from  the  expansion 

k  n 

^(p,9)=  p"  (a  cos^9 -r  6  sin‘0),  (^  -•3) 

^  nt  nl 

n  =  0  /  =  o 

where  the  coefficients  a„/  and  b„i  are  commonly  called  the  Seidel  coefficients  The 
difference  between  this  expansion  and  the  Zernike  expansion  of  (2.32)  is  that  the 
wavefront  is  exoressed  in  terms  of  powers  of  cos0  and  sin9  while  in  the  Zernike 
expansion  the  wavefront  is  expressed  m  terms  of  multiples  of  the  polar  angle  9 

The  entries  for  x-tilt,  y-tilt,  and  piston  are  relevant  for  interferometric 
measurement  and  interpretation  of  wavefronts  Piston  Is  a  term  referring  to  amounts  of 
uniform  phase  across  the  exit  pupil. 


Varioui  repiesenjaiiofis  of  the  first  and  third  order  wavefront  aberrations^*. 


2.3  PREVIOUS  RESEARCH  •  APODISATION  TO  IMPROVE  SYSTEM  PERFORMANCE 


The  technique  of  suppressing  the  sidelobes  in  the  irradiance  distribution  of  the 
far-field  diffraction  oattern  of  a  finite  aperture  was  the  original  definition  of 
apodisation.  Me  idly  defined  now,  apodisation  refers  to  any  method  which 
changes  the  performance  of  an  optical  system  by  modifying  the  amplitude  and  phase 
transmittance  of  the  exit  pupil.  It  has  been  shown  by  many  researchers  that  apodisation 
can  indeed  improve  the  imaging  properties  of  an  optical  system.  Most  of  the  work  has 
assumed  an  aberration-free  system,  but  recently  there  has  been  some  interest  in 
aberrated  systems  as  well.  However,  there  are  apparently  no  studies  concerning  the  use 
of  apodisation  to  improve  the  performance  of  aberrated  coherent  optical  systems. 

2.3.1  Studies  of  Systems  without  Aberrations 

2.3. 1.1  Incoherent  Optical  Systems 

In  retrospect,  the  most  common  apodisation  seen  m  optical  systems  is  the  central 
obscuration.  Lord  Rayleigh^i  was  the  first  to  point  out  that  a  circular  stop  in  the  center 
of  a  circular  aperture  causes  the  central  maximum  of  the  Airy  pattern  to  become 
narrower  and  also  increases  the  deoth  of  focus.  Straube!--  aooears  to  have  been  the  ^>s; 
to  study  effeas  of  an  apodiser  other  than  the  central  obscuration.  Most  studies  since 
toe''  nave  attemoteo  'o  aotimize  some  asoec  ot  oerrormance  of  an  aoticai  svstem  os 
finding  an  appropriate  apodismg  filter  Toraldo  di  rrancia-"*  originated  the  concept  of 
superresoiution  over  a  limited  field  This  occurs  when  the  diffraction  pattern  is  forced  to 
assume  a  senes  of  zeroes,  the  ‘irst  one  ceing  such  that  tne  central  core  s  narrower  than 
that  for  the  uniform  pupil,  while  the  remaining  zeroes  are  spaced  around  the  core  in 
such  a  way  'hat  a  darK  zone  of  limited  extent  and  axhibiting  'ow  secondary  lobes  s 


produced.  Luneberg25  used  the  calculus  of  variations  to  determine  the  pupil  function 
given  a  prespecified  distribution  of  image  irradiance  In  two  comprehensive  review 
articles,  Wolf26(l951)  and  later  Jacquinot  and  Roizen-Dossier27(i964)  described  the 
previous  /.ork.  Since  then,  many  workers  (see  for  example  references  28-31)  have 
coni.-  I .  . :  :o  this  field.  Barakata-32.33  has  been  especially  productive,  solving  three  of 
the  four  apodisation  problems  originally  posed  by  LunebergiS. 

2. 3. 1.2  Partially  Coherent  Optical  Systems 

Studies  using  partially  coherent  illumination  are  motivated  by  the  physical  fact 
that  real  systems  are  seldom  either  fully  coherent  or  fully  incoherent.  McKechnie^^  has 
shown  that,  by  obstucting  the  central  part  of  the  condenser,  the  resolution  of  a 
microscope  may  be  improved.  In  a  similar  study,  Nayyar  and  Verma^®  have  calculated  the 
limit  of  resolution  (a  modified  Rayleigh  criterion)  for  opaque  annular  and  n-phase  semi¬ 
transparent  annular  apertures  of  two  partially  coherent  point  objects.  Other  studies 
concerning  two  point  objects,  using  different  forms  of  apodisation,  were  conducted  by 
Magiera  and  Pietraszkiewicz^s,  Nayyar  and  Verma^?,  and  Mehta^B  39.  Kmtne'’  and 
Sillitto'’^’  have  found  a  general  condition  for  aoodisers  which  insure  suppression  of  edge¬ 
ringing.  Som  and  Biswas**’  have  studied  the  performance,  m  partially  coherent  light,  of 
apcdisers  wnicn  we''e  optimized  for  use  m  ccnerent  illumination. 

2. 3. 1.3  Coherent  Ootical  Systems 

The  most  noticeable  defect  m  the  coherent  image  of  an  edge  object  is  the 
phenomena  of  edge-ringing.  The  use  of  apodisers  to  suppress  edge-ringing  has  been 
shown  to  oe  successful  oy  Araki  ano  Asakura**-  ;<,o  .^ao.et  ai-^  ■*B,  Smith-B,  .eaver  and 
Smith*^,  and  Thompson  and  Krisl^s  However,  it  is  noted  in  several  of  these  studies  that 
3s  edae  ringing  s  suppressed  the  amount  of  eoge  shift  s  increased,  ~he  A/orv  of  Atom 


and  Asakura*^  is  particularly  interesting  because  they  have  fournj  an  optimum  apodiser 
which  suppresses  edge-ringing  but  introduces  a  minimum  amount  of  additional  edge 
shift.  In  related  studies.  Leaver^?  has  investigated  the  image  of  a  line  objea  and  Rao  et 
also  have  studied  the  images  of  trapazoidal  objects. 

The  improvement  of  the  two-point  resolution  of  a  coherent  optical  system  has 
been  approached  from  several  viewpoints.  Barakat^  found  apodisers  which  improved 
the  two-point  resolution  according  to  the  Sparrow  criterion.  Thompsons'  has 
investigated  the  diffraaion  by  annular  apeaures  which  have  central  regions  that  are 
semi-transparent  and  which  add  a  uniform  phase  Such  apertures  will  result  in  an 
improved  two-point  resolution.  Wilkins52  has  solved  the  apodisation  problem  of 
determining  the  diffraaion  pattern  which  has  a  specified  Sparrow  limit  of  resolution  and 
the  maximum  possible  Strehl  ratio.  Clements  and  WilkinsSS  have  solved  the  similar 
problem  of  maximizing  the  encircled  energy  with  a  specified  Rayleigh  resolution  limit. 

The  work  of  Wilkins  and  ClementsS2.53  applies  not  only  to  image  forming  systems 
but  also  to  beam  propagation  systems.  In  a  related  study,  Hazra^^  also  has  found  a  class 
of  oDtimum  apodisers  which  will  maximize  the  encircled  energy.  Rao  et  al55  have 
investigated  the  effect  of  a  specific  type  of  apodiser  on  several  measures  of  far-fie'c 
oerformance. 

The  performance  of  apodisers  in  the  near  field  has  been  studied  by  Hadley56  as  well  as 
"hompson  and  Krisl'*^  These  workers  have  shown  that  aoodisation  can  decrease  ’he 
peaks  of  the  Fresnel  diffraaion  pattern  which  arise  from  apertures  within  the  beam 
propagation  system  with  very  little  transmission  loss. 

The  paoe*’  bv  Thompson  and  Krisl'’8  oresents  a  comoaa  overview  of  the  subject  of 
apodisation  as  applied  to  coherent  imaging  and  beam  propagation  systems  The 
problem  of  reducing  the  oscillations  of  the  near  field  pattern  in  an  aberration-free 
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system  is  solved.  A  comparison  of  five  apodisers  relative  to  several  far-field  performance 
parameters  is  also  presented. 

2.3.2  Studies  of  Systems  with  Aberrations 

Most  studies  of  systems  with  aberrations  have  assumed  incoherent  illumination, 
while  only  a  few  have  considered  the  use  of  partially-coherent  illumination.  There  have 
apparently  been  no  studies  specifically  concerning  the  use  of  apodisation  to  improve  the 
performance  of  aberrated  coherent  imaging  systems.  Some  work  has  been  done  by 
considering  the  Gaussian  nature  of  laser  beams  m  aberrated  beam  propagation  systems. 

2.3. 2.1  Incoherent  Optical  Systems 

The  majority  of  the  work  aimed  at  using  apodisation  to  improve  the  performance 
of  aberrated  imaging  systems  has  occurred  in  the  last  few  years.  However,  Tsujiuchi'sS^ 
review  in  1963,  discussed  the  derivation  of  optimum  amplitude  and  phase  filters  to 
partially  correct  for  defocusing  and  spherical  aberration  in  systems  with  circular  pupils. 

A  group  at  the  University  of  Laval,  has  been  studying  the  effect  of  apodisers  m  the 
presence  o^  aberration  Biswas  and  Boivm^o  investigated  ootimum  apodisers  wh'cn  we'^e 
optamed  by  maximizing  the  encircled  energy  in  circular,  aberration-free  systems.  Then, 
*nev  '.neoretically  nvestigated  the  oerformance  of  these  apodisers  :n  the  oresence  of 
primary  and  secondary  spherical  aberrations.  Analytic  expressions  for  various  important 
properties  af  the  diffracted  field  m  the  presence  of  apodisers  and  spherical  aberration 
were  presented,  including  amplitude  point  spread  function,  irradiance  distribution, 
fractional  encircled  energy,  Strehl  ratio,  and  the  two-pomt  resolution.  The  authors 
•o.-d  *hat  ’he  ^se  of  these  filters  in  circular  aoe^Cated  systems  aiso  .ncreaseo  tne 
encircled  energy.  These  same  authors  also  investigated  the  influence  of  primary 
3s*:gmat:sm,  on  the  performance  of  these  optimum  aoodisei-si? 


Attacking  the  problem  directly,  Hazra  et  al^o  found  optimum  apodisers  which 
maximize  encircled  energy  in  the  presence  of  spherical  aberration  and  defocus.  To  make 
the  filter  manufacturing  process  easier,  they  designed  optimum  ladder-step  apodisers 
having  a  finite  number  of  uniformly  transmitting  zones,  each  with  a  different 
transmittance.  They  illustra;  ..  ..•:e  method  with  an  example  of  an  optimum  8-step 
apodiser  used  with  different  amounts  of  defocus,  primary,  and  secondary  spherical 
aberrations.  They  calculated  the  image  irradiance  distribution,  encircled  energy,  and 
Strehl  ratio.  They  predict  that  the  performance  of  the  system  as  measured  by  the 
encircled  energy  will  be  improved. 

The  problem  was  approached  in  a  slightly  different  manner  by  Yzuel  and  Calvo^V 
They  optimized  the  Strehl  ratio  by  varying  the  amplitude  and  phase  of  filters,  in  the 
circular  pupil  of  an  aberrated  optical  system.  The  optimization  process  was  halted  when 
the  Strehl  ratio  exceeded  0.8.  The  point  spread  function,  MTF,  and  the  irradiance 
distribution  along  the  axis  were  then  computed.  This  method  was  shown,  theoretically, 
to  improve  both  the  Strehl  ratio  and  the  MTF.  Mints  and  Prilepskii62  also  chose  to 
maximize  the  Strehl  ratio  using  phase  and  amplitude  filters.  They  considered  a  system 
raving  third-order  sonencal  aberration  ana  found  an  apodlsation  which  theoretically 
improved  the  Strehl  ratio. 

A  Polisn  group  searched  for  apodisers  which  minimized  the  second  moment  of  the 
r’-adiance  distribution  Magiera  et  ai63  did  this,  n  the  presence  of  four  orders  of 
soherical  aberration  by  varying  the  coefficients  of  a  polynomial  describing  the 
transmittance  of  the  apodiser  This  group,  m  a  seoarate  papers^,  then  examined  the 
effects  of  their  optimum  filter  on  the  second  moment  of  the  imaoe  irradiance 
distribution  in  the  presence  of  three  orders  of  spherical  aberrations.  They  showed 


numerical  results  indicating  the  predicted  improvement  in  performance.  They  also 
compared  their  filter  to  an  optimum  Bessel  filter  which  was  derived  by  Asakura  and 
Ueno65  for  the  aberration-free  case. 

2.3. 2.2  ?■;  •  Coherent  Optical  Systems 

Using  the  polynomiel  filter  previously  mentioned,  Magiera  and  Pluta^^  derived 
formulae  for  the  amplitude  impulse  response,  second  moment  of  the  point  spread 
function,  the  two-point  irradiance  distribution  in  partially  coherent  light,  the  point 
image  contrast,  and  the  equations  which  must  be  satisfied  by  the  Sparrow  and  Rayleigh 
resolution  criteria  in  the  presence  of  third,  fifth,  and  seventh  order  spherical  aberrations. 

Gupta  and  SmghSt  have  investigated  the  influence  of  primary  astigmatism  on  the 
irradiance  distribution  in  the  Fraunhofer  diffraction  patterns  formed  by  an  optical  system 
with  a  circular  aperture  under  partially  coherent  illumination  when  the  mutual 
coherence  function  contains  a  spatially  non-stationary  quadratic  phase  term.  Results  are 
presented  which  illustrate  the  degrading  effects  of  various  amounts  of  astigmatism  and 
the  extent  of  improvement  obtained  when  using  a  apodising  filter 

2. 3.2. 3  Coherent  Optical  Systems 

AUhougn  there  have  oeen  no  stuaies  specifically  addressing  apooisation  m 
aberrated  coherent  imaging  systems,  portions  of  several  of  the  studies  previously 
outlined  have  some  relevance  to  the  proposed  research,  in  several  cases-'. 58, 66_  the 
■amplitude  impulse  response  was  calculated.  However,  these  studies  did  not  address  the 
imaging  properties  of  the  aberrated  and  apodised  systems.  Magiera  and  Pluta^^  derived 
a  formula  for  the  two-point  irradiance  distribution  in  partially-coherent  light.  Taking  the 
fully  coherent  limit,  the  formula  presented  will  provide  a  useful  comparison  for  one  of 
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the  results  of  the  Chapter  IV.  Barakat^s  has  calculated  and  plotted  the  images  of  edges, 
discs,  and  bars  in  the  presence  of  defocus,  coma,  and  spherical  aberration.  In  aberrated 
beam  propagation  systems,  the  effect  of  a  truncated  Gaussian  apodisation  has  been 
investigated  by  several  researchers.  Holmes  et  al69  have  considered  the  effect  of 
Gaussian  as  well  as  sinusoidal  phase  aberrations  on  Gaussian  beams  in  the  near  and  far 
fields  of  an  annular  system.  LowenthaPO.Ti  has  examined  the  far-field  diffraction 
patterns  for  Gaussian  beams  with  spherical  aberrations;  and  has  developed  an 
aberration  tolerance  criterion  for  Gaussian  beams  For  these  systems  the  apodisation  was 
due  to  the  inherent  nature  of  the  laser  beam 

2.3.3  Previous  Research  Conclusions 

An  extensive  amount  of  work  has  been  done  using  apcxlisation  to  improve  the 
performance  of  aberrated  optical  systems.  Some  euthon  have  used  filters  obtained  m 
aberration-free  studies  while  others  have  derived  optimal  filteis  assuming  known 
amounts  of  fixed  aberration.  Generally,  Such  apodisers  have  improved  performance.  Out 
in  a  few  cases  they  have  actually  degraded  performance 

~o  summarize  tne  puDitsh.eo  research,  most  of  tne  'esearc'cs  “a-e  st^c  ec 
the  effects  of  apodisation  on  aberrated  systems  have  assumed  the  use  of  incoherent 
.ilumination.  Those  who  studied  coherent  systems  nave  mostly  concentratec  on  the 
image  of  a  single  point  source.  They  have  been  interested  in  measures  of  performance 
Such  as  the  irradiance  point  spread  function,  encircled  energy,  and  Strehl  ratio  Only  a 
few  have  considered  the  image  of  two  pomt  objects  m  partially  coherent  light.  The  issue 
of  generalized  imaging  in  apodised  and  aberrated  coherent  imaging  systems  has  not 


been  addressed 


2.4  THE  GAUSSIAN  APODISER 


The  Gaussian  apodiser  is  the  only  type  of  apodiser  considered  in  this  research. 
Only  one  type  of  apodiser  was  used  because  most  of  the  apodisers  which  were  found  to 
be  optimum  for  various  other  applications  differ  little  from  each  other  or  from  a 
Gaussian.  Additionally,  the  Gaussian  was  chosen  because  it  has  some  intuitively  pleasing 
effects  when  used  in  a  coherent  optical  imaging  system. 

The  first  consideration  when  chosing  an  apodiser  is  that  it  produce  an  amplitude 
impulse  response  which  is  real  and  positive's.  Many  of  the  deleterious  effects  seen  m 
coherent  imaging  (edge  ringing  for  example)  occur  because  the  amplitude  impulse 
response  has  negative  regions.  By  contrast,  an  incoherent  system  has  an  impulse 
response  which  is  everywhere  real  and  positive.  This  realization  leads  to  the  conclusion 
that  it  would  be  appropriate  to  make  the  amplitude  impulse  response  for  the  coherent 
case  also  real  and  positive.  The  amplitude  impulse  response  of  a  coherent  imaging 
system  is  proportional  to  the  Fourier  transform  of  the  exit  pupil  The  Fourier  transform 
of  a  Gaussian  function  is  itself  a  Gaussian  function  So  a  Gaussian  apodiser  m  the  exit 
pupil  will  produce  an  amplitude  impulse  response  which  is  real  and  oositive,  ‘'his  .s  not 
strictly  true  because  the  Gaussian  will  always  be  truncated  by  the  finite  aperture  of  the 
exit  pupii.  .However,  if  the  Gaussian  is  truncated  far  from  its  center,  this  statement  wiii  be 
approximately  true. 

Some  of  the  apodisation  functions  suggesteo  for  use  m  conerent  optical  systems 
are  shown  in  Fig.  2.6,  where  the  amplitude  transmittance  is  plotted  versus  the  radius  of 
the  aperture.  The  radius  of  each  function  has  been  scaled  so  that  the  truncation  point 
corresponds  to  a,  the  radius  of  the  aperture  being  apodised.  With  the  exception  of  the 
Lukosz  limit,  the  funaions  are  all  very  similar. 


An  important  consideration  when  a  real  and  positive  amplitude  impulse  response 
is  desired  was  presented  in  a  strict  mathematical  sense  by  Boas  and  Kac^^  and  later 
independently  from  a  signal-processing  point  of  view  by  Lukosz^s.  They  showed  that 
there  is  an  upper  limit  on  the  value  of  a  space-limited  funaion  if  the  funaion  is  to  have 
an  everywhere  positive  Fourier  transform.  This  limit  is  shown  in  Fig.  2.6. 

The  Gaussian  was  chosen  for  use  in  this  research  because,  as  well  as  producing  a 
real  and  positive  amplitude  impulse  response  and  not  differing  significantly  from  many 
of  the  apodisers  which  have  been  found  useful,  it  has  some  nice  effects  on  the  system. 
One  effect  is  that  the  transform  of  a  Gaussian  is  a  Gaussian.  Also,  when  a  Gaussian  is 
used,  the  imaging  system  has  some  similarities  to  a  laser  propagation  system.  Both 
effects  aid  the  intuitive  process  when  thinking  about  the  effects  of  the  apodiser  on  the 
system 


CHAPTER  3  THE  AMPUTUDE  IMPULSE  RESPONSE  •  THEORY 

3.1  INTRODUCTION 

The  amplitude  impulse  response,  as  pointed  out  in  the  last  chapter,  is  central  to 
the  analysis  of  optical  imaging  systems.  In  such  systems,  the  image  amplitude  can  be 
found  by  performing  a  convolution  of  the  geometric  image  with  the  amplitude  impulse 
response. 

The  amplitude  impulse  response  has  been  calculated  for  a  simple  system  having 
various  combinations  of  aberrations  and  apodisations.  The  physical  situation  is  depicted 
in  Fig.  3.1.  A  unit-amplitude  plane-polarized  coherent  plane  wavefront  is  incident  on  a 
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Fig. 3. 1  The  di  ffraction  geometry  used  in  calculated  the  amplitude,  impulse  response. 

test  lens  of  radius  a  and  focal  length  f .  This  is  equivalent  to  an  imaging  system  with  a 
point  object  at  infinity.  The  Vo  and  v,  axes  which  are  perpendicular  to  both  the  and  x. 
axes  respectively  and  to  '.he  optical  axis  are  net  shown  "he  'ens  n  general  "as 

aberrations  and  an  apodiser  could  be  adjacent  to  the  lens.  The  (io,Vo)  plane  is  chosen  to 


optical  axis 


coherent 

illumination 


be  coincident  with  the  lens. 


In  this  case,  (2.11)  is  used  directly  as 


The  field  amplitude  within  the  exit  pupil  is  described  by  (2.7).  Combining  (2.7)  and  (2.1 1) 
yields 


where  the  notation  f'{  }  refers  to  the  Fourier  transform  of  the  quantity  within  the 
brackets,  the  spatial  frequencies  of  the  transform  are  f,  =  x,  Af  and  fysy,  Af,  and  the 
terms  m  front  of  the  integral  have  been  dropped.  The  term  U,„(io,>'o)  (2-7),  is  the 

amplitude  of  the  input  wavefront,  which  in  this  case  equals  one. 

The  aberration  term  ♦  Cxo^o)  ‘f'  The  exponent  of  (3.2)  can  be  expressed  in  terms  of 
the  Seidel  coefficients  by  describing  the  aberrated  wavefront  with  the  expression 

-P^p.G  j  ^  ^  ’  (2  32) 

n  =0 ( *  0 

where  the  coefficients  a^i  and  b„i  have  units  of  wavelength  and  p  has  been  normalized  to 
unity  at  the  edge  of  the  aperture.  This  is  the  form  of  the  aberration  function  assumed 
when  referring  to  the  amount  of  aberration  present  in  a  particular  system. 

When  calculations  are  performed  however,  the  monomial  form  of  the  Zernme 
polynomial 

7*0  J  =  0 

iS  used  instead.  This  is  because  the  Fourier  transTorm  operation  .noicateo  :n  (3.2)  s 
performed  on  the  computer  as  a  discrete  fast  Fourier  transform  on  a  regular  rectangular 
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The  apodisation  used  throughout  this  dissertation  is  in  terms  of  the  real  amplitude 
transmittance.  It  has  the  form 


(3,3b) 

where  the  constant  in  the  exponent  was  chosen  so  that  the  value  of  A(r)  at  the  edge  of 
the  aperture,  r=  1,  was  equal  to  0.050  (see  Fig  3.2).  Thus  there  was  very  little  amplitude 
m  the  transmitted  wavefront  at  the  edge  of  the  aperture.  This  should  allow  the 
separation  of  the  effects  of  the  aberrations  from  the  effects  of  the  hard  aperture 

A  VAX-1 1/750  computer  was  used  to  compute  (3  2)  using  a  program  written  in 
FORTRAN-77.  The  Fourier  transform  was  accomplished  by  using  a  fast-Fourier 
transform73  subroutine  from  the  IMSL^*  package.  This  program,  labeled  PSF.FOR,  is 
documented  in  Appendix  1 . 

The  quantity  U|(i,.y,)  is.  in  general,  complex  Consequently,  the  output  of  the 
orogram  is  m  terms  of  real  and  imaginary  coefficients  That  is. 
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(3.4) 


where  a(xi,y0  af^d  b(x„yj  are  the  just  mentioned  coefficients  and  m(Xi,yi)  and  p(x„y^)  are 
•.he  modulus  and  onase,  ,'espectiveiy,  of  U(i.,y..). 

Coherently  illuminated  systems  are  linear  in  complex  amplitude.  However,  if 
measurements  are  to  be  made  m  an  expenm.ent.  it  is  the  irradiance  I,(x.,y.)  which  is 


usually  measured  .  The  irradiancc  is  related  to  the  complex  amplitude  by 

I  -  <3.5) 

where  the  asterisk  denotes  the  complex  conjugate. 

3.2  COMPARISONS  WITH  KNOWN  RESULTS 

A  typical  output  of  this  program  is  displayed  in  Fig.  3.3.  The  system,  in  this  case 
was  unaberrated  and  unapodised  and  the  pupil  function  was  a  circular  aperture.  The 
distance  coordinates  used  in  Fig.  3.3  ,  as  well  as  in  many  other  figures  in  the  thesis,  are 
the  canonical  distance  coordinates  u  and  u.  They  are  canonical  in  the  sense  that  the 
pattern  displayed  in  these  figures  is  independent  of  the  exit  pupil  size,  the  wavelength  of 
illumination,  and  the  distance  from  the  exit  pupil  to  the  image  plane.  Specifically,  these 

new  coordinates  are  defined  by 

2na  2na 

u  —  - X.  and  u  =  - y  ,  ,, 

Ad  ‘  Ad  (3  6) 

t  { 

where  a  is  the  radius  of  the  exit  pupil,  X  is  the  -vavelength  of  the  illumination,  d,  is  the 
distance  from  the  exit  pupil  to  the  image  plane,  and  x,  and  y,  are  the  real  coordinates  in 
the  image  plane. 

A  coherent  optical  system  is  linear  in  field  amplitude,  which  is  a  complex  quantity 
5o  unlike  an  incoherent  system,  the  phase  in  the  optical  field  is  critically  imponant  m 
determining  the  final  image  irradiance  distribution.  For  this  reason,  the  phase  has  been 
calculated  for  the  unaberrated  case  m  Fig.  3.3  and  for  many  of  the  aberratea  cases 
considered  later  in  this  chapter. 

The  phase  distribution  in  the  amplitude  impulse  response  of  an  unaberrated 
system  (see  Fig.  3.3b)  Is  uniform  everywhere  except  along  concentric  circles  where  the 
phase  jumps  discontinuously  by  an  amount  equal  to  it  radians.  These  jumps  occur  at  the 
same  spatial  locations  as  do  the  zero  values  in  the  modulus  distribution  of  Fig.  3.3a.  This 
implies  that  every  other  ring  in  the  amplitude  impulse  response  is  composed  of  negative 
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values.  It  is  these  negative  regions  which  cause  the  ringing  phenomena  seen  in  the 
Image  of  an  edge  (see  Pig.  1.2).  Thus  the  knowledge  of  the  phase  distribution  in  the 
amplitude  impulse  response  is  very  important  in  determining  the  behavior  of  the  system 

The  irradiance  cist'’'!:  this  impulse  response  is  shown  in  Fig.  3.3c.  It  is 

obtained  by  squaring  the  ” : -  distribution.  Since  the  optical  system  in  this  case  is 
unaberrated  and  has  a  circular  exit  pupil,  it  is  expected  that  the  irradiance  distribution 
would  approximate  an  Airy^  pattern  A  careful  comparison  of  Fig.  3.3c  with  a  theoretical 
Airy  distribution  reveals  the  error  between  the  two  patterns  to  be  less  than  1%.  This 
error  arises  mainly  from  the  problem  of  inadequately  representing  a  circular  aperture 
with  a  rectangular  array  of  samples.  Nevertheless,  this  agreement  is  quite  close  and  thus 
lends  credibility  to  the  theoretical  development  and  the  computer  program  used  to 
realize  that  development. 

The  program  gains  additional  credibility  when  the  irradiance  patterns  in  the 
presence  of  simple  aberrations  are  compared  to  similar  results  in  the  literature.  Fig  3  4 
displays  the  theoretical  irradiance  impulse  responses  for  two  optical  systems  with 
different  amounts  of  aberrations  The  impulse  response  in  Fig.  3  4a  Is  from  an  optical 
system  having  0.48  waves  of  coma;  i.e..  $  (  p ,  9 )  =  O.48Ap3cos0  The  Quantity  X  refers  to 
the  wavelength  of  the  coherent  illumination  '''he  '■ing  structure  has  been  distorted  by 
the  presence  of  the  aberration  such  that  the  ring  has  been  flattened  on  one  side  of  the 
central  lobe  and  Heightened  on  the  oooosite  side.  In  i^ig.  3.4b  the  aberration  m  the 
system  is  0.16  waves  of  astigmatism  ($  (  p ,  9 )  =  0  1  SXp^cos^G  ) .  These  particular 
aoerration  values  were  chosen  because  another  author  (Nijboer  75,76)  has  calculated 
tirr'lar  'esults  using  these  values.  Niiboer's  results  are  -eoroduced  'n  Fig.  3.5  r/here  the 
plots  show  the  irradiance  distribution  for  (a)  0.48  wave  of  coma  and  (b)  0.16  waves  of 
astigmatism. 
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Fig  3  5  Irradiance  contours  in  the  images  of  point  sources  through  an  ootical  system 
having  ^aj  0.48  waves  of  coma  ana  (b<  0. 16  waves  of  astigmatism.  These  plots  are 

^rom  Niiboer^S.I’S 
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These  plots  display  the  irradiance  impulse  responses  in  terms  of  the  contours  of  equal 
irradiance.  The  presentation  formats  for  Figs.  3.4  and  3.5  differ  but  the  qualitative 
agreement  is  good.  The  comatic  patterns  both  show  a  disappearance  of  the  first  ring  on 
one  side  of  the  central  lobe  and  an  enhancement  of  that  nn-  -  opposite  side;  and 
the  astigmatic  patterns  both  show  a  distinct  four-lobed  character  in  the  first  ring.  As  a 
quantitative  measure,  the  Strehl  ratio  was  computed  for  the  two  cases  and  compared  to 
Nijboer's  results.  The  Strehl  ratio  for  each  plot  of  the  two  figures  is  shown  beside  the 
plot.  The  agreement  here  is  very  good. 

Similar  comparisons  were  made  for  other  amounts  of  aberrations.  In  each  case  the 
qualitative  agreement  between  the  model  developed  here  and  the  data  published  by 
other  authors  was  quite  good.  The  quantitative  agreement  was  typically  on  the  order  of 

1%. 

3.3  THE  EFFEa  OF  ABERRATIONS  ON  THE  IMPULSE  RESPONSE 


3.3,1  Defocus 

The  aberration  of  defocus  has  the  functional  form  $(p,9)  =a2iD^,  where  the 
coefficient  a2i  is  the  amount  of  aberration.  The  subscripted  numbers  2  and  1  come  from 
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3  6  The  amplitude  imoulse  response  (modulus  and  phase)  in  the  presence  of  0.5A 
defocus  and  for  the  case  of  an  unapodised  ana  Gaussian  apooised  aperture.  The 
top  two  plots  arc  for  the  unapodised  case  while  the  bottom  two  are  for  the  case 
of  a  Gaussian  apodiser.  The  vertical  scales  for  the  modulus  plots  (left  hand 
column)  and  the  phase  plots  (right  hand  column)  are  indicated  by  the  top  two 
plots.  This  same  scaling  is  used  m  .-igs  3.9,3.  f  2.  ana  3  '.5. 


The  modulus  and  phase  of  the  unapodised  amplitude  impulse  response  (top  two 
plots)  should  be  compared  to  the  analogous  plots  of  Fig.  3.3  where  there  are  no 
aberrations  in  the  system.  The  peak  value  of  the  modulus  in  the  aberrated  case  has 
decreased  relative  to  the  unaberrated  case.  The  zero  values  in  the  modulus  pattern  for 
the  unaberrated  case  have  evolved  to  relative  minimums  which  do  not  go  to  zero.  The 
phase  of  the  aberrated  amplitude  impulse  response  (upper  right  plot  of  Fig.  3.6)  no 
longer  has  the  discontinuities  evident  in  the  unaberrated  case. 

When  the  apodiser  described  by  (3.3)  and  plotted  in  Fig.  3.2  is  applied  to  this 
aberrated  system,  the  modulus  of  the  amplitude  impulse  response  (lower  left  plot  of  Fig. 
3  6)  IS  considerably  smoothed,  as  is  the  phase.  The  phase  varies  by  less  than  it  radians 
over  the  region  of  the  modulus  plot  where  the  modulus  is  greater  than  10%  of  its  peak 
value.  So  this  amplitude  impulse  response  does  not  change  sign  until  the  absolute  value 
of  the  amplitude  is  quite  small.  Thus  the  impulse  response  is  almost  real  and  positive. 

This  has  important  implications  for  the  imaging  performance  of  this  optical 
system.  For  instance,  the  ringing  in  the  coherent  image  of  an  edge  is  caused  by  the 
negative  regions  of  the  impulse  response.  In  this  case  the  apodiser  has  smoothed  the 
amoiitude  impulse  response  such  that  it  has  very  little  amplitude  in  the  regions  where 
thei'e  are  negative  values  of  amoiitude.  It  can  be  expected  then,  that  the  image  of  an 
edge  througn  this  system  would  be  free  from  ringing.  The  results  in  the  next  chapter 
'OP'  rrrt  -ms  exoec’.ation 

The  amplitude  impulse  response,  both  unapodised  and  apodised,  for  other  values 
of  defocus  are  shown  in  Fig  3.7  Here  the  amount  of  aberration  is  different  for  each 
■plump  of  plots,  .ars'lng  ^rom  0.'*  X  on  the  left  to  1.0  X  on  the  right.  From  •‘■'ese  plots  we 
can  see  the  evolution  of  the  modulus  and  phase  as  more  defocus  is  added  to  the  system. 


The  modulus  and  phase  along  slices  through  the  center  of  some  of  these  impulse 
responses  are  shown  in  Fig.  3.8.  The  relationship  of  phase  to  modulus  is  clearly  seen  in 


this  figure. 

Apodisation,  in  each  case  of  the  last  three  figures,  smoothes  both  the  modulus  and 
phase.  In  each  case,  the  amplitude  impulse  response  becomes  almost  real  and  positive 
when  the  apodiser  is  applied.  There  are,  however,  limits  to  this  process.  As  the  amount 
of  aberration  increases  the  apodiser  becomes  less  effective  in  making  the  amplitude 
impulse  response  almost  real  and  positive.  For  the  case  of  one  wave  of  defocus,  it 
appears  that  the  phase  changes  by  more  than  it  radians  over  the  region  where  the 
modulus  is  still  relatively  large. 

3. 3. 2  Spherical  Aberration 

The  aberration  of  spherical  has  the  functional  form  ♦  ( p ,  0 )  =  a42p4.  Spherical  is  a 
radially  symmetric  aberration.  Due  to  its  fourth  power  dependence  on  the  radial 
distance  parameter  p,  spherical  aberration  describes  a  wavefront  having  the  largest 
deviation  from  the  spherical  reference  wavefront  of  any  of  the  aberrations  considered 
here  A  calculation  of  the  impulse  response  (modulus  and  phase  when  aa2  =  0  5  X  is 
shown  in  Fig.  3.9.  The  top  two  plots  (the  unapodised  case)  should  be  compared  to  the 
'jpabe'’rated  'moulse  resoonse  of  Fig,  3  3  ~he  presence  of  sohencal  aberration  causes  a 
decrease  m  the  value  of  the  central  peak  (Strehl  ratio)  and  an  increase  in  the  energy  m 
tne  side  lobes,  particularly  the  first  side  iobe.  The  position  of  the  ring  of  minimum  values 
between  the  central  lobe  and  the  first  side  lobe  remains  unchanged. 

The  impulse  response  when  an  apodiser  is  used  is  shown  in  the  bottom  two  plots 
of  =ig.  3.9.  The  aoodiser  used  is  the  one  described  bv  13.3)  and  plotted  in  ^ig.  3.2.  As  m 
the  case  of  defocus,  the  use  of  the  apodiser  has  resulted  in  a  much  smoother  impulse 
response  An  e.rammation  of  the  phase  distcbution  shows  that  the  phase  is  near'v 
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The  amplitude  impulse  response  (modulus  and  phase)  with  varying  amounts 
of  spherical  aberration  for  the  case  of  an  unapodised  and  Gaussian  apodised 
exit  pupil.  The  amount  of  aberration  for  each  column  is  irsdicated  at  the 
bottom  of  that  column. 


uniform  over  the  region  of  the  impulse  response  having  significant  amounts  of  energy. 
This  impulse  response  can  also  be  described  as  being  almost  real  and  positive. 

The  evolution  of  the  unapodised  and  apodised  amplitude  impulse  responses  as 
more  spherical  aberration  is  added  to  the  system  is  shown  in  Fig.  3.10.  Central  slices 
through  some  of  these  plots  are  shown  in  Fig.  3.11.  The  same  general  phenomena  seen 
in  the  case  of  defocus  are  seen  here  as  well.  The  use  of  the  apodiser  results  in  an  impulse 
response  which  is  free  from  side  lobes  in  the  modulus  pattern  and  which  has  a  relatively 
flat  phase  over  the  region  where  there  is  a  significant  amount  of  energy. 

Again  there  are  limits  to  this  process.  When  the  amount  of  spherical  aberration  is 
about  one  wave,  the  impulse  response  has  significant  amounts  of  energy  in  regions 
where  the  phase  has  changed  by  X/2.  So  the  apodiser  is  not  totally  effective,  although 
the  apodised  impulse  response  is  still  much  smoother  than  the  unapodised  one. 

3.3.3  Coma 

The  aberration  of  coma  can  be  described  by  ♦(pjS)  sa3ip*cos0.  Coma  is  the 
first  aberration  considered  for  which  the  wavefront  in  the  exit  pupil  depends  on  the 
polar  angle  9  as  well  as  the  radial  distance  p.  This  aberration  therefore  produces  an 
unsymmetncal  amplitude  impulse  response  as  seen  m  Fig.  3.12.  The  use  of  the  apodiser 
n  this  case  appears  to  be  less  effective  than  m  the  orevious  cases,  because  the  first  side 
ioDe  IS  still  evident  in  the  modulus  of  the  apodised  impulse  response  (lower  left  plot  of 
“ig.  3.1 2).  Also  .  there  is  a  tr  ohase  change  m  tne  'egion  of  *his  first  side  lobe. 

The  evolution  of  these  impulse  responses  with  increasing  amounts  of  aberration  is 
snown  m  Fig.  3.13.  Central  slices  of  some  of  these  data  are  shown  in  rig.  3.14.  The  slices 
are  along  the  axis  showing  the  minimum  amount  of  svmmetrv  fhe  same  oene'’al 
conclusions  that  were  drawn  for  the  cases  of  defocus  and  spherical  aberration  can  be 
drawn  for  this  case  as  well.  First,  the  aoodiser 's  e^fec*ive  m  transformino  *he  aberf'ated 
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The  amplitude  impulse  response  (modulus  and  phase)  in  the  presence  of  -O.SX 
y-axi$  coma  and  for  the  case  of  an  unapodised  and  Gaussian  apodised 
aperture.  The  top  two  olou  are  for  the  unaoodised  case  while  the  bottom 
two  are  for  the  case  of  a  Gaussian  apodiser.  See  Fig.  3.6  for  scaling. 
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impulse  response  into  a  much  smoother  function.  Second,  there  are  limits  to  the 
effectiveness  of  the  apodiser.  The  limit  in  this  case  appears  at  a  lower  value  of  aberration 
than  in  the  other  cases.  For  as  little  as  half  a  wavelength  there  is  a  ii  phase  change  in  a 
region  of  the  impulse  response  where  *  -  -  -  dulus  shows  a  side  lobe. 

3.3.4  Astigmatism 

The  aberration  of  astigmatism  is  described  by  ♦  (  p  ,  6  )  =  a2oP^cos20.  Like  coma, 
astigmatism  is  an  unsymmetric  aberration.  The  appearance  of  the  unapodised  and 
apodised  impulse  responses  when  a2o  =  0.5X  is  shown  in  Fig.  3.15.  Both  phase  plots 
shoNv  a  saddle  shace  which  is  characteristic  of  astigmatism.  Astigmatism  results  in  a 
wavefront  which  has  different  radii  of  curvature  along  two  orthogonal  directions  in  the 
plane  of  the  exit  pupil.  This  behavior  is  evident  also  in  the  focal  plane  as  seen  in  Fig. 
3.15. 

The  effect  of  apodisation  is  to  smooth  both  the  amplitude  and  phase  of  the 
impulse  response.  This  behavior  holds  as  the  amount  of  aberration  is  increased  from  0.1 
X  to  1.0  X,  as  seen  in  Figs.  3.16  and  3.17.  The  general  shapes  of  functions  in  Fig.  3  17 
differ  little  from  the  analogous  plot  for  the  case  of  defocus  (Fig.  3.8).  This  is  because,  m 
one  dimension,  astigmatism  results  in  a  wavefront  wtiich  is  spherical  but  having  a  radius 
of  curvature  different  from  the  reference  wavefront. 

The  same  conclusions  as  before  can  be  drawn.  That  is,  the  apodiser  smooths  the 
aberrated  impulse  response  for  a  limited  amount  of  aberration.  The  limit  in  this  case 
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appears  to  be  about  one  wavelength. 
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0'  ASTIGMATISM 


The  amplitude  impulse  response  (modulus  and  phase)  with  varying  amounts 
of  0°  astigmatism  for  the  case  of  an  unapodised  and  Gaussian  apodised  exit 
pupil.  The  amount  of  aberration  for  each  column  is  indicated  at  the  bottom 
ot  that  column. 
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3.3.5  On>axis  Calculations 

A  common  feature  in  all  of  the  impulse  data  is  that  the  addition  of  aberrations  to 
an  unaberrated  optical  system  always  results  in  a  decrease  in  the  peak  value  of  the 
central  lobe  of  the  modulus.  This  behavior  is  plotted  in  Fig.  3.18  for  various  values  of  the 
four  aberrations  and  the  two  cases  of  apodisation.  The  peak  values  have  been 
normalized  to  unity  for  the  unaberrated  system.  From  this  plot,  it  can  be  seen  that  the 
peak  value  of  modulus  occurs  when  there  are  no  aberrations.  The  addition  of  any 
amount  of  the  four  aberrations  to  an  unaberrated  system  results  in  a  decrease  in  the 
peak  value  Notice  that  in  the  case  of  defocus,  the  central  peak  completely  disappears 
for  $  =  10  wave  This  is  a  so-called  fringe  of  defocus. 

Figure  3.19  shows  the  behavior  of  the  central  value  of  irradiance  as  the  same 
aberrations  were  added.  Fig.  3.19  can  be  obtained  by  squaring  each  value  in  Fig.  3.18. 
When  the  irradiance  is  used  as  is  Fig.  3.19,  the  plots  could  also  be  labeled  as  the  Strehl 
ratio.  Note  that  the  use  of  the  apodiser  immediately  decreases  the  Strehl  ratio  by  a 
factor  of  10  for  the  unaberrated  case. 

The  last  figure  in  this  chaDte^  F:g  3  20,  shows  the  behavior  of  the  phase  at  the 
central  peak  of  the  irr, pulse  resoonse  as  a  function  of  the  amount  of  aoerration  (i)  ana 
aoodisation  present.  For  most  of  t.he  abe'rations  tne  ohase  is  relatively  unchanged  as 
aberrations  are  added.  For  the  case  of  defocus,  however,  there  is  a  it  phase  jump  m  the 
unapodised  case.  The  location  of  this  ju.Tio  corresconas  to  the  location  of  the  zero  m  the 
modulus  of  the  unaberrated  impulse  response.  This  is  consistent  with  the  earlier 
observation  that  zeros  in  the  modulus  oistnoution  of  an  unaberrated  amplitude  impulse 
-esconse  are  iccomoomeo  bv  ••  ohase  ;u  res 
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3.4  Impulst  Rtsponst  Conclusions 


In  this  chapter  the  form  of  the  impulse  response  has  been  determined  for  various 
amounts  of  defocus,  spherical,  coma,  and  astigmatism  and  for  two  cases  of  apodisation. 
The  two  cases  were  no  apodisation  and  the  Gaussian  apodiser  described  by  (3.3).  Since 
coherent  optical  systems  are  linear  in  the  complex  field  amplitude,  the  impulse  response 
was  displayed  in  terms  of  both  its  modulus  and  phase. 

It  was  seen  that,  in  the  absence  of  the  apodiser,  the  addition  of  an  aberration 
significantly  modified  the  form  of  the  impulse  response.  However,  when  the  apodiser 
was  m  place,  the  form  of  th*e  impulse  response  did  not  vary  much  as  an  aberration  was 
increased;  the  modulus,  for  the  most  part,  smoothly  fell  from  its  peak  value  to  zero  and 
the  phase  was  relatively  flat  over  the  region  where  there  was  a  significant  amount  of 
modulus. 

The  apodiser  was  effective  in  removing  the  'feet'  of  the  aberrated  impulse 
response  within  certain  limits.  Those  limits  were  approximately  1.0  wave  for  defocus, 
spherical,  and  astigmatism  and  0.5  waves  for  coma.  The  implication  is  that  within  these 
limits  the  apodiser  will  improve  certain  aspects  of  the  performance  m  a  coherent  imaging 
system;  e.g.,  removing  the  edge  ringing.  Even  beyond  these  limits  the  performance 
should  be  improved:  for  example,  the  edge  ringing  should  be  greatly  diminished  if  not 
eliminated. 
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CHAPTER  4  IMAGING  SYSTEMS  ANALYSIS  •  THEORY 


The  effea  of  apodisation  on  the  performance  of  aberrated  coherent  optical 
systems  has  been  theoretically  investigated.  The  objects  considered  in  this  analysis  were 
an  edge,  a  slit,  and  two  closely  spaced  points;  and  the  performance  of  the  optical  system 
was  measured  against  criteria  which  were  SF>ecific  for  each  of  these  test  objeas.  The 
performance  predictions  detailed  below  are  based  on  the  theory  described  in  Chapter  2. 
Experimental  verification  of  key  parts  of  the  theory  are  contained  in  Chapter  5. 

The  optical  system  was  modeled  as  shown  in  Fig.  2.2  which,  for  convenience,  is 


reproduced  here  and  labeled  as  Fig.  4, 1 .  Under  the  conditions  detailed  in  Section  2. 1 .2, 
Xo  Xi 


-ig.4  1  The  geometry  of  a  centered  ooticai  lystem  wmcn  gives  rise  to  an  aberrated 
wavefront  W.  The  aberration  function  4>  is  in  terms  of  the  distance  QQ*  between 
W  .and  the  reference  sohere  S. 


the  optical  field  in  the  image  plane  -s  related  to  the  optical  field  in  the  object 

3iane  by  the  '2. IS); 
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(x,,,,)  =  Flk(/',,)F-' 


(2.28) 


wnere  F{  }  and  F-^{  }  represent,  respectively,  the  forward  and  inverse  Fourier  transforms 
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of  the  quantities  within  the  brackets  and  k(f.g)  is  the  scaled  exit  pupil  transmittance. 
Equation  (2.28)  will  be  the  basis  for  the  predicted  results  within  this  chapter. 

The  results  are  plotted  in  terms  of  both  the  image  irradiance  and  the  square  root  of 
image  irradiance.  As  stated  before,  the  square  root  of  irradiance  is  a  good  measure  of  an 
image  because  the  nonlinearity  of  the  detection  process  is  removed.  However  the  usual 
presentation  seen  in  the  literature  is  in  terms  of  the  irradiance.  Thus,  both  scales  are  used 
in  this  chapter. 


4.1  EDGE  OBJECTS 

(4.1) 


(4  2) 

where  6  is  the  Dirac  delta  function  and  /"and  g  are  the  spatial  frequency  coordinates  m 
the  plane  of  the  exit  pupil: 

^  n 

'  -  Jjr  3r.a  d-—  l4  3 

0  0 

’’he  coordinates  (E.q)  refer  to  the  space  coordinates  n  »he  exit  pupil  plane  ana  is  the 
distance  from  the  object  plane  to  the  entrance  pupil 

Combining  (4  i)  -  (4,3)  into  (2.28)  results  m  the  optical  field 

f  ^  1  /'  M  /  n  N  1 1 

U  i  ,y  ;  =  I  -  K  .n  ijoi  —  I  ~  -  '  -) 
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An  edge  in  the  object  (xo,yo)  plane  can  be  described  by 


ify^SO 


[0  otherwise  * 

where  the  edge  has  been  aligned  to  be  coincident  with  the  Xq  axis. 
The  inverse  Fourier  transform  of  (4.1)  is  given  by^^ 
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where  the  problem  has  been  reduced  to  one  dimension  because  of  the  8(/)  term  in  (4.2). 


The  term  k(n)  describes  the  exit  pupil  transmittance  in  a  manner  analogous  to  (2.7);  i.e., 


where  A{r\)  and  (|>(ri)  represent  the  amplitude  and  phase  transmittance,  respeaively,  of 
the  exit  pupil  along  the  line  ^  =  0  and  B(n)  represents  the  finite  extent  of  the  exit  pupil 
along  that  same  line. 

The  calculation  of  (4.4)  for  various  values  of  aberration  and  apodisation  was 
accomplished  using  the  University  of  Rochester  Computer  Center's  DEC  10  computer. 
The  program  controlling  these  calculations  was  written  in  FORTRAN  and  employed  a  FFT 
routine  from  the  IMSL^^  library.  A  documented  copy  of  the  program,  labeled  EDGE. FOR, 
IS  in  Appendix  1.  The  program  is  a  simple  modification  of  the  program  PSF.FOR  which 
was  used  in  Chapter  3.  The  output  of  the  program  is  in  terms  of  the  irradiance  of  the 
image  which  is  propoaional  to  the  product  of  (4.4)  times  its  complex  conjugate.  The 
square  root  of  the  irradiance  is  also  a  program  output. 

Typical  results  from  this  program  are  depiaed  in  Figs.  4.2  through  4,4,  These 
figures  display  both  the  irradiance  and  the  square  root  of  the  irradiance  of  a 
perpendicular  slice  through  the  'mage  of  an  edge  when  the  optical  system  has  various 
amounts  (0  0,  0,5,  and  1.0  waves)  of  aefocus,  y-axis  coma,  or  spherical  aberration 
Separate  plots  ot  the  effects  of  astigmatism  are  not  shown  oecause  m  one  dimension 
astigmatism  has  the  same  functional  form  as  defocus  The  top  plots  in  each  figure  are  for 
an  unappdised  optical  system,  while  the  bottom  plots  are  for  a  system  with  a  Gaussian 
apodiser  having  an  amplitude  transmittance  of 


la  F' 


In  each  plot,  the  ordinate  is  m  terms  of  the  relative  irradiance  (or  square  root  of 
irradiance)  and  the  abscissa  is  in  terms  of  the  normalized  distance  v.  For  reference,  the 
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^igA  3  The  coherent  image  of  an  edge  through  an  ootical  svstem  having  0.0.  0  5  and  ’  0 
waves  01  y-axi$  coma.  The  top  plots  are  for  an  unapodised  optical  system  whne 
the  bottom  plots  are  for  a  system  having  a  Gaussian  apodiser  The  ordinates  of 
the  two  plots  on  the  left  are  in  terms  of  image  irradiance  (I)  while  the  ordinates 
for  the  others  are  in  terms  of  the  souare  root  of  image  irradiancefViV  Legend- 
_ jnaoerrateo, _ J  5  waves, _  _ i  .Owave 
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Fig. 4  4  The  coherent  image  of  an  edge  through  an  optical  system  having  0.0,  0.5,  and  i  0 
waves  oi  spnericai.  The  top  piots  are  lor  an  unapodised  optical  system  wniie  tne 
bottom  plots  are  for  a  system  having  a  Gaussian  apodiser.  The  ordinates  of  the 
two  plots  on  the  left  are  in  terms  of  image  irradiance  (I)  while  the  ordinates  for 
the  others  are  in  terms  of  the  square  root  of  image  irradiance  (Vl).  Legend; 
_ jnaoerratea, _ O.Swaves, _  _  _ l.Owave. 
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position  of  the  geometrically  imaged  edge  is  included  on  each  plot  as  a  line  of  short 
dashes. 

For  the  apodiser  used  here,  the  amplitude  transmittance  has  a  value  of  0.05  at  a 
distance  from  the  center  of  ri  =  1.0.  The  hard  aperture  also  has  a  radius  of  q  a  1.0. 
Because  the  Gaussian  apodiser  is  truncated  at  a  very  low  value,  the  discontinuity  in 
transmittance  at  the  edge  of  the  aperture  is  almost  removed.  A  plot  of  the  apodiser 
described  by  (4.6)  is  shown  in  Fig,  3.3. 

The  measures  of  importance  when  evaluating  coherent  edge  images  are  edge- 
ringing,  edge-shifting,  and  acutance. 

4.1.1  Edge  Ringing 

Edge  ringing  is  the  periodic  deviation  of  the  edge  image  from  the  geometrically 
ideal  edge  image.  This  phenomena  is  most  pronounced  in  the  top  plots  of  Fig.  4.2.  It  can 
be  quantified'*^  in  terms  of  the  ratio  of  the  maximum  positive  deviation  to  the  geometric 
exoected  value  of  unity.  This  ratio  for  various  values  of  the  three  aberrations  considered 
IS  snown  m  Fig,  4  5  (a-c).  From  these  plots  it  is  clear  that  the  use  of  apodisation  wiii 
reauce  tne  edge-rmging  effect,  especially  for  the  case  of  coma  shown  This  beravior  of 
an  aoodised  system  is  well  known  in  the  case  of  aberration-free  systems.  It  is  now 
aooarent  that  apodisation  will  reduce  eoge  ringing  even  m  the  presence  ot  aoerraiions. 
This  result  is  not  surprising  because  the  impulse  resoonse  of  the  system  has  been 
Significantly  smoothed,  as  seen  m  Chapter  3.  Essentially,  the  apodiser  acts  as  a  we'gnting 
'unction  in  the  exit  pupil  It  deemphasiies.  oy  attenuation,  the  regions  of  the  exit  oupi' 
near  the  edge  of  the  aperture,  it  is  in  these  same  regions  when  aberrations  are  present, 
that  the  phase  of  the  wavefront  deviates  most  from  the  ideal  shape  of  the  spherical 
reference  sphere.  The  effect  of  using  the  square  root  of  the  image  irradiance  is 


also  illustrated  in  this  figure.  As  expected,  the  amount  of  edge  ringing  is  reduced  in 
these  cases.  Although,  as  a  consequence  of  taking  the  square  root  of  the  irradiance,  the 
ringing  on  the  dark  side  of  the  edge  has  increased. 

4.1.2  Edge  Shift 

The  edge  shift  is  the  distance  from  the  theoretical  position  of  the  geometric  edge 
to  the  edge  defined  by  the  real  image.  The  location  of  the  imaged  edge  is  defined  here 
as  the  position  where  the  irradiance  is  half  the  value  of  the  first  peak.  The  dependence 
of  edge  shih  on  the  amount  and  type  of  aberration  as  well  as  the  presence  or  absence  of 
apodisation  is  illustrated  in  Fig.  4.6  (a  -  c).  As  seen  in  this  figure,  the  edge  shift  is  not  a 
strong  function  of  tne  presence  or  absence  of  the  apodiser.  In  general,  the  application  of 
the  apodiser  increases  the  amount  of  edge  shift.  However,  in  the  case  of  y-axis  coma,  the 
edge  shift  is  reduced  when  the  apodiser  is  used. 

The  amount  of  edge  shift  is  a  funrtion  of  how  the  edge  is  defined.  A  common 
choice  (one  not  made  here  however)  for  the  location  of  the  edge  is  the  point  where  the 
irraoiance  is  equal  to  0.5,  when  the  value  of  the  irradiance  far  from  the  edge  on  the 
Illuminated  side  is  equal  to  1.0.  If  this  choice  is  made  and  the  square  root  of  the 
irradiance  is  taken  then  there  will  be  no  edge  shift  for  the  aberrations  which  are 
symmetric  aoout  the  optical  axis,  e.g.,  astigmatism,  oefocus,  and  spherical.  The  presence 
of  coma  will  still  cause  the  edge  to  shift(see  Fig.  4.3).  The  problem  with  this  choice  is  that 
experimentally  it  is  often  difficult  if  not  impossible  to  determine  the  value  of  the  the 
brignt  side  of  the  edge  at  a  great  enougn  distance  f’-om  t^e  edge.  I’he  presence  of  noise 
and  other  edges  m  the  system  very  often  precludes  this 

It  IS  for  these  reasons  that  the  edge  is  defined  here  as  the  point  where  the 
irradiance  is  half  of  the  value  of  the  irradiance  at  the  maximum  of  the  peak  closest  to  the 
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Fig  4  6  The  amount  of  edge  shift  is  shown  for  various  amounts  of  the  aberrations  of  (a) 
sone'^icji,  {bl  yaxis  coma,  and  jc)  ae’ocjs  ”he  sond  curves  are  ‘or  ‘."e 

unapodised  cases  while  the  dashed  curves  are  for  the  apodised  cases.  The  edge 
shift  is  defined  as  the  distance  of  the  half  peak  irradiance  point  from  the 
geometric  edge.  The  plots  in  the  left  column  are  calculated  from  the  edge  image 
irradiance  while  the  others  are  from  the  square  root  of  the  edge  image 
irradiance. 


edge.  The  edge  is  similarly  defined  when  the  square  root  of  the  irradiance  is  considered. 
The  advantage  of  this  definition  of  the  edge  is  that  it  is  always  possible  to  find  the  first 
peak. 

4.1.3  EdgeAcutance 

Acutance  is  a  measure  of  the  the  slope  of  the  edge.  Specifically,  as  a  definition, 
acutance  is  defined  here  as  the  slope  of  the  irradiance  (or  square  root  of  irradiance)  at 
the  location  of  the  edge.  The  behavior  of  the  acutance  as  a  function  of  the  amount  and 
type  of  aberration  as  well  as  the  presence  or  absence  of  apodisation  is  shown  in  Fig  4,7  (a 
-  c)  In  each  case,  acutance  decreases  with  both  an  increase  of  the  amount  of  aberration 
and  the  application  of  apodisation. 

Using  the  square  root  of  the  irradiance  distribution  seems  to  have  some 
advantages  in  the  analysis  of  edge  objects.  Relative  to  the  use  of  the  irradiance 
distribution,  the  edge  ringing  and  shift  are  reduced  and  the  acutance  is  not  affected 

Significantly  The  edge  shift  should  approach  zero  if  it  is  possible  to  find  the  limiting 
O’  .rradiance  (or  square  root  of  irradiance)  on  the  bright  side  of  the  edge  rrage, 
except  for  axially  unsymmetnc  aberrations  like  coma.  The  ma)Or  disadvantage  of  using 

the  square  root  is  that  the  apparent  ringing  on  the  dark  side  of  the  edge  image  is 
ennanced. 

Overall,  it  can  be  concluded  that  aoodisation  is  effective  in  imorovmg  the 
performance  of  an  aoerrated  coherent  optical  system  for  edge  objects.  Specifically,  eoge 
n-igmg  is  greatly  reduced,  while  the  edge  'S  not  significantly  snitted  ••elative  to  the 
unaoodised  situation  The  major  drawbacks  of  aoodisation  are  ’hat  ’he  total  ‘lux  m  ’"e 

image  plane  is  reduced,  as  is  the  acutance  of  the  edge 
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Fig. 4. 7  The  acutance  of  the  edge  is  shown  for  various  amounts  of  the  aberrations  of  (a), 
spherical,  (b)  y-axis  coma,  and  (c)  defocus.  The  solid  curves  are  for  the 
unapodised  cases  while  the  dashed  curves  are  for  the  apodised  cases.  The  plots  in 

•.he  !en  column  are  caicuiateo  ‘rom  the  eoge  'mage  rraaiance  wmie  the  otne’’*; 
are  ^rom  the  square  root  of  the  edge  image  irradiance. 


4.2  SLIT  OBJECTS 


A  slit  in  the  object  plane  can  be  described  by 


'(v.)='«(7)= 


1  if  -a  S  y  £  a 
0 

0  otktrwist 


(4.7) 


where  a  is  the  half-width  of  the  slit.  Following  (2.28),  the  inverse  Fourier  transform  of 
(4  7)  is  first  Obtained,  then  multiplied  by  the  scaled  exit  pupil  function  (4.S),  and  finally, 
Fourier-transformed  to  yield  the  amplitude  of  the  coherent  image.  The  image  irradiance 
IS  then  found  as  the  product  of  the  amplitude  times  its  complex  conjugate. 


A  program  called  SLIT. FOR  was  written  in  FORTRAN  on  the  DEC10  computer  of  the 
University  of  Rochester.  An  annotated  copy  of  this  program  is  contained  in  Appendix  1 . 

Typical  results  are  shown  in  Figs.  4.8  through  4.13.  Figures  4.8,  4.10,  and  4.12 
display  the  irradiance  of  a  slice  through  the  image  of  a  slit  when  the  optical  system  has 
various  amounts  (0  0,  0.5,  and  1 .0  waves)  of  defocus,  y-axis  coma,  or  spherical  aberration 

Figures  4  9,4  1' ,  ana  4  1 3  give  the  same  information,  resoectively,  m  terms  of  the  scuare 
root  of  the  'mage  rradiance  Once  agam  separate  plots  of  astigmatism  are  not  shown 
because,  m  one  aimension,  astigmatism  has  'ne  same  functional  form  as  defocus.  The 
•00  olot  m  each  *'Qure  s  ^zr  an  unaoooised  octical  svstem.  while  the  bottom  olot  :s  *or  a 
system  having  a  Gaussian  apodiser  with  an  amplitude  transmittance  described  by  (4  6) 
"  eac"'  z'O’  -e  o'O.nate  'S  n  terms  of  -he  '•eiative  'maoiance  and  the  abscissa  s  ;n  terms 
O’  ”'e  -'o"" a  .ceo  Oistance  For  re'e'-e''ce.  the  position  of  the  geometncaily-imaged  sht 
,  .aea  '  r  jc- a  as  a  ne  of  shor  aatnas 


As  .an  oe  ,ee"  tom  these  ugures,  tne  .mage  ot  a  snt,  like  the  image  of  an  ecge, 

exhibits  the  phenomenon  of  edge  ringing  Since  a  slit  is  composed  of  two  edges,  this  is 
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Fig  4.8  The  coherent  image  of  a  slit  through  an  optical  system  having  0.0,  0.5,  and  1.0 
waves  of  defocus.  The  too  plot  is  for  an  unapodised  ootical  system  while  the 
oonom  plot  i  for  a  system  naving  a  Gaussian  apooiser.  ,.egena; 

■ _ unaberrated. _ _ 0.5  waves. _  _ t  0  waves 
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4.9  The  square  root  of  the  irraoiance  oi  the  conerent  image  of  a  siit  Ihrougn  an 
optical  system  having  0  0,  0  5,  and  10  waves  of  defocus  The  top  plot  is  for  an 
unapodised  optical  system  while  the  bottom  plot  is  for  a  system  having  a 
Gaussian  aoodiser.  Legend: 

_ unaoerrated, _ 0.5  waves, _  _  1.0  waves. 
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The  coherent  image  of  a  slit  through  an  optical  system  having  0.0,  0.5,  and 
1.0  waves  of  y*axis  coma.  The  lop  plot  is  for  an  unapodised  optical  system 
while  the  bottom  plot  is  for  a  system  hoving  a  Gaussian  aoodiser.  Legend 
_ unaberrated _ 0.5  waves,  1.0  waves. 
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"he  square  root  of  tne  irraoiance  gt  ;ne  ^onerent  irrage  of  a  yu  inrcugn  an 

optical  system  having  0  0,  0  5,  and  10  waves  of  y-axis  coma.  The  top  plot  is 
for  an  unapodised  optical  system  while  the  bottom  plot  is  for  a  system  having 
a  Gaussian  apodiser.  Legend: 

_ unaoerratea, _ 0.5  waves.  1.0  waves. 
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The  coherent  image  of  a  slit  through  an  optical  system  having  0.0,  0.5,  and 
1.0  waves  of  spherical.  The  top  plot  is  for  an  unapodised  optical  system  while 
the  bottom  olot  is  for  a  system  naving  a  Gaussian  apodiser  ,sgeno; 

_ unaberrated, _ ; _ 0.5  waves, _  1.0  waves 


SQUARE  ROOT  OF  RELATIVE  IRRADIANCE 


not  surprising.  Also  the  individual  location  of  each  edge  in  the  slit  image  shifts  as  a 
function  of  the  amount  of  aberration.  Likewise,  the  acutance  of  each  edge  varies  with 
the  amount  of  aberration. 

The  effect  of  apodisation  on  the  image  of  a  slit  is  also  similar  to  the  case  of  the 
edge.  The  edge  ringing  is  reduced,  the  edge  shift  is  not  affected  much,  and  the  actuance 
is  decreased  when  an  apodiser  is  used. 

The  only  difference  between  a  single  edge  in  a  slit  image  and  an  isolated  edge  is 
the  effect  of  the  presence  of  the  second  edge.  Depending  on  the  width  of  the  slit 
relative  to  the  width  of  the  impulse  response,  this  effect  can  be  large  or  small.  In  the  case 
illustrated  here  the  effect  is  small,  for  examples  of  smaller  slits,  refer  to  the  experimental 
results  section  of  Chapter  5,  Figs.  5.18  through  5.23. 

For  the  image  of  a  slit,  there  are  two  other  important  measures  of  performance: 
the  width  of  the  slit  image  and  the  translation  of  the  slit  image  as  a  whole. 

The  width  of  the  slit  image  is  defined  as  the  separation  of  the  two  individual  eoges 
of  the  slit  image  (images  here  were  in  terms  of  irradiance)  The  behavior  of  the  width  of 
the  irradiance  of  the  slit  image,  both  with  ana  without  apodisation  ,  is  shown  m  Fig.  4. 14 
(a  ■  c)  for  the  three  different  aberrations,  "^he  image  of  the  slit  ideally  would  have  a 
width  of  80  in  units  of  tne  dimensionless  variable  ■.•.  It  can  be  seen  from  this  figure  that  m 
3:1  rases  the  measured  width  is  less  than  .deai  Aaditionallv.  the  wiath  of  the  :iit  image 
decreases  as  the  amount  of  aberration  increases  If  the  width  were  calculated  from  the 
image  in  terms  of  the  square  'oot  of  -mage  irradiance,  t.ne  width  m  all  cases  woulc  be 


AMOUNT  OF  ABERRATION 


’he  viarn  ‘.ne  rr3aiance  distribution  ot  the  mage  of  a  coherently 

illuminated  slit  is  shown  for  various  amounts  of  the  aberrations  of  (a) 
spherical,  (b)  y-axis  coma,  and  (c)  defocus.  The  solid  curves  are  for  the 
unapodised  cases  while  the  dashed  curves  are  for  the  apodised  cases.  Only 
slit  'mages  that  were  .n  -erms  of  the  square  root  of  rradiance  were 
considered. 


With  the  application  of  the  Gaussian  apodiser,  the  width  of  the  slit  image  will  also 
decrease  with  increasing  amounts  of  aberration.  However,  for  the  cases  of  spherical 
aberration,  defocus,  and  astigmatism,  the  presence  of  the  apodiser  results  in  more 
accurate  width  when  the  amount  of  aberration  exceeds  about  one  wave.  For  the  case  of 
y-axis  coma,  the  presence  of  the  apodiser  doesn't  affect  the  width  of  the  image 
significantly. 

The  position  of  the  slit  image  in  the  image  plane,  besides  being  a  function  of  the 
system's  transverse  magnification,  depends  on  the  type  and  magnitude  of  aberration 
present.  An  examination  of  Figs.  4.9,  4.11.  and  4.13  confirms  that  for  the  symmetric 
aoerrations  (defocus,  astigmatism,  and  spherical)  the  position  of  the  slit  is  correct  (square 
root  of  irradiance  considered);  while  for  the  unsymmetric  aberration  of  coma  the  center 
of  the  slit  image  is  a  function  of  the  amount  of  aberration  present.  This  is  in  agreement 
with  the  definition  of  coma;  i.e.,  the  aberration  of  coma  is  a  variation  of  transverse 
magnification  with  radial  zone  of  the  exit  pupil. 

As  in  the  case  of  an  edge,  the  use  of  apodisation  is  beneficial  when  imagmg 
through  an  .ooticai  system  which  has  thiro-order  aberrations.  The  amount  or  edge- 
ringing  is  reduced  and  the  width  is  not  greatly  affected  when  the  apodiser  is  applied.  As 
oerore,  the  arawDacxs  are  a  .oss  ol  total  rlux  tnrougn  the  mage  piane  and  a  decrease  n 
the  acutance  of  the  edges 
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4.3  TWO-POINT  OBJECTS 

Two-point  sources  in  the  object  plane  can  be  described  by 

‘  ( Vo)  =  ®  ®  ^ 

where  5  is  the  Dirac  delta  function  and  b  is  the  displacement  of  each  point  source  from 
the  optical  axis.  Once  again,  (2.28)  is  employed.  The  inverse  Fourier  transform  of  (4  8)  is 

where  g  =  q/.V/o  The  image  field  amplitude  is  then,  according  to  (2.28), 


A  program  labeled  TPNT  FOR,  which  is  a  derivation  of  PSF  FOR,  was  used  to 
calculate  the  irradiance  distribution  associated  with  (4.10).  An  annotated  copy  of 
TPnT  FOR  IS  in  Appendix  1. 

Typical  results  from  this  program  are  displayed  in  Figs.  4.1S  and  4.16.  These  figures 
disolay  the  irradiance  distribution  of  the  coherent  image  of  two  point  sources  through 
an  jnaoerrated  optical  system  In  each  figure  the  two  points  are  separated  by  (a)  2o  = 
2C,  I'b)  2b  =  30,  and  (c)  2b  =  40.  Fig.  4  i5  shows  the  results  wnen  the  svstem  s 
unapooised  ana  Fig.  4. 1 6  shows  the  results  when  the  system  has  the  Gaussian  apodiser  of 
i-t  5) 

In  the  top  plot  of  each  figure  the  two  points  are  unresolved,  while  for  greater 
separations  m  the  otner  plots,  they  are  resoivea  Comparing  the  miadle  plots  of  the  two 
figures,  it  can  be  seen  that  the  limit  of  resolution,  by  the  Sparrow  limit  of  resolution  is 
greater  for  the  unapodised  case  and  than  for  the  apodised  case,  it  follows  that  the 
apodiseo  case  should  have  a  larger  Sparrow  nmit  of  resolution  because  the  apocisea 
impulse  response  is  wider  than  the  unapodised  one.  The  limit  of  resolution  by  the 


The  irradiance  distribution  in  the  image  of  two  point  sources  separated  m  the 
object  plane  by  (a)  2b  s  20,  (b)  2b  «  30).  and  (c)  2b  •  40.  The  optical  system  Is 
unaberrated  and  unapodised.  The  vertical  axes  are  Kaled  in  the  same 
relative  units  of  irradiance,  while  the  horizontal  axes  are  scaled  n  the  same 
relative  units  of  dimensionless  distance  coordinates. 


The  irradiance  distribution  in  the  image  of  two  point  sources  separated  in  the 
object  plane  by  (a)  2b  «  20,  (b)  2b  «  30),  and  (c)  2b  s  40.  The  optical  system  is 
unaberrated  and  has  a  Gaussian  apodiser.  The  vertical  axes  are  scaled  in  the 
same  relative  units  of  irradiance,  while  the  horizontal  axes  are  Kaieo  m  the 
same  relative  units  of  dimensionless  distance  coordinates. 


Sparrow  criterion  is  charted  in  Fig.  4.17  for  the  third-order  aberrations  and  defocus  both 
with  and  without  a  Gaussian  apodiser.  Data  points  are  not  included  for  higher  amount 
of  y-axis  coma  and  defocus  because  the  point  image  is  not  sharply  peaked  (see  Chapter  3) 
at  these  higher  values. 

In  the  image  of  two  points,  the  centers  of  the  point  images  are  not  necessarily  in 
the  geometrically  predicted  locations.  The  difference  between  their  actual  separation 
ana  their  predicted  separation  is  the  mensuration  error  The  mensuration  error  for  the 
unaberrated  case  is  shown  in  Fig.  4  18  ,  where  for  the  solid  curve  no  apodiser  was  used 
and  for  the  dashed  curve  the  Gaussian  apodiser  was  used.  The  abscissa  of  the  plot  is  the 
expected  separation  of  the  two  points  based  on  geometric  considerations.  The  scale  is  in 
terms  of  the  canoical  coordinate  v.  The  ordinate  is  in  terms  of  the  mensuration  error 
which  is  the  difference  between  the  geometrically  expected  and  the  measured 
separations.  At  low  values  of  image  plane  separation  the  points  are  unresolved,-  as  in  Fig. 
4  15a;  so  that  the  mensuration  error  is  equal  to  the  expected  separation.  At  the  point 
where  the  curves  star  to  turn  up,  the  points  are  just  resolvable  by  the  Sparrow  cnte'^ion 
The  mensuration  error  then  oscillates  about  the  zero  value  before  damping  to  the 
steaay-state  value  of  zero.  It  can  be  seen  that  the  use  of  apoOisation  decreases  tne 
mensuration  error  in  the  regime  where  the  points  are  well-resolved.  However,  the  limit 
of  resolution  has  been  degraded. 

It  can  be  argued  that  the  better  f-esolution  m  the  unaoodised  case  'S  ''o? 
"^■eaningful.  ^his  is  because  the  mensuration  eTor  is  mucn  greater  for  the  jrapodised 
case  So  even  though  the  resolution  of  two  closely-spaced  points  :s  better  in  the 

unapodised  case,  the  apparent  position  of  the  two  points  is  more  likely  to  be  in  error 

As  an  example  of  what  happens  when  aberrations  are  present,  Fig.  4  19  shows  the 
mensuration  error  when  the  optical  system  has  0  5  waves  of  third-order  spherical 


1^1 


AMOUNT  OF  ABERRATION 


The  Soarrow  limit  of  resolution  for  various  values  of  (a)  spherical,  (bl  y-axi$ 
coma,  and  (c)  cefocus.  Data  points  are  not  included  for  higher  amounu  of  y- 
axis  coma  and  defocus  because  the  Ooint  image  is  not  sharoly  peaked  at  these 
higher  values. 


IMAGE  PLANE  SEPARATION 


The  mensuration  error  for  a  coherent  optical  system  having  no  aberrations 
The  solid  curve  is  for  the  unapodised  case  while  the  dashed  curve  is  for  the 
apodised  'ases.  The  mensuration  er'or  s  tne  measured  '.eoaration  of  the  two 
points  minus  the  exoeaed  seoaration  of  the  Gaussian  image  points. 
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aberration.  Both  cases  of  apodisation  are  shown;  the  two  sets  of  curves  of  Figs.  4. 18  and 
4.19  are  nearly  identical.  This  lack  of  difference  is  because  the  effect  of  small  amounts  of 
spherical  aberration  on  the  amplitude  impulse  response  is  to  transfer  energy  from  the 

central  lobe  to  the  first  side  lobe  while  leaving  the  position  of  the  first  dark  ring 

% 

unchanged  (see  for  instance  Fig.  3. IS).  In  this  case,  the  aberration  in  the  otical  system 
has  not  affected  the  general  shape  of  the  curves  very  much.  However,  a'  isation  has 
the  same  general  effect  as  in  Fig.  1.14.  it  decreases  the  mensuration  .  or  where  the 
points  are  resolved,  but  degrades  the  limit  of  resolution 


4.4  IMAGING  THEORY  CONCLUSIONS 


The  utility  of  apodisation  in  improving  the  imaging  performance  of  an  aberratec 
coherent  optical  system  has  been  investigated.  The  imaging  performance  has  been 
evaluated  in  terms  of  the  properties  of  the  images  of  simple  obieas 

The  image  of  an  edge  was  evaluated  ;n  terms  of  the  edge  rirg  ng,  sn  ^v  arc 
acutance  in  general,  as  the  amount  of  a  paaic«iar  aber'atio''  re-ease:  -'e  e:;e 
':ng;ng  ncreased.  the  edge  shift  mcreaseo.  ana  tne  eoge  acutance  aec-ease:  "-e 
apodiser  was  ver/  effective  m  decreasing  the  amount  of  eage  ringing  r  -.re  aoer-aiec 
"rage  of  an  eoge  "he  -inging  couiQ  oe  suoc'esseo  -e  v  'cweve-  s-"  ,  •-e 

optical  system  was  free  from  aberrations. 
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"he  evaluation  criteria  for  tne  mage  or  a  siit  we'-e  '•'ng  -'g  or  :-'e  ro.vicua^  eeges 

of  the  silt,  the  apparent  width  of  the  slit,  anc  the  amour'  s..'  apcea-ec  tc  oe  s~  "ec 

The  ringing  oenavior  of  the  edges  in  the  slit  was  quite  simi'ar  to  the  rmgmg  n  •.'re  s.noie 

eoge.  This  .s  to  oe  expecteo  since  me  edge  separation  m  tne  sut  cons.ce'ec  .s  moc- 


larger  than  the  period  of  the  ringing.  The  width  is  a  funaion  of  the  individual  shifts  of 
me  two  eeges.  The  ocation  of  the  sm  image  is  correct  tor  tne  symmetnea.  aoerrat.ons 


and  is  shihed  for  the  unsymmelrical  ones  like  coma.  The  effect  of  the  apodiser  is  very 
similar  to  the  effect  seen  in  the  case  of  the  edge. 

When  the  image  of  two  closely-spaced  pomts  is  evaluated,  the  important 
parameter  is  the  mensuration  error.  The  mensuration  error  for  an  unaberrated  system  is 
an  oscillatory  function  which  decreases  to  zero  as  the  point  objects  are  separated  in  the 
object  space.  The  addition  of  spherical  aberration  to  the  system  did  not  significantly 
affect  this  behavior  The  apodiser  was  quite  effective  m  damping  the  oscillatory  nature 
of  the  mensuration  error  The  resolution  limit  however  suffered  as  a  consequence. 

It  was  not  intended  that  this  chapter  contain  an  exhaustive  catalog  of  theoretical 
results  The  intention  was  to  present  enough  data  to  gain  an  understanding  of  the 
physical  processes  It  was  not  known  previously  that  apodisation  would  be  as  effective  m 
improving  the  performance  of  aberrated  systems  as  it  was  in  improving  unaberrated 
ones  Now  however,  m  light  of  the  model  developed  and  the  results  presented,  it  can  be 
cone  uoed  tinat  the  apodiser  is  effective  m  improving  the  imaging  performance  of 
j0r"3’ec  cone'ent  ootical  systems.  The  reason  for  this  is  that  the  apodiser  attenuates  bv 
ao'  D'c’  or^  ’nose  -eg  ens  o''  ’.he  exit  pupil  where  the  aberrations  most  affect  the  phase 
•’ars’”  tteo  waver'cnt 

■^he  results  presented  here  were  generated  from  a  theoretical  model  which  has 
—  an,  ass,,rnptions  ana  approximations  There’ore,  the  model  which  generateo  them 
-ooos  *0  be  ’estec  encer  m, entail v  before  ’hev  can  pe  ’  jilv  acceptea 


CHAPTER  5  EXPERIMENTS 


The  purpose  of  conducting  these  experiments  was  to  test  the  theory  developed  m 
Chapter  2  and  applied  in  Chapters  3  and  4.  The  experimental  setup  and  the  data 
acquisition  procedure  is  fully  described  below  before  the  results  are  given.  The  results 
are  given  m  two  separate  sections  which  parallel  the  subject  matter  of  Chapters  3  and  4 
That  is,  the  data  sections  deal  separately  with  the  impulse  response  and  the  imaging  of 
Simple  objeas 

5.1  EXPERIMENTAL  DESIGN  AND  DATA  ACQUISITION 

The  basic  experimental  configuration  is  snown  m  Fig.  5.1.  This  configuration  was 
designed  to  measure  the  irradiance  impulse  response  of  the  optical  system  formed  by 
lens  LI,  the  apodiser,  the  iris,  and  lens  L2.  Simple  modifications  of  this  setup  made  it 
possible  to  measure  the  image  irradiance  distribution  of  simple  objects.  The 
soecifications  of  the  equiament  used  ’n  these  configurations  are  m  Taoie  5  1 

5.1.1  The  Impulse  Response  Configuration 

In  the  basic  optical  system  of  F-g  5  :  cor'e''er't  Mlummation  originated  from  a  C  5 
mW,  linearly  polarized  HeNe  laser.  T"e  aser  beam  -vas  .n  the  lowest  oraer  fansverse 
mode  (TEMaq)  and  m,  at  most,  three  'ongitudmal  modes  with  a  center  wavelength  of 
0  6328  microns.  The  beam  had  a  diameter,  between  the  l/e^  irradiance  points,  of  0.64 
mm  and  a  divergence  of  1.3  mrad. 
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Fig. 5.1  The  experimental  configuration  used  to  measure  the  irradiance  impulse 
response  of  the  ootical  system  formed  py.Ll,L2,  the  ins,  and  the  aoodiser 


K  r 


Designation 


Source 


■'.n  »t'’  KT-  »-v  n.' 


Description 


Laser 

Absorbing  filter 


Hughes  H3022P  5  mW,  linearly  polarized 


Wratten  neutral  density,  various 
values 


Ml, M2 

Newport 

10D10 

1 "  dia.  pyrex  1«  surface  mirrors,  A/10 

Spatial  filter 

Jodon 

LPSF-100 

10X  microscope  objective,  lOprn 

pinhole 

LI 

Jaegers 

Telescope  doublet,  various  focal 
lengths 

X\':\ 

Apodiser 

See  App.  3 

52,6 

Gaussian  amplitude  apodiser 

*•  , 

Iris 

Newport 

ID- 1.5 

18  curved  leaves 

Is 

’  .  “  Si 

C.-w* 

L2 

Depended  on  particular  test 

configuration 

L3 

J.E.A. 

780094  . 

20X,  0.40  N.A.  microscope  objertive 

Detector  Array 

Fairchild 

I-SCAN 

2S6-element  CCD  linear  array  with 

electronics  ' 

PDi 

Ealing 

Point  Diffraction  Interferomete’’,  *.  i  2  j 

pinhole  ^ 

XV'' 

1^ 
^ "  * " 

Camera 

Pentax 

KIOOO 

Body  used  without  lens 

Table  5.1  The  equipment  used  m  the  experiment. 


The  first  optical  element  the  beam  encountered  was  the  absorbing  filter.  This 
filter  was  actually  one  or  several  Wratten  neutral  density  filters  of  various  values.  The 
filters  was  used  to  adjust  the  irradiance  level  in  the  plane  of  the  detector  array. 

Two  first-surface  mirrors.  Ml  and  M2,  directed  the  beam  into  the  spatial  filter.  The 
mirrors  had  a  surface  flatness  of  X/10.  The  spatial  filter  was  formed  by  a  lOX,  0.25  N.A. 
microscope  objective  and  a  10  micron  diameter  pinhole. 

The  pinhole  acted  as  a  point  source  for  the  optical  imaging  system  formed  by  LI, 
the  apodiser,  the  ins,  and  L2.  The  image  of  the  point  source  through  this  system  was  the 
irradiance  impulse  response  of  the  system  Lens  LI  was  a  telescope  objeaive  with  a 
diameter  of  78  7  mm  and  a  focal  length  of  495  mm.  This  iris,  while  adjustable,  was  used 
with  a  nominal  diameter  of  25mm.  If  LI  was  used  to  focus  a  plane  wave  limited  by  the 
ins  at  the  nominal  diameter,  the  diameter  of  the  resulting  Airy  spot  would  be 
approximately  30  microns.  Thus,  it  is  clear  that  the  10  micron  pinhole  was  indeed  acting 
as  a  point  source  for  this  system.  Furthermore,  a  simple  calculation  reveals  that  the 
irradiance  of  the  wavefront  from  the  pinhole  at  the  lens  LI  varied  by  no  more  than  2 
oercent  over  a  central  circular  region  of  25  mm  m  ciameter  Leris  LI  was  a  hign  auality 
telescope  oojective  purchased  from  jaeger  Optical  Company  When  used  at  this  recuced 
acemure  the  ens  would  oe  e.xceccec:  to  .'e  es'.ent.a.Iy  aoerration-iree  A  measure'^?''’ 
of  the  wavefront  after  the  lens  LI  confirmed  this  expectation 

'he  apodiser  was  one  of  several  wnose  amplitude  transmittance  was  .nominally 
Gaussian  The  design,  construction,  ano  characterization  of  these  aoodisers  is  detailed  in 
Appendix  2 

Lens  L2  was  the  final  element  in  tne  optical  imaging  system.  During  tne 
experiments  different  lenses  were  used  m  this  position.  The  particular  lens  used 
aepenoeo  upon  the  aoerrations  desired. 
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The  image  of  the  point  source  was  formed  in  the  image  plane.  Two  different 
configurations  were  used  to  characterize  this  image. 

The  first  is  shown  on  the  optical  axis  below  the  image  plane  in  Fig.  5.1.  It  consisted 
of  a  20X  microscope  objective  (L3)  and  a  linear  detector  array.  The  microscope  objective 
acted  as  a  magnifier  Magnification  was  needed  because  the  scale  of  the  image  of  the 
point  in  the  image  plane  was  too  small  to  be  adequately  sampled  by  the  detector  array 
Lens  L2  generally  had  a  focal  length  of  about  300  mm  Thus,  the  diameter  of  the  Airy 
soot  in  the  image  plane  was  approximately  19  urn.  The  detector  array  was  a  Fairchild 
CCDl  1 1  linear  array  of  charge  coupled  detectors  (CCD),  each  of  area  17x13  pm,  on  13 
pm  centers.  Hence  the  magnification  provided  by  the  microscope  objective  L3  was 
needed.  The  array  was  part  of  a  card  of  electronics  marketed  by  Fairchild  under  the 
name  l-SCAN.  The  data  flowing  from  this  array  was  collected  by  a  microcomputer 
Appendix  3  contains  a  complete  description  of  the  data  collection  equipment.  This 
configuration  resulted  in  a  linear  array  of  samples  of  the  irradiance  through  the  center  of 
the  image  of  the  point  source. 

The  second  configuration  used  to  measure  the  'mage  'S  shown  shifted  to  the  right 
of  the  first  configuration  in  Fig.  5.1.  it  consistea  of  a  ooint-diffraction  interferometer'3 
(POD,  lens  L-.  and  a  camera  back,  the  PDI  s  a  oaruallv  transmitting  olate  with  a  ‘otallv 
transmitting  circular  region  in  its  center  The  circular  region  was  on  the  order  of  the  size 
of  the  -Lirv  soot  generated  bv  LZ.  '.Vhen  'he  ’01  -vas  oiaceo  n  -he  mage  oiare.  “'e 
wavefront  from  l2  passed  through  the  i^Dl  and  was  split  mto  two  wavefronts  One  was  a 
reoiica  of  the  incident  wavefront  with  .'educeo  rradiance,  wm.e  the  otner  was  a 
diverging  toherical  wavefront  generated  by  -he  email  :ircuiar  -egion.  "'■'e'te  ■.vo 
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wavefronts  interfered  and  the  resultant  interference  pattern  contained  information 

about  the  aberrations  in  the  ootical  system  Lens  Lb  was  adjusted  to  '.hat  'he  exir  ;.jcil 


108 


was  imaged  orto  the  film  plane  in  the  camera.  Thus  the  recorded  interferogram 
contained  information  about  the  aberrations  in  the  exit  pupil. 

S.1.2  The  Imaging  Configurations 

For  the  experiments  involving  the  image  of  two  point  sources,  the  configuration 
of  Fig.  5  1  was  used  with  little  modification.  However,  when  imaging  an  extended  object 
such  as  an  edge  or  a  slit  this  configuration  had  to  be  extensively  modified  as  shown  in  Fig 
5.2. 

5. 1.2.1  Imaging  of  Two-point  Objects 

A  set  of  five  double  pinholes  was  purchased  for  use  in  the  two-point  imaging 
experiments.  They  were  placed,  one  at  a  time,  m  the  plane  of  the  pinhole  of  the  spatial 
filter  of  the  configuration  of  Fig.  5.1.  The  double  pinholes  were  mounted  on  the  same 
type  of  circular  magnetic  holders  used  for  the  spatial  filter.  The  spatial  filter  lens  was 
moved  longitudinally  away  from  the  filter  location  so  that  each  pinhole  was  illuminated 
uniformly  The  detector  array  was  centereo  on  the  resulting  image  irradiance  pattei-n 
’he  separation  of  “he  oeaKS  in  the  image  was  determined  by  counting  the  num.oer  of 
pixels  between  the  oeaks. 

The  double  pinholes  were  manufactured  by  burning  two  closely-spaced  sm.all 
noies  witn  a  focusseo  laser  oeam  ’he  "loies  wei'e  scecitieo  to  oe  10  microns  in  oiameter 
’he  spacmgs  wei-e  soecified  to  oe  •12.35.23.2’.  and  ’<1  microns.  The  actual  diameters  and 
soacmgs  were  oetermineo  experimentally. 

The  diameters  of  the  individual  holes  were  determined  by  observing  the  far-field 
diffraction  oattern  when  the  hole  was  iMummated  f'om  behind  Essentiallv.  each  hole 


was  set  up  as  a  slightly  defocused  spatial  frequency  filter.  The  incident  wavefront  was 
approximately  flat  in  phase  and  uniform  in  irradiance.  The  far-field  pattern  was  then  an 
Airy  diffraction  pattern.  The  width  of  the  first  dark  ring  (diameter  of  the  Airy  spot)  was 
measured  and  the  width  of  the  diffracting  aperture  (pinhole)  was  inferred  using  Fourier 
optics  techniques.  The  pinhole  diameters  were  found  to  be  9.5  i  0.5  pm. 

The  separations  were  found  in  a  similar  manner.  For  this  measurement,  both  holes 
of  each  double  pinhole  were  illuminated  by  a  beam  which  was  approximately  uniform  in 
phase  and  modulus.  The  resulting  far-field  pattern  was  a  system  of  Young’s  interference 
fringes.  The  spacings  of  these  fringes  were  related  to  the  separation  of  the  holes. 
Assuming  that  one  pinhole  separation  was  known,  the  others  could  be  determined  by 
the  ratios  of  the  fringe  spacings  in  the  various  Young's  interference  patterns. 

The  separation  of  the  two  holes  in  the  double  pinhole  with  the  largest  labeled 
separation  was  measured  using  a  powerful  optical  microscope.  For  this  purpose,  a  Leitz 

Or.hoiux  microscope  was  used  It  had  a  SOX  objective  (N  A  =  0  65)  and  a  lOX  eyepiece 
‘or  an  overall  visual  .magnification  of  SCO.  "he  s.malles:  'esolvable  linear  dimens.on  m  '.re 
oDiect  plane,  by  ’.ne  Sayieign  criterion  was  c  =  O.si.'v'J.A  .with  .‘,  =  0.3  microns  ’.n.'S  gives 
;  :  .i7  mic.'ons.  '-e  measureo  .eDaraiio.n  '.he  'wo  'oies  n  rms  rase  -vas  -Z  :  ' 

micron.  Notice  that  the  measured  separation  is  much  greater  than  the  resolution  hm.it  so 
•.hat  tne  mensuration  error  discussed  m  Chaoters  2  ana  -  was  negligible. 

5. 1 .2.2  Imaging  of  Edge  and  Slit  Objects 

The  .oniiguration  or  rig,  5  2  was  .leeaea  ’or  me  .maging  oi  the  eage  ana  s..: 
objects  because  of  the  presence  of  speckle.  Speckle  is  the  random  irradiance  fluctuations 
seen  .n  me  .mage  irraoiance  distribution,  t  :s  causea  oy  ranaom  variations  n  the 


nominal  surface  shapes  of  mirrors  and  lenses,  as  well  as  dust  particles  on  these  surfaces. 
Speckle  was  reduced  in  this  experiment  primarily  by  using  fewer  optical  elements,  one 
lens  instead  of  three,  between  the  objea  and  the  image  planes.  Thus  there  were  fewer 
sources  of  speckle. 


The  effects  of  speckle  were  also  controlled  by  carefully  choosing  the  placement  of 
the  iris  The  characteristic  size  of  the  speckle  pattern  was  determined  by  the  limiting 
aperture  in  that  part  of  the  system  which  followed  the  source  of  the  speckle.  For 
instance,  if  the  source  of  the  speckle  was  m  the  object  plane  of  Fig.  5.2,  the  iris  would  be 
the  limiting  aperture  and  the  characteristic  size  of  the  speckle  in  the  image  plane  would 
be  the  Airy  spot  diameter  of  the  system  in  the  prcse''t  case  however,  the  object  was  one 
or  a  pair  of  razor  blades  (no  source  of  specl"e  except  perhaps  along  the  edges),  and  the 
incident  wavefront  had  just  been  spatially  filtered  so  it  was  relatively,  if  not  completely, 
free  from  random  phase  flucuations.  This  then  left  as  a  major  source  of  speckle,  the  lens 
Li.  But  the  iris  was  not  the  aperture  stop  for  the  optical  sytem  following  the  lens.  The 
system  following  the  lens  was  simply  a  free-space  oropagation  system  Thus  the 
characteristic  size  of  the  soeckle  'rom  sources  m  or  on  the  lens  L '  was  much  .arger  than 
the  speckle  size  from  ooject  sources.  The  effect  was  a  'eduction  of  the  spec*,  e  wn:cn 
-•jr-jrar;  - 3 r  — ' a  .  .vith  ' .0 e  jesTeo  ..mages  ''■'e  :cecv.e  "  f^e  mage  oia.^e  .-.'■i  e  'o' 
completely  removed,  was  considerably  reduced  by  this  techmoue. 


5.1.3  Alignment 

■^he  proper  alignment  of  the  various  experimental  configurations  was  crucial  to 
•no  cuccets  of  ’he  experiment.  A  misaligned  tvstem  .vould  nave  unwanted  aberr.ation'.. 
The  alignment  procedure  was  based  on  a  technique  developed  by  Taylor  and 
■'nomosont9  The  essence  of  the  'dea  was  'o  first  establish  unambiguously  an  optical  .axis 


and  then  align  each  element  to  that  axis.  The  procedure  outlined  here  is,  for  clarity, 
specific  to  Fig.  5.1.  However,  the  techniques  are  perfectly  general  and  were  used  to  align 
all  of  the  experimental  configurations. 

Initially,  the  optical  axis  was  established  with  '  t  laser,  the  turning  mirrors  Ml 
and  M2,  and  a  target  placed  in  the  detector  piano  .  g.  5.1).  The  target  was  a  black 
cross  inked  on  a  white  card.  The  card  was  mounted  so  as  to  be  perpendicular  to  the 
optical  table  surface.  The  center  of  the  cross  was  4  44  inches  above  the  surface.  Placing 
the  cross  m  front  of  the  laser,  the  output  aperture  of  the  laser  was  adjusted  to  be  4  44 
inches  above  the  surface.  Then  the  cross  was  moved  to  a  position  in  front  of  Ml  and  the 
laser  was  tilted,  about  the  outout  aperture,  to  center  the  laser  beam  on  the  cross.  Next, 
the  target  was  moved  just  in  front  of  M2  and  then  Ml  was  rotationally  adjusted  to  the 
center  of  the  laser  beam  on  the  cross.  Then  the  target  was  moved  to  the  detector  plane 
and  M2  similarly  adjusted. 

Next  a  white  painted  sheet  of  metal  with  an  aperture  that  was  smaller  than  the 
laser  beam  was  placed  between  M*  and  the  spatial  filter  location.  The  small  hole  was 
centered  on  the  beam  by  observing  and  centering  the  Circular  diffraction  pattern  from 
•he  hole  on  the  cross  This  hole  was  used  as  a  reference  for  aligning  the  lens  m  the  spatial 
•'liter 

Before  the  spatial  filter  was  aligned,  a  crosshair  formed  by  two  thin  metal  wires 
.vas  aiacecJ  between  'he  position  of  the  soatiai  'ii’er  and  position  of  iens  "he 
C'bssnair  was  centered  on  the  optical  axis  by  observing  and  centering  the  shadow  of  the 
crosshair  on  the  target  cross.  ”he  center  of  ’he  hole  n  the  metal  sheet  and  ^he  crosshairs 
•hen  defined  ‘he  optical  axis. 

The  spatial  filter  was  modified  so  that  the  lens  and  the  pinhole  were  separately 
adiustable.  After  mounting  the  spatial  ‘liter  lens  'n  its  approximate  location  it  was 


adjusted  in  tilt  by  observing  the  refleaed  light  as  it  impinged  upon  the  white-painted 
metal.  The  reflected  light  formed  a  circular  interference  pattern  which  was  centered  on 
the  hole  by  adjusting  the  tilt  of  the  lens  Simultaneously,  the  transmitted  beam  which 
again  cast  a  shadow  of  the  metal  crosshairs  onto  the  target  was  aligned  to  the  target 
cross  by  adjusting  the  lens  in  translation.  When  the  reflected  circular  diffraction  pattern 
was  aligned  on  the  hole  and,  at  the  same  time,  the  shadow  of  the  crosshair  was  aligned 
on  the  target  cross,  the  alignment  of  the  spatial  filter  lens  was  complete. 

The  spatial  filter  pinhole  was  then  pul  m  place  and  aligned  in  the  usual  manner  by 
maximizing  the  power  throughput  while  minimizing  the  transmission  of  higher  spatial 
frequencies.  When  done  correctly,  the  shadow  cast  by  the  metal  crosshairs  remained 
centered  on  the  target  cross.  The  metal  sheet  with  the  small  hole  was  then  removed  At 
this  ooint  the  optical  axis  was  defined  by  the  pmhole  and  the  metal  crosshairs 

Next,  the  mam  elements  of  the  optical  system  were  put  in  place  and  aligned 
These  elements  were  lenses  L1  and  L2,  the  iris  and  the  apodiser. 

first,  l2  was  put  in  place.  It  was  adjusted  in  translation  by  observing  and  centering 
the  shadow  of  the  metal  crosshairs  on  the  target  cross  The  rotational  alignment  was 
acco.molisned  by  jsmg  the  Boys  points  Boys  points  are  images  of  the  soaiiai  filter 
oinhole  formed  by  reflections  from  the  surfaces  o^  the  lens  These  ooints  occur  both 
oe’ore  and  after  the  lens  The  position  of  these  points  which  occur  before  the  tens  are 
oar-  cuiariy  sensitive  to  'otationai  jdiustr'e"t'?  sio  sligning  "nese  oomts  a  f.-e 
metal  mesh  screen  was  interposed  between  the  lens  12  and  the  pmhole.  This  screen  was 
placed  m  tne  plane  of  one  of  the  reflected  Boys  points.  Thus,  a  small  bngnt  spot 
aooeared  on  the  screen.  The  shadow  of  the  metal  c-osshair  also  aoceared  on  -he  screen 
The  lens  was  properly  aligned,  in  a  rotational  sense,  when  the  Boys  point  was  centered 


on  the  shadow  of  the  crosshair  and  simultaneously  the  transmitted  shadow  of  the 
crosshair  was  centered  on  the  target  cross.. 

Lens  LI  was  then  put  in  place  and  aligned  in  an  analogous  manner.  Additionally, 
lens  L2  was  adjusted  along  the  optical  axis  so  that  it  was  exactly  its  focal  lerigt^  r;  ■ 
from  the  pinhole.  This  was  done  by  placing  a  shear  plate  between  lens  LI  ar...  ...  .^nd 
observing  the  reflected  interference  pattern.  The  lens  was  adjusted  correctly  in  a 
longitudinal  sense  when  the  interference  fringes  were  completely  fluffed  out. 

Next  the  iris  was  put  in  place  quite  close  to  lens  L2.  It  was  aligned  in  translation  by 
observing  its  shadow  cast  on  the  target.  The  apodiser,  when  used,  was  aligned  similarly 
and  was  placed  as  close  as  possible  to  the  iris. 

At  this  point  the  optical  system  was  completely  aligned  Depending  on  the 
particular  need,  either  L3  and  the  detector  array  or  the  PDI,  L4,  and  camera  was  then 
aligned  to  the  optical  axis. 

As  a  test  of  the  accuracy  and  repeatability  of  alignment,  the  same  optical  system 
was  aligned  from  scratch  several  times.  Each  time  the  aberrations  were  measured.  The 
optical  system  (see  Fig.  5. 1 )  consisted  of  a  Jaegers  #1158  lens  (f  =  495  mm)  as  LI,  a  Jaegers 
#  953  lens  (f  =  390  mm)  as  L2,  and  an  ;ris  with  a  ciam.eter  of  0  525  ncnes.  The  measui  eo 
third  order  aberrations,  m  terms  of  Siedel  aberration  coefficients,  for  each  of  these  tests 
are  taouiated  n  Taoie  5.2.  Two  conclusions  can  be  drawn  from  tt'.e  data  m  this  taoie. 
^I'-st.  because  of  'he  'ow  magnitude  of  each  aberration,  it  can  be  concluded  that  the 
lenses  were  of  high  quality  and  the  system  was  aiigned  well  each  time.  5econd,  because 
t’^'e  different  alignments  produced  very  similar  results,  it  can  be  concluded  that  the 
alignment  procedure  was  repeatable. 


Astigmatism 

0.0 

0.0 

0.0 

Coma 

0.1 

0.1 

0.1 

Spherical 

-0.1 

-0.0 

-0.1 

Table  5  2  The  measured  third-order  Seidel  aberrration  coefficients  after  several 
augments  of  the  same  optical  system.  The  units  of  the  tabulated  values  are 
wavelengths  (X  =  0.6328  microns) 

5.1.4  Data  Acquisition 

After  a  particular  optical  system  had  been  set  up  and  aligned  the  data  acquisition 
qrocess  could  proceed  The  raw  data  were  ;n  two  fo.'ms.  photograohic  and  electro-optic 
"‘'ese  will  be  discussed  separately.  To  clarify  this  orocess  a  carjcuiar  ootical  system  will 
oe  analyzed  as  an  examole.  For  this  ourpose.  the  ooticai  system  of  ►•g.  3  t  A/as 
"?''s:''jc:ed  whe’-e  '.t  was  a  .'aeqer's  ^  '138  le-'s  (f  =  -i9S  m.m).  L2  was  a  ,'aeqe’- s  #  938 
ens  (f  =  390  mm),  and  the  apodiser  was  #6  These  comoonents  were  selected  anc 
oriented  to  .acmeve  a  system  having  small  amounts  of  aberrations  wnicn  would  'esuit  n 
an  Airy  diffraction  pattern  in  the  detector  plane 

5. 1 .4. 1  Electro-optic  Data 

Electro-optic  data  were  gene’ated  by  operating  the  linear  CCD  detector  array 
wnen  that  array  was  in  the  piane  ot  the  diffraction  pattern  or  interest.  The  detector 


array  and  the  associated  20X  microscope  objective  were  installed  and  aligned  as 
described  in  section  5.1.1  and  5.1.3. 


The  data  were  in  the  form  of  256  eight-bit  words;  one  word  for  each  element  m 
the  detector  array  Such  a  grouping  was  called  a  frame  of  data.  Eight  frames  of  data 
were  taken  each  time  data  was  collected. 

A  single  frame  of  data  from  a  test  of  the  example  system  is  shown  in  Fig.  5.3a.  The 
optical  configuration  was  that  of  the  example  system.  The  vertical  axis  is  relative 
irradiance  plotted  logarithmically  and  the  horizontal  axis  is  normalized  distance  v.  The 
scales  are  the  same  for  all  the  plots  in  Fig.'S  3. 

The  digitized  nature  of  the  data  is  evident  from  Fig  5.3a  There  are  256  (28) 
possible  values  of  relative  irradiance  because  the  word  size  is  eight  bits.  Thus,  the 
minimum  resolvable  signal  difference  is  1/256  =  0.0039  relative  to  a  maximum  signal 
value  of  1.0.  This  minimum  resolution  is  seen  as  discrete  steps  in  the  data  of  Fig.  5.3a 
The  steps  are  of  different  vertical  size  because  the  vertical  axis  was  plotted 
logarithmically. 

When  eight  consecutive  frames  of  data  were  ave’^aged  the  result  was  the  curve  m 
f’c.  5  3b,  The  integration  time  for  eacn  'ramie  was  2  msec,  'he  data  after  this 
averaging  orocess  were  considerably  smoother. 

In  Fig.  5.3c.  the  averaged  data  (pictted  as  tnangies)  are  shown  relative  to  an  Airy 
oattern.  The  difference  between  the  data  and  the  Airv  oattern  never  exceeds  the 
.minimum  resolution  of  the  system, 

in  Fig  5  3d,  the  averaged  data  are  snown  relative  to  a  central  slice  of  the 
irradiance  imoulse  resoonse  of  an  ootical  svstem  having  the  aberrations  measured  usmc 
the  PDI.  The  theoretical  plot  was  generated  using  the  model  described  in  Chapter  3 


normalized  distance  V 

3  A  comparison  of  theoretical  and  experimental  results,  (a)  A  single  frame  of  data, 
'bi  "he  overage  of  .»igni  frames -if -ata  [c)  ^  comoarison  of ‘.ne  oata  in  :b)  *o  an 
Airy  pattern,  (d)  A  comparison  of  the  data  in  ((j)  to  a  central  slice  through  the 
impulse  response  predirted  using  the  measured  aberrations.  The  vertical  axis  for 
each  plot  is  relative  irradiance(I)  plotted  on  a  logarithmic  scale.  The  horizontal 
axis  for  aacn  aiot  s  normaiizsd  distanca  v.  "lie  :caies  are  ‘.ne  same  on  ail  of  ‘he 
plots. 


A  photograph  of  the  irradiance  impulse  response  through  this  optical  system  is 
shown  in  the  top  of  Fig.  5.4.  In  this  case  there  was  no  apodiser  and  essentially  no 
aberrations  (only  those  of  column  3  of  Table  5.2).  The  photograph  should  be  and  is 
approximately  ah  Airy  diffraction  pattern.  The  plot  in  the  bottom  of  this  figure  is  the 
calculated  modulus  of  the  amplitude  impulse  response  with  the  above  aberrations  and 
no  apodiser.  The  plot  was  generated  by  the  model  described  in  Chapter  3  (PSF.FOR)  with 
the  measured  aberrations  as  inputs. 

Two  different  apodisers  were  used  at  different  times  during  these  experiments. 
When  they  were  put  into  the  system,  the  impulse  response  was  significantly  modified. 
Experimental  data  in  the  form  of  a  central  slice  through  the  impulse  response  for  one  of 
these  apodisers  (#6)  is  shown  in  Fig.  5.5.  The  broken  curves  represent  the  experimental 
data  from  two  data  collections  when  the  detector  array  was  rotated  about  the  optical 
axis  by  90^  for  one  of  the  collections  relative  to  the  other.  The  theoretically  prediaed 
irradiance  impulse  response  is  shown  as  a  solid  line  in  this  figure.  The  predictions  were 
generated  by  the  model  described  in  Chapter  3  with  the  aberrations  measured  with  the 
?0;  and  with  the  apodiser  parameter  descnoed  by  (4  6)  The  apodiser  can  be  descnoea 
oy  'he  function  Ti ",]  =  exD(-G r;^), where,  when  G  =  3.  the  apooiser  oescrioea  ay  this 
function  was  m  place  in  the  system  and,  when  G  =  0,  there  was  no  apodiser  in  place  This 
shorthan  notation  will  be  used  m  future  references  to  the  existence  of  an  apooiser  For 
these  data  the  iris  diameter  was  set  to  be  0.525  inches.  Details  about  how  the  apodisers 
were  designed,  constructed,  and  tested  are  contained  in  Appendix  2. 

A  photograph  of  an  apodised  (G  =  3)  ana  essentially  unaberrated  irraCiance 
imouise  resDonse  is  shown  in  Fig.  5.6.  The  optical  system  had  the  third'Order  aberrations 


.he  pnotograph  in  the  top  o»  this  rigure  is  of  the  irradiance  impulse  response  ot 
an  unapodised  and  essentially  unaberrated  optical  system.  The  aberrations  of 
this  system  are  listed  in  column  3  of  Table  6.2.  The  plot  on  the  bottom  of  this 
figure  is  the  calculated  modulus  (vertical  axis)  of  the  impulse  response  through 
the  same  optical  system.  The  horizontal  axes  have  the  same  units  of  normalized 
distance. 


normalized  distance  V 


Fig.S.S  Theoretical  (solid  curve)  and  experimental  (broken  curves)  plots  of  central  slices 
of  the  irradiance  impulse  resoonse  for  an  aoodised  (G  =  3)  and  unaberraied 
optical  system. 


Fig. 5. 6  The  photograph  in  the  top  of  this  figure  it  of  the  irradiance  impulse  respor>se  of 
an  apodised  (G  *3)  and  essentially  unaberrated  optical  system.  The  aberrations 
of  this  system  are  listed  in  column  3  of  Table  6.2.  The  apodiser  is  described  by 
(4  6)  and  shown  in  Fig.  3.9.  The  olot  on  the  bottom  of  this  figure  is  the  calculated 
modulus  (vertical  axes)  of  the  impulse  response  through  the  same  optical  system 
The  horiiontal  axes  have  the  same  uni'ts  of  normalized  distance. 
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listed  in  column  3  of  Table  5.2.  The  apodiser  was  #6  and  the  iris  diameter  was  0.525 
inches  The  plot  in  the  bottom  of  this  figure  was  generated  with  the  same  aberrations 
and  apodisation  by  the  computer  model  PSF  FOR. 

At  this  point,  it  is  clear  that  the  experimental  apparatus  was  capable  of  making 
measurements  which  were  fine  enough  to  measure  the  important  details  of  the 
diffraction  patterns  produced  by  these  optical  systems.  The  measurements  were 
accurate,  precise,  and  replicable.  The  measured  patterns  agreed  well  with  the  prediaed 
patterns . 

5. 1.4. 2  Photographic  Data 

The  data  recorded  by  photographic  means  were  the  interference  patterns 
generated  from  the  point-diffraction  interferometer  (PDI)  and  analog  data  from  the  I- 
SCAN  detector  displayed  on  the  oscilloscope  and  photographed. 

The  PDI  was  used  to  measure  the  aberrations  in  the  optical  system.  The  PDI,  lens 
L4,  and  the  camera  back  were  installed  and  aligned  as  described  in  the  previous  sections 
■^his  was  done  on  every  data  run.  The  interference  pattern  was  recordeo  on  Kodak  Plus-X 
Pan  35  mm  film,  which  was  developed  according  to  tne  manufacturer  s  instructions.  Tre 
resulting  negative  was  used  to  print  a  positive  .mage  onto  an  8"  x  iO"  sheet  of  KooaK 
’oiycontrast  RC  i  ‘lim,  wnicn  also  was  oeveiooec  accoroing  to  the  manufacture!’’ i 
instructions. 

An  interterogram  from  the  example  system  ,s  dispiayeo  .n  Fig.  3.7.  These  oaia 
wei-e  colleaed  with  the  iris  diameter  set  to  be  0.525  inches.  The  shape  of  the  fringes  anc 
their  spacing  are  related  to  the  aberrations  m  the  optical  system  These  aberrations  we’-e 
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found  by  digitizing  the  coordinates  of  central  points  along  the  fringes  The 
interferogram  as  well  as  the  digitized  points  for  this  example  are  shown  on  the  in  Fig 
5.7.  The  digitized  points  were  the  input  data  to  a  fringe  analysis  program  wrinen  by 
Professor  James  Wyant  of  the  University  of  Arizona.  One  of  the  outputs  of  the  program 
is  the  aberration  values  in  terms  of  the  Zernike  as  well  as  Siedel  coefficients.  The  third 
order  Siedel  coefficients  for  the  interferogram  of  Fig.  5.7  are  listed  in  column  3  of  Table 
5  2 

At  times,  instead  of  collecting  the  l-SCAN  data,  digitizing  it,  and  processing  it  with 
the  mic'ocomputer,  the  analog  data  from  the  i-SCAN  detector  was  displayed  on  the 
oscilloscope  where  it  was  photographed.  This  analog  data  was  then  compared  to  the 
theoretical  predictions.  This  method  of  data  collection  was  considerably  more  simple  to 
accomplish  than  the  totally  electro-optic  method  of  section  5. 1.4.1. 
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S.2  EXPERIMENTAL  EVIDENCE 

As  stated  in  the  beginning  of  this  chapter,  the  purpose  of  the  present  experiments 
was  to  adequately  test  the  theory  developed  in  Chapters  2,  3,  and  4.  The  development  of 
that  theory  involved  several  assumptions  and  approximations.  Extensive  predictions 


about  the  behavior  of  aberrated  coherent  imaging  systems  and  the  utility  of  apodisation 
were  made  based  on  this  theoretical  development.  Thus  it  was  necessary  to  perform 
credible  experiments  to  establish  the  validity  of  the  model  used  to  make  those 
predictions. 

In  the  previous  sections  some  experimental  results  were  presented  The 
experimental  data  on  the  irradiance  impulse  response  of  the  nearly  unaberrated  system 
agreed  very  well  with  the  theoretical  predictions. 

In  the  following  sections,  the  theory  will  be  tested  under  a  wider  range  of 
conditions  Aberrations  will  be  introduced  into  the  optical  system  and  imaging  of  a 
point,  two  closely-spaced  points,  an  edge,  and  a  slit  will  be  performed.  It  will  be  seen 
that  in  these  cases  also  the  theory  is  sufficient. 

5.2.1  Aberrated  Impulse  Response  Data 

•  Aberrations  were  introduced  by  a  number  of  methods.  Spherical  aberration  was 
infoouced,  for  instance,  by  reversing  the  orientation  of  lens  L2  in  the  configurations  of 
rigs  5  1  and  5  2  The  lens  L2  was  usually  a  well-corrected  telescope  doublet  Such 
e'ements  are  designed  assuming  that  one  of  ‘ne  conjugates  is  at  nfinitv  ~his  'esuits  n 
the  surface  of  the  doublet  having  the  largest  curvature  being  on  the  infinite  conjugate 
side  For  example,  the  lenses  in  Fig.  3.1  are  drawn  with  the  infinite  conjugate  sides 
‘owards  each  other  In  this  orientation  ’he  system  will  have  a  minimum  amount  of 
spherical  aberration.  This  is  indeed  how  the  system  was  oriented  when  the  data 
tabulated  in  Table  5.2  were  gathered  if  the  orientation  of  lens  L2  were  reversed  so  ‘hat 


the  infinite  conjugate  side  faced  the  image  plane,  it  could  be  expected  that  the 
dominant  aberration  in  the  system  would  be  third  order  spherical. 

This  is  indeed  the  case  for  the  data  in  Figs.  5.8  through  5.11.  The  data  in  these 
figures  were  collected  using  the  configuration  of  Fig.  5.1,  where  different  lenses, 
oriented  for  maximum  spherical  aberration,  were  used  as  L2. 

In  Figs.  5.8  and  5.9,  photographs  of  the  irradiance  impulse  response  of  an 
aberrated  system  are  displayed  along  with  calculated  3-D  plots  of  the  modulus  of  the 
amplitude  impulse  response  generated  by  the  model  PSF.FOR  using  the  measured 
aoe'-rations  of  the  photographed  system.  Fig.  5.8  is  for  an  unapodised  system  and  Fig.  5  9 
IS  for  an  apodised  (G  =  3)  system.  The  measured  third  order  5eidel  coefficients  were 
astigmatism  =  0.1  X,  coma  =  O.IX,  and  spherical  =  -0.8X.  There  were  negligible 
amounts  of  higher  order  aberrations.  In  the  theoretical  plots  shown  just  below  the 
photographs,  the  modulus  is  used  instead  of  the  irradiance  to  emphasize  the  ring 
structure. 

The  agreement  between  the  experimental  and  theoretical  data  is  qualitatively 
ve'>  good  m  the  case  of  the  unapodised  impulse  response  (Fig.  5  8).  tne  mag-'it.^ae  o' 
t.’'e  •irsi  ri.ng  grew  whiie  the  magnitudes  of  the  central  peais  ana  tne  seccna  "i.ng 
decreased  relative  to  the  unaberrated  case  (compare  Fig.  5.8  to  Fig.  5.4)  This  behavior 
was  preaicteo  m  Figs.  3. 1 1  and  3. 15.  As  seen  in  Fig.  5.9,  the  apooiser  was  quite  effective 
in  removing  the  ring  structure.  The  imoulse  response  has  been  considerably  smoothed 
wnicp  s  m  agreement  with  theory.  This  smoothing  implies  that  tne  apodiser  will  be 
e”’ect:'.e  n  reducing  the  ringing  in  the  image  of  an  edge  or  a  slit,  as  -veil  as  reoucmg  :‘',e 
mensuration  error  m  the  image  of  two  closely-spaced  oomts.  Data  from  these  types  of 
images  through  a  system  with  the  present  impulse  response  will  be  displayed  in  the  next 


section. 


When  other  lenses  of  various  focal  lengths  and  orientations  were  substituted  for 
L2,  different  impulse  responses  were  obtained.  Data  collected  by  the  electro-optic 
technique  of  Section  5. 1.4.1  are  shown  in  Figs.  5  10  and  5.11.  In  each  figure  the 
theoretical  data  are  plotted  as  solid  lines  and  the  experimental  data  are  plotted  as 
dashed  lines.  The  unapodised  (GsO)  cases  are  on  the  left  side  of  each  figure  and  the 
apodised  (G  =  3)  cases  are  on  the  right.  The  measured  amounts  of  third  order  spherical 
aberration  for  the  four  cases  was  -0.2X,  -0.6X,  -1.2X,  and  -1.9X  for  Figs.  5.10a,  5  10b, 
5  11a,  and  5. 1  lb  respectively.  Smaller  amounts  of  the  other  third  order  aberrations  were 
present  and  are  noted  in  the  figure  captions.  Higher  order  aberrations  were  negligible 

These  data  show  quantitatively  the  same  effects  seen  qualitatively  in  Figs  5  8  and 
5  9  There  is,  first  of  all,  good  agreement  between  the  experimental  and  theoretical  data 
in  each  of  the  figures.  The  only  exception  is  that  the  theoretical  and  experimental  data 
differ  by  about  10%  in  some  portions  of  Fig,  5.1  lb. 

Aberrations  other  than  spherical  could  be  generated  by  misaligning  an  originally 
unaberrated  system.  For  instance,  the  optical  system  of  Fig,  5,1  with  the  elements 
oriented  for  minimum  spherical  aberration  was  deliberately  misaligned  until  the  imouise 
response  in  the  top  of  Fig.  5.12  was  obtained  The  aberrations  were  measurec  and  ‘ounc 
’0  be:  astigmatism  =  0.2X,  coma  =  0.2X.  and  soherical  =  -O.'X  "^he  3-D  olot  n  ‘he 
oottom  of  this  figure  is  the  calculated  modulus  of  me  amplitude  impulse  response  of  a 
ivstem  having  the  measured  aoerrations  ana  no  oooaiser  Similar  aata  ‘or  ‘ne  :a.me 
svstem  with  an  apodiser  (G  =  3)  as  part  of  the  system  is  shown  m  Fig.  5  1 3 
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Fig  5.10  Theoretical  (solid  lines)  and  experimental  (dashed  lines)  results  are  compared 
for  two  systems.  In  the  first  system  (a)  astigmatism  s  O.lX.coma  s  0.  IX,  and 
spherical  =  -0.2X;  and  in  the  second  system  (b)  astigmatism  =  O.lX.coma  = 
0.1a,  spherical  =  -0.5X.  The  piots  on  the  left  are  for  the  unapooiseo  case 
while  the  others  are  for  the  apodised  (G  =  3)  case. 
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Fig  5  1 1  Theoretical  (solid  lines)  and  experimental  (dashed  lines)  results  are  compared  ./-V 

for  two  systems.  In  the  first  system  (a)  astigmatism  =  0. OX,  coma  =  0  lX.and  jT 

soherical  =  -1.2X;  and  in  the  second  system  (b)  astigmatism  =  0.2X,coma  = 

0.3 X,  spnerical  =  -1.3X.  The  piots  on  tne  eri  are  for  the  unapooiseo  ;ase  v'/.- 

while  the  others  are  for  the  aoodised  (G  =  3)  case.  v*;- 


The  pnotograpn  is  of  the  irradiance  impulses  response  of  an  unapooiseO 
optical  system  with!  he  aberrations  of  astigmatism  s  0.2X,coma  =  0.2X.and 
spherical  =  >0.1  X.  The  3-D  plot  is  the  theoretically  predicted  modulus 
(vertical  axis)  based  on  these  aberrations.  The  horizontal  axes  are  'n  the  same 
units  of  normalized  distance. 
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f-ig  5  13  The  photograph  is  of  the  irradiance  impulses  response  of  an  apodised  (G  =  3) 
optical  system  witht  he  aberrations  of  astigmatism  s  0.2X,coma  =  0.2X,and 
soherical  =  -0.1  X.  The  3-0  olot  is  the  theoretically  predicted  modulus 
(veatcal  axis)  based  on  these  aoerrations.  The  horizontal  axes  are  m  the  same 
units  of  normalized  distance. 
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A  quantitative  comparison  between  theory  and  experiment  for  the  impulse 
responses  of  Figs.  5.12  and  5.13  is  provided  by  the  data  shown  in  Figs.  5.14  and  5  15.  The 
data  in  Fig.  5. 14  were  generated  by  placing  the  l-SCAN  linear  detector  array  in  the  center 
of  t^-?  -^pulse  response.  The  orientation  of  the  linear  array  corresponds  to  a  horizontal 
.  rough  the  center  of  the  photographs  of  Fig.  5.12  and  5.13,  with  the  sense  of  left 
and  right  reversed.  The  analog  output  of  the  l-SCAN  detector  array  was  displayed  on  the 
CRT  of  an  oscilloscope  where  it  was  photographed. 

These  photographs  are  displayed  m  Fig.  5.14  The  horizontal  axes  are  the  same  for 
both  photographs  and  correspond  to  a  linear  distance  in  the  plane  of  the  impulse 
response.  The  veritcal  axes,  respresenting  image  irradiance,  also  have  the  same  units, 
except  that  for  the  top  photograph  a  neutral  density  filter  (ND  =  0.7)  was  inserted  into 
the  system. 

The  data  of  Fig.  5.14  are  displayed  as  black  dots  in  the  plots  of  Fig.  5.15.  The  solid 
curves  in  this  figure  are  theoretical  plots  of  the  irradiance  along  a  line  through  the  center 
of  the  impulse  resoonse  for  a  system  having  the  measured  aberrations  of  Fig  5  12.  The 
too  olots  ir  both  =ig  5  14  and  5. '  5  are  for  an  unapodised  system,  while  the  bottom  plots 
are  ‘or  an  apodised  IG  =  3)  system. 

From  the  several  different  imoulse  responses  considered  in  this  section  t  is  dear 
that  the  model  is  accurate  m  predicting  the  experimental  results.  The  largest  error  is 
3Dout  '0°/a  m  the  case  wnere  the  cpnencal  aoerration  eaualeo  ‘  9X.  n  the  otner  cases, 
involving  aberration  with  magnitudes  less  than  this,  the  agreement  between  theory  anc 
e<cenme'''  was  ess  than  5'-i 
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Fig.  5. 14  Exp«rtm«ntal  data  through  tha  cantor  of  tha  unapodisad  (top)  and  apodisad 
(bottom)  irradianca  impulsa  rasponsa  of  a  systam  with  astigmatism  «  0.2X, 
coma  *  0.2X,andspharical  «  -0.1  X.  Tha  ap^isar  had  a  valua  of  G  >  3. 


S.2.2  Abtrrattd  Imaging  Data 

Imaging  experiments  were  performed  where  the  objects  were  an  edge,  a  slit,  and 
a  set  of  two  closely-spaced  points.  The  experimental  configuration  is  shown  in  Fig.  S.2. 
The  specific  optical  system  used  here  was  the  one  which  generated  the  impulse  response 
data  of  Figs.  5.8  and  5.9.  That  is,  the  measured  aberrations  were:  astigmatism  s  O.lX, 
coma  =  0.1X,  and  spherical  s  -O.SX. 

5. 2. 2.1  Edge  and  Slit  Images 

After  the  aberrations  were  measured  and  the  impulse  data  were  taken,  the 
pinhole  in  the  spatial  filter  was  replaced  with  the  edge  object.  Both  the  pinhole  and  the 
edge  object  were  on  identical  magnetic  mounts  so  that  the  edge  object  went  into  the 
system  at  the  same  axial  location  formerly  occupied  by  the  pinhole.  The  10X  microscope 
objective  was  moved  axially  back  from  the  pinhole  so  as  to  provide  a  relatively  uniform 
illumination  at  the  edge  object.  The  10  micron  pinhole  was  reinserted  at  the  focal  point 
of  the  laser  beam  after  the  objective.  This  was  done  to  spatially  filter  out  the  high 
frequency  noise  on  the  beam.  The  edge  was  a  section  of  a  razor  blade. 

Data  were  collected  by  placing  the  t-SCAN  detector  array  in  the  image  plane  with 
the  axis  of  the  detector  array  perpendicular  to  the  edge  image.  The  analog  output  from 
the  detector  array  was  connected  to  an  oscilloscope  and  photographs  were  taken  of  the 
patterns  displayed  on  the  CRT. 

Data  obtained  in  this  manner  are  shown  in  Fig.  5. 16  where  the  coherent  image  of 
an  edge  with  no  apodisation  is  displayed  in  the  top  photograph  and  the  same  image 
with  an  apodiser  (G  s  3)  in  the  system  is  displayed  below  it. 

These  data  (in  the  form  of  black  dots)  are  shown  superimposed  on  theoretical 
curves  (solid  lines)  in  the  plots  of  Fig.  5. 17.  Two  conclusions  can  be  drawn  from  these 


Fig.  S.16  Experimental  data  from  the  analog  output  of  the  USCAN  linear  detector 
oriented  perpendicular  to  the  edge  image  through  an  unapodised  (top  plot) 
ar>d  apodised  (bottom  plot,  G  «  3)  optical  system  with  measured  third  order 
aberrations:  astigmatism  «  O.lX.coma  •  0. IX. and  spherical  «  -0.8X. 


data.  First,  the  apodiser  has  been  effective  in  removing  the  edge  ringing  in  the  image  of 
the  edge  through  this  aberrated  system.  This  behavior  is  a  direct  result  of  the  smoothing 
of  the  impulse  response  seen  when  comparing  Figs.  5.8  and  5.9.  The  second  conclusion 
that  can  be  drawn  is  that  the  theoretical  model  is  quite  accurate  in  predicting  the 
experimentally  observed  behavior. 

When  the  object  was  two  razor  blade  edges  instead  of  one,  a  slit  was  formed.  The 
width  of  the  slit  was  varied  and  image  data  were  collected  in  the  same  manner  as  with 
the  edge.  The  photographic  data  for  three  different  slit  widths  are  shown  in  Figs.  5.18, 
5.20,  and  5  22  These  data  (in  the  form  of  black  dots)  are  shown  relative  to  theoretical 
predictions  in  Figs.  5  19,  5.21,  and  5  23  respectively.  The  predictions  were  based  on  the 
measured  system  aberrations  (same  as  for  the  edge). 

The  same  conclusions  that  were  drawn  for  the  edge  can  be  drawn  here.  Namely, 
the  apodiser  is  effective  in  removing  the  edge  ringing  and  t^'e  theoretical  model  does  a 
good  job  in  predicting  the  experimental  outcome. 

The  small  deviations  between  theory  and  experiment  in  the  top  plots  of  Figs.  5  19 


and  5.20  are  caused  by  speckle.  It  is  interesting  to  note  that  the  apodiser  has  diminished 
this  difference  (compare  the  top  and  bottom  plots  of  Figs  5  19  and  5.20) 


Fig.  S.20  Experimental  data  from  the  analog  output  of  the  l-SCAN  linear  detector 
oriented  perpendicular  to  the  slit  image  through  an  unapodised  (top  plot) 
and  apodised  (bottom  plot,  G  s  3)  optical  system  with  measured  third  order 
aberrations:  astigmatism  b  O.lX.coma  x  0. 1 X,  and  spherical  x  •0.8X. 


Fig.  5.22  Exp«rim*ntal  data  from  th*  analog  output  of  the  l-SCAN  linear  detector 
oriented  perpendicular  to  the  slit  image  through  an  unapodised  (top  plot) 
and  apodised  (bottom  plot,  C  *  3)  optical  system  with  measured  third  order 
aberrations:  astigmatism  «  O.lX.coma  «  0. 1 X.  and  spherical  =  -0.8X. 


5.2.2.2  Two-point  Images 

The  data  collected  in  the  two-point  imaging  experiments  were  from  two  different 


optical  configurations.  The  first  configuration  (Fig.  5.1  was  used)  was  designed  to  have 
small  amounts  of  aberrations.  A  measurement  of  the  aberraticn-  ?.^ter  this  design  was 
configured  gave:  astigmatism  =  0.1  X,  coma  =  O.OX,  and  =  O.OX. 

The  five  double-pinholes  described  in  Section  5. 1.2.1  were  individually  inserted 
into  the  system  and  the  resulting  image  plane  separation  of  peaks  was  measured  An 
example  of  the  data  obtained  from  this  process  is  shown  in  Fig.  5.24.  In  this  case  the 
object  was  the  double  pinhole  with  a  pinhole  separation  of  42  microns.  In  the  top  plot  of 
this  figure,  the  experimental  data  with  no  apodiser  was  drawn  as  a  dashed  line  The 
theoretically  predicted  image  irradiance  distribution  for  the  same  condition  was  drawn 
as  a  solid  line.  When  the  apodiser  (G  =  3)  was  added  to  the  system  the  curves  in  the 
bottom  plot  of  Fig.  5.24  are  obtained.  The  geometrically  expected  locations  of  the 
pinhole  images  are  indicated  by  vertical  dashed  lines. 

This  procedure  was  repeated  for  the  other  double  pinholes  as  well  as  for  all  five 
double  pinholes  in  an  essentially  unaberrated  system  having  twice  the  magnification  All 
of  these  data  are  combined  into  Fig.  5.25  The  horizontal  axis  in  this  figure  is  the 
geometrically  exoected  image  ooint  separation  while  the  vertical  axis  is  the  mensuration 
error.  The  mensuration  error  is  the  difference  between  the  expected  and  the  actual 
seoaration.  "he  ’heoreticallv  oredicted  mensuration  e'^ror  ;$  snown  as  a  solid  ‘me  ‘or  tne 
unapodised  case  and  as  a  dashed  line  for  the  apodised  (G  =  3)  case.  The  experimental 
data  points  are  shown  as  boxes  for  the  jnaoodised  situation  and  as  X's  for  the  apodised 
(Q  -3)  situation., 

A  similar  set  of  data  is  displayed  in  Fig.  5.26.  In  this  case  the  optical  system  was 
modeled  after  Fig.  5  2  and  the  measured  aberrations  were  astigmatism  =  0.1  V  coma  = 
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The  irradiance  distribution  in  the  image  of  a  two-point  object  through  an 
essentially  unaberrated  system  which  is  unapodised  in  the  top  plot  and 
apodised  (G  =  3)  in  the  bottom  plot.  In  each  plot,  the  solid  curve  represents  a 
theoretical  preoiaion  ana  the  oroKen  curve  represents  experimental  aata. 


151 


O.lX,  and  spherical  =  -O.SX. 

In  both  Figs.  5.25  and  5.26  the  data  show  that  a  mensuration  error  exists  when  the 
separation  is  close  to  the  limit  of  resolution.  This  is  true  for  both  cases  of  apodisation  and 
for  both  figures.  However,  the  data  also  reveal  that  the  mensuration  error  when  the 
points  are  resolved  is  less  when  the  sytem  is  apodised.  The  unapodised  case  has  a  better 
resolution  limit,  but  this  is  misleading  because  the  mensuration  error  is  greatest  in  the 
region  between  the  resolution  limits  for  the  two  cases  of  apodisation.  This  conclusion 
applies  equally  to  the  unaberrated  and  aberrated  situations. 

5.3  EXPERIMENTAL  SUMMARY 

Experiments  were  performed  involving  the  imaging  of  simple  objects  through 
coherent  optical  systems  having  various  amounts  of  optical  aberrations.  The  effect  on 
imaging  performance  of  adding  a  Gaussian  (G  =  3)  apodiser  to  the  optical  system  was 
determined.  The  objects  were  a  point,  two  closely-spaced  points,  an  edge,  and  a  slit. 

The  goal  of  these  experiments  was  to  adequately  test  the  theory  develooed  m 
Chapters  2,  3,  and  4  These  chapters  contain  a  large  amou.nt  of  theoretical  data  'he 
experiments  were  not  intended  to  test  eacn  individual  theoret.cai  result:  rather 
representative  theoretical  results  were  tested.  The  experimental  results  were  very  close 
to  the  theoretical  predictions  for  those  cases  which  were  tested.  The  experiments 
confirm  that  'he  theorv  is  valid  at  least  under  the  range  of  conditions  tested. 


CHAPTERS  CONCLUSIONS 


This  study  was  performed  to  investigate  both  theoretically  and  experimentally  the 
effects  of  aberrations  on  the  performance  of  coherent  imaging  systems,  and  the  utility  of 
apodisation  in  improving  that  performance.  It  has  been  established  by  other  workers 
that  apodisation  is  useful  in  improving  the  performance  of  unaberrated  coherent 
imaging  systems.  Any  real  optical  system,  of  course,  has  some  amount  of  residual 
aberrations.  So  this  study,  by  considering  the  reality  of  aberrations,  has  extended  the 
body  of  knowledge  in  the  field  of  apodisation 

It  IS  worth  emphasizing  the  ubiquity  of  aberrations  in  optical  systems.  The.'e  are  a 
few  optical  systems  which  are  essentially  free  of  aberrations;  for  example,  the  projection 
system  used  to  make  masks  for  the  semiconductor  industry.  Such  systems  are  extremely 
expensive  however.  It  is  more  likely  that  an  optical  system  will  be  considered  well 
corrected  if  it  meets  the  Rayleigh  quarter  wavelength  rule.  This  rule  of  thumb  states  that 
an  element  is  well  corrected  if  the  absolute  value  of  its  dominant  aberration  does  not 
exceed  a  quarter  of  a  wavelength  In  buying  ootical  elements  for  this  experiment,  optics 
of  this  quality  were  found  to  be  the  norm,  initially  in  the  experimental  phase  it  was  very 
difficult  to  acheive  a  system  having  low  amounts  of  aberrations  (say  less  than  0  1  waves 
of  each  of  the  third  order  aberrations.)  The  problem  was  that  each  of  the  individual 
elements  ;ust  met  'he  quarter  wavelength  ruie.  Eventually  high  quality,  large  ’14 
telescope  objectives  were  employed  to  acheive  a  low  aberration  system  But  even  a 
system  like  '.his  could  have  seve'^al  tentrs  of  a  wave  o'  aoemation  un’ess  *.ne  e‘e'~e'"'s 
.ve'e  aligned  -virh  great  'are  ’’I'e  'esson  ne'e  s  *".at  ■  s  /e'’-  d:“lcjlt  ‘o  ac*'e''-e  'm 
effeaively  Jnaberrated  system.  This  difficulty  translates  into  extra  weight,  size,  and 
expense  in  any  oractical  ootical  system 
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In  the  theoretical  development  of  Chapter  2  and  its  application  to  specific  imaging 
situations  in  Chapters  3  and  4,  several  simplifying  assumptions  were  made.  These  were 
clearly  noted  when  made  but  for  convenience  are  asserribled  as  a  group  here.  The  major 


assumptions  are  .  ..^^g. 

1 .  Scaia;  r  ^s.'.ei-Kirchoff  diffraction  theory  (see  for  example  Born  and  Wolf^)  is 
adequate  to  describe  the  systems  investigated  here. 

2.  The  system  parameters  are  limited  such  that  the  complex  amplitude  Uu  in  the 
(xu..Vo)  plane  is  related  to  the  complex  amplitude  in  the  (z,,y,)  plane  by  the  Fresnel 
propagation  formula 


.  2n 

‘V  ‘  i\z  J  J_.  n  "  0/ 


dx  dy  , 
0  "^0 


where  X  wavelength  of  light  and  z  is  the  distance  between  planes.  This  equation  will  be 
valid  if  the  largest  aperture  dimension  is  much  greater  than  X  and  less  than  z/10. 

3.  The  light  source  is  monochromatic  and  produces  a  field  in  the  first  plane  of 
interest  which  is  completely  spatially  coherent 

4  The  imaging  process  for  the  magnitude  of  aberrations  considered  and  the 
size  of  the  objects  is  isoplanatic.  There  are  aberrations  m  the  optical  systems  considered 
n  this  Study,  so  the  ‘orm  of  the  imoulse  resoonse  changes  as  an  object  ooint  exciores  'he 
object  plane.  However,  for  the  small  region  of  the  objects  considered  here  it  was 
assumed  that  the  form  of  the  impulse  resoonse  changes  negiigibiv  n  this  tase.  the 
imaging  process  was  expressed  as  the  convolution  integral 

u  ( z  O'. )  =  I  f  ^  .-y  ~y^)d^.dy,^ . 


where  K  is  the  amplitude  impulse  response  of  the  system. 


The  assumption  of  isoplanitism  for  the  systems  considered  in  this  study  were  checked 
and  found  to  be  valid.  The  largest  object  considered  was  a  slit  which  had  an  image  width 
of  V  =  80.  This  corresponds  to  an  object  width  of  approximately  300  microns  for  the 


imac:  -  n  of  Fig.  5.2.  Theoretical  imaging  calculations  were  performed  using  this 
ODjeci  mzs  through  an  optical  system  having  as  much  as  four  waves  of  defocus  or  one  of 
the  third  order  aberrations. 

The  isoplanatic  assumption  implies  that  an  off-axis  amplitude  impulse  response  in  the 
vicinity  of  the  edge  of  the  largest  slit  differs  only  slightly  from  an  on-axis  impulse 
response.  These  two  amplitude  impulse  response  functions  were  calculated  using  the 
program  PSF.FOR  described  in  Chapter  3.  The  off-axis  amplitude  impulse  response  was 
then  translated  until  i  peak  was  aligned  with  the  peak  of  the  on-axis  impulse  response. 
The  difference  in  the  complex  values  of  the  two  patterns  was  obtained  and  the  intensity 
of  this  difference  was  found.  The  area  under  the  resulting  intensity  pattern  was 
calculated  and  compared  to  the  area  under  the  on-axis  intensity  impulse  response.  The 
comparison  was  in  the  form  of  a  ratio.  If  this  ratio  was  less  than  0.01  then  the  assumption 
of  isoplanatism  was  considered  justified. 

This  calculation  was  performed  for  the  off-axis  point  being  150  microns  from  the  axis 
and  with  4  0  waves  of  one  of  the  aber'ations  'deiocus.  astigmatism,  tom.a.  and  sche'”';3!! 
A  different  calculation  for  each  aberraton  was  performed.  The  effect  of  apodisation  was 
also  determined,  .n  each  case,  the  'atio  was  .much  less  man  0.01 

Since  the  agreement  between  the  experiments  and  the  theory  is  so  good,  the 
assumptions  outlined  here  clearly  are  valid  for  the  optical  systems  considereo  .n  this 
•■'“SIS.  n  -rat  .ase,  remain  general  -ronciusicns  can  he  drawn  -irst  of  aii.  the  cresense  of 
aberrations  in  a  coherent  imaging  system  degrades  the  imaging  performance  of  that 
system  m  a  oredictaole 'way  ana  apocisatio''  :s  effecf  ve  n  imorovinq  selected  ascec's  o’ 


system  performance.  The  experiments  conduaed  on  both  the  unapodised  and  apodised 
systems  agree  well  with  the  theoretical  predictions.  This  agreement  implies  that  the 
other  theoretical  predictions  of  Chapters  3  and  4  are  valid  as  well. 

It  can  also  be  concluded  that  the  use  of  apodisation  is  beneficial  in  improving  certain 
aspects  of  system  performance.  Specifically,  apodisation  of  aberrated  coherent  imaging 
systems  reduces  the  mensuration  error  in  the  image  of  two  closely-spaced  points  and 
reduces  the  edge  ringing  in  the  images  of  both  an  edge  and  a  slit. 

It  appears  that  expressing  the  images  in  terms  of  the  square  root  of  irradiance  is 
beneficial.  This  effectively  removes  the  nonlinearity  associated  with  the  detect. on 
process  As  a  result,  the  edge  ringing  and  shift  are  reduced.  If  the  asymptotic  value  of 
the  irradiance  on  the  bright  side  of  the  edge  can  be  determined,  then  the  edge  shift  can 
be  reduced  to  zero,  except  when  unsymmetrical  aberrations  are  present. 

Finally,  this  research  has  indicated  that  there  is  still  some  work  to  be  done  m  order  to 
fully  understand  the  effect  of  apodisation  on  aberrated  coherent  optical  systems 
Specifically,  the  effect  of  aberrations  and  the  utility  of  apodisation  on  beam  prooagat 
systems  should  be  studied.  The  theoretical  effect  of  various  aoodisers  on  aoe^'a;  o-.'-oe 
beam  propagation  systems  has  been  studied,  but  the  benefit  of  using  aooc  so;  on  ~ 
aberrated  beam  orcoaoation  systems  has  .not  been,  mvestioated. 
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APPENDIX  1  COMPUTER  PROGRAMS 


The  computer  programs  which  are  documented  in  this  appendix  were  developed 
to  perform  the  calculations  leading  to  the  various  theoretical  resulu  of  Chapters  3  and  4. 
The  program  PSF.FOR  is  the  parent  program;  all  the  others  are  derived  from  it. 

The  program  PSF.FOR  calculates  the  modulus,  phase,  and  irradiance  of  the  point 
spread  function  of  an  optical  system.  The  inputs  to  the  program  are  the  size  of  the  array 
representing  the  exit  pupil,  the  third  order  aberrations  in  terms  of  Zernike  coefficients, 
and  the  radial  width  factor  of  the  Gaussian  apodiser.  The  point  spread  function  is  found 
by  performing  a  fast  Fourier  transform  of  the  exit  pupil  array. 

The  program  TPNT.FOR  calculates  the  irradiance  in  the  image  of  two  point  sources. 
The  point  sources  are  of  equal  amplitude  and  have  zero  relative  phase.  TPNT.FOR  is  an 
only  slightly  modified  version  of  PSF.FOR.  The  modification  is  an  additional 
multiplicative  term  in  the  exit  pupil: 


ARG2=F 


6  z 


„  f  2niS’  \ 

=  2  cos  -  , 

'  V  > 


(A'  l) 


where  b  is  the  separation  of  each  object  point  from  the  optical  axis,  y  is  the  linear 
dimension  of  the  exit  pupil  parallel  to  a  line  joining  the  two  points,  is  the  wavelength 
of  the  illumination,  and  f  \%  the  distance  from  the  object  plane  to  the  entrance  pupil  The 
modification  occurs  m  the  section  of  '.ne  crogram  iaoeled  ’Circular  Aoerture.  "  Vluch 
of  tne  program  after  that  point  is  not  included  because  it  is  identical  to  PSF  FOR 

The  programs  EDGE. 'OR  and  SLIT  FOR  are  also  modifications  of  PSF  FOR  Both 
orograms  have  been  soecialized  to  one  dimension  because  of  the  inherent  one- 
dimensional  nature  of  the  edge  and  slit  objects. 
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PROGRAM  PSr 


C 

C 

C  THIS  PROGRAM  CALCULATES  THE  PHASEr  MODULUS  AND 
C  INTENSITY  OF  THE  TWO-DIMENSIONAL  IMAGE  OF  AN 
C  ON-AXIS  POINT  SOURCE  THROUGH  AN  OPTICAL 
C  SYSTEM  WITH  ARBITRARY  ABERRATIONS  AND 
C  APOniSATION 
C 
C 

INTEGER  IWK(2000>  rIAl » I A2 r N1 » N2 r N3 » IJOB 
REAL  RWK<2000) 

REAL  P(48f4B)»M(48f4B)>E(4B)»INT(48f48) 
COMPLEX  A ( 256 » 256) »CUK(256) » D (48^48) »C(48) 
COMPLEX  D<4B) 

C 

aPEN(UNIT=22.STATUS='NEU' .FILE='PH.DAT  '  ) 
0PEN(UNIT=23fSTATUS='N£W' » F 1LE= ' MOD . DAT ' ) 

OPEN ( UNI T=24fSTATUS=' NEW' i FILE= ' INT . DAT ' ) 

OPEN (UNIT=25»STATUS=' NEW ' fFILE='BB. DAT'  ) 

C 

C  READ  PARAMETER  UALUES  FROM  MIC  FILE 
C 

READ(5»20)N 
READ<5f40)A20 
READ(5»40)A21 
READ(5f 40)A22 
READ(5»40)A31 
READ<5f40)A32 
READ(5»40)A42 
READ(5f 40)GG 
20  F0RMAT<I4) 

40  FaRMAT(F10,4) 

41F0RMAT( I3»3X>I3f3XfF10,4»3X7F10.4) 

42F0RMAT< I3»3X» I3r3XrF10.4) 

43FGRMAT':3(2XrI3)  .5(3X;F3.4;  , 

44F0RMAT(I3f 3XFF10.4»3XfF10. 4) 

C 

C  N=NUMBER  OF  SAMPLE  POINTS  IN  ONE  DIMENSION 

C  A20=45  des  ASTIGMATISM 
C  A21=F0CAL  SHIFT 

C  A22=0  dea  ASTIGMATISM 

C  A31=X  COMA 

C  A32=Y  COMA 

C  A42=SPHERICAL 

0  GG^'WIDTH  OF  GAUSSIAN 


nnno  nnn  oonn 
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CIRCULAR  APERTURE  WITH  VARIABLE  TRANSMITTANCE 
AND  ABERRATIONS 

DO  220  1=1 »N 
DO  200  J=lrN 

RAD=SQRT( (I-<l8.+0.5) >**2+(J-(18,+0.5) )**2) 
X=<I-(18.+0.5) )/17.0 
Y=( J-(18.+0.5) )/17.0 
CALCULATE  CONTRIBUTIONS  OF  THIRD  ORDER- 
ABERRATIONS 

USING  ZERNIKE  MONOMIAL  REPRESENTATION 
ARG20=A20»2.*X»Y 
ARG21=A21* (-1 .+2.»Y*y+2.»X*X )/2. 

ARG22=A22*( Y*Y-X»X) 

ARG31=A31»(-2.*X+3.*X*Y*Y+3.*X*»3)/3. 

ARG32  =  A32*(-2.»Yt3.*Y**3  +  3  .  ♦X!t:X*Y  )/3. 

ARG4  2A=A42*(  1  . -6  .  »  Y*Y-6  .  *X»X  +  6  .  ♦  Y!*!*4  ) 
ARG42=ARG42A+( 12.*X*X*Y»Y+6.»X»»4 )*A42 
ARG42=ARG42/6. 

ARGl  =  (  ARG20  +  ARG21  +  ARG22  +  ARG3 1  +ARG32  +  ARG42  )  *2  .  *F- 1 
CALCULATE  CONTRIBUTIONS  OF  VARIOUS  APODISERS 
GAUSSIAN  APODISER 
TT=EXP<-GG*RAD»RAD/17.0/17,0) 

SUPERGAUSSIAN  APODISER 

F  =  1  . 

TT  =  EXF  (  -  (  GG*RAD*F<AD/ 1  7  .  /  I  7  .  )  **F  ) 
ir (RAD.LT, 16.5)G0  TO  60 
IF'CRAD.GT.  17.S)G0  TO  100 
Sr.  =  -RAD+17.5 
GO  TO  80 

60  A(  I  »  J  )  =C,MPLX  (  TT*COS  (  arg:  *SIN  (  arg:  : 

3  C  TO  120 

80  A( I , J)=CMPLX(SC*TT*CCS( ARCl ) » 5C» T T*S I N ( ARG 1 )  i 
GO  TO  120 

1  'j  0  A  ( I ,  j :  =  ( 0 .  f  0 . 

120  CONTINUE 
200  CONTINUE 
220  CONTINUE 
DO  230  1=1 f 48 
DO  225  J=1 f 48 
M  ■:  I  •  J  )  ^CAPS  (A(I»w); 

CWRITE(24f 42)1 r Jf M( I f J) 

225C0NTINUE 

230C0NTINUE 
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C  PERFORM  FAST  FOURIER  TRANSFORM  USING  IMSL 
C 

240  IA1=N 
IA2  =  N 
N1  =  N 
N2  =  N 
N3=l  • 

IJOEi  =  -l 

CALL  FFT3D  ( A p I A1 » I A2 » N1 r N2 » N3 r I  JOB r I WK r RUK » CUK ) 
DO  242  J=1»N 
1  =  1 

CURITE(22f41)I»JfA(If J) 

242C0NTINUE 

C 

C  REDUCE  SIZE  OF  ARRAY 
C 

DO  250  I=lf24 
DO  245  J=1 »24 
B( I f J)=A( I f J) 

&( I »24rJ)=A( I »N-24+J) 

B(24+I , J)=A(N-24+I f J) 

B(24-t  I  f  24  +  J)=A(N-24+I  f  N-  24  +  J> 

245C0NTINUE 
250C0NTINUE 
DO  252  I=l»48 
DO  251  J=l»48 
CUF;ITE(23f41)I»J»E<(Ir  J) 

251C0NTINUE 

252C0NTINUE 

C 

C  SHIFT  ZERO  FREQUENCY  TO  ARRAY  CENTER 
C 

DC  260  1  =  1  »  48 
DO  255  J=1 .24 

A'  I  .24tJ)=&<.  I  .49  —  J) 

255C0NTINUE 
260C0NTINUE 
DO  270  1=1.24 
DO  265  J=1.4e 
B(  I .  J)=A(25-I . J) 

&(  24-rI  .  J  )  =A  (  49-1  .  J  ) 

265C0NTINUE 

270C0NTINUE 

DU  272  1  =  1.48 
C  TYPE  41 . I . J.B( I . J) 

272C0NTINUE 


c 

C  FIND  PHASE  OF  ARRAY 
C 

248  DO  360  I=1»4B 
DO  340  J=1.48 
AA=AIMAG<P( I f J)  ) 

BD=REAL(B(I»J> ) 

IF  (BB,EQ.O.OOOOOOOO)BB=0. 000001 
Pdf  J)=ATAN2(AA.B&) 

C  WRITE(25f42)I» JfPCIf J) 

340  CONTINUE 
360  CONTINUE 
C 

C  FIND  MODULUS  AND  INTENSITY  OF  ARRAY 
C 

DO  400  I=lf4B 
DO  300  J=lf48 
M( I f J)=CABS(&( I f J) )/.0139 
lNTd,J)=M(I,J)*»2. 

380  CONTINUE 
400  CONTINUE 
DO  404  1=1,48 
DO  402  J=l,48 

C  URITE(23f42) I , J,M( I , J) 

URITE(24 ,42) I , J, INT( I , J) 

402C0NTINUE 

404C0NTINUE 

C  WRITE  CENTRAL  SLICE  OF  INTENSITY 
DO  415  1=1,48 
CURITE(23,44)I ,M(I,24) 

WRITE (25,44)1, INT ( 1,24) 

4l5CaNTINUE 

(-• 

u 

C  F'EMOUC  TILT  FROM  EACH  ROW 
C 

r,n  420  1  =  1 

DC  410  0=2,48 

P(  I  ,  J)=F ( I , J)-( J-1 )» . 4  295 

C  URITE(23,42) I , J,F (  I  ,  J) 

JIOCONTINUE 

420C0NTINUE 

C 


C  HANDLE  PHASE  JUMPS  ON  EACH  ROU 
C 

DO  500  I«=l»48 
DEL=0, 

ICNT=0 

DO  490  J=l»48 
ICNT*ICNT+1 
P<I»J)=P<I»J)+DEL 
F1=P<I* J) 

G1=P<I.J+1)+DEL 
H1=ABS(G1-F1 ) 

CWRITE(22»43)I r J»ICNTfDELfFl rGl fHl 
IF(H1.LE.01)G0  TO  490 

IF( (HI .GE.DL) . AND. (G1 .GT.Fl > )G0  TO  450 
IF( (HI .GE.DL) .AND. (Gl.LT.Fl > )G0  TO  460 
IF( ICNT.LT.5)G0  TO  411 
ICNT=4 

C  DO  LINEAR  FIT  TO  LAST  ICNT  PO I  NTS ( &  1 =SLOFE ) 
411IF(ICNT.E0.1 )G0  TO  431 
SXY=0. 

SX  =  0. 

SY  =  0. 

SX2=0. 

DO  421  K=1»ICNT 
X=ICNT+1-K 

SXY=SXY+X)|rP(  If  J+l-K) 

SX=SX+K 

SY=SY+P(I»J+1-K) 

SX2=SX2+K**2. 

421C0NTINUE 

E<1  =  (  SXY--SX*SY/ICN7  )  /  (  5X2-SX*SX/ 1  CNT  ) 

GO  TO  429 

431&1=P(  I  f  J  +  2)-P(  I  f  J+1 

42?Ir'i'-Bl*LE.0.).AND.'.5l.G7.'-l;;GG  TO  430 

IF( (PI  .GT.O. )  .AND.  (G1 .LT.Fl ) )G0  TO  440 

ICNT  =  0 

GO  TO  4?0 

430INCT=0 

GO  TO  450 

440ICNT=0 

30  TO  460 

450D['L  =  DEL-2.*PI 

GO  TO  490 

460DEL=DEL+2 . *PI 

490C0NTINUE 

500C0NTINUE 

TO  500  1=1.40 

DO  849  J=l»48 

C  lJRITE(24f  42.(  I  f  J»P(  I  »  J) 

849C0NTINUE 

350C3NT:NUE 


C  CORRECT  PHASE  ON  ENTIRE  ROUS  RELATIVE  TO  ROU  1=1 
C 

DO  1160  I«l»48 
J=24 

C  URITE<25»42)I» J»P(I » J> 

1160C0NTINUE 

C  REMOVE  TILT  BETWEEN  ROWS 

DO  1200  I»2»48 

DO  1150  J=lf48 

P( I . J)=P( I » J)-( I-l )».4295 

1150C0NTINUE 

1200C0NTINUE 

[lO  1210  I  =  lf4a 

J=24 

C  URITE(22.42) I f J»P( I f J) 

1210C0NTINUE 

C 

C  HANDLE  PHASE  JUMPS  BETWEEN  EACH  ROW 
C 

DEL=0. 

ICNT=0 

DU  1500  I*lf48 
ICNT=ICNT+1 
DO  1510  J=1 »48 
P( I » J)=P< I , J)+DEL 
1510C0NTINUE 
J  =  24 

F1=P( I , J) 

G1=F(  IPl  jJ)-“EiEL 
H1=ABS<G1-F1 ) 

CWRITE (24  f  43) I F Jr ICNTr  DEL»F1 r G1 rHl 
IF(H1 .LE .01 )G0  TO  1500 

IF.'.  Hl.GE.E■L/>.^i^JI..3i.G~.■~l/)G0  TO  1^50 
IF< (HI .GE.DL) .AND. (G1 .LT.Fl ) )G0  TO  1460 
IF( ICNT.LT.4)G0  TO  1411 
ICNT=3 

C  DO  LINEAR  FIT  TO  LAST  ICNT  POINTS 
1411IF(ICNT.E0.1)GO  TO  1426 
SXY=0. 

S  X  =  0  . 

SY  =  0. 
s  Y  r  ■=  0’ . 

Du  1421  ^=irICN^ 

X=ICNT+1-K 


SXY«SXY+X»P(I+l-Kf J) 

SX=SX+K 

SY=SY+P( I+l-K. J) 

SX2=SX2+K**2. 

1421C0NTINUE 

P1=(SXY-SX*SY/ICNT)/(SX2-SX»SX/ICNT) 

GO  TO  1427 

1426P1=F( 1+2. -  -I+1»J) 

1427  IF(  (B1  .L-w.O.  )  .AND.  (G1  .GT.Fl  >  >GO  TO  1430 

IF(<B1.GT.0.).AND.(G1.LT.F1>>G0  TO  1440 

ICNT=0 

GO  TO  1500 

1430ICNT=0 

GO  TO  1450 

1440ICNT=0 

GO  TQ  1460 

1450DEL=DEL-2.*PI 

GO  TO  1500 

1460riEL=DEL  +  2.*PI 

1500C0NTINUE 

DO  2470  J=l»48 

1=24 

C  URITE<25»42)I» J»P(I» J) 

24  70C'0NTINUE 

C  CENTER  PHASE  AROUND  ZERO  VALUE 
DO  2501  1=1 »48 
DO  2491  J=l»48 
P(I»J)=P(I»J)-9. 

C  URITE(22,42) I f J.P( I . J) 

2491C0NTINUE 

2501C0NTINUE 


C 

:  JF.'ITE  CENTRAL  SLICES  OF  MODULUS  AND  -MASE 

C 

1  =  24 

DO  2550  J=1.48 
C  URITE(22.44) J.P( I . J) 

c  URITE(23f44)J.M(I.J) 

2550C0NTINUE 

CL0SE(UNIT=22) 

CL0SE(UNIT=23) 

CL0SE(UNn=24) 

CLOSE  ;  u'NlT  =  25  ; 

1100  STOP 

END 

C;  IMSL  ROUTINE  NAME  -  FFT3D 


PROGRAM  TPNT 


THIS  PROGRAM  CALCULATES  THE  PHASE r  MODULUS  AND 
INTENSITY  OF  THE  TUO-DIMENBIONAL  IMAGE  OF  TUO 
POINT  SOURCES  THROUGH  AN  OPTICAL  SYSTEM  UITH 
ARDIT^  '-Y  ABERRATIONS  AND  APODISATION 


INTEGER  IWK(2000) tIAlf IA2.N1 »N2»N3rIJ0B 
REAL  RUK(2000) 

REAL  P(48t48) »M(48»48>  fE(48>  > INT(4Bf 48) 
COMPLEX  A(256»256) »CUK(256) »B(48»48) »C(48) 
COMPLEX  D<48) 


0PEN<UNn=22fSTATUS  =  'NEU'  ,FILE='PH.DAT'  ) 
0FEN(UNIT=23f STATUS='NEU' r F ILE= ' MOD . DAT ' > 
0PEN(UNIT  =  24»STATUS='NEU' f F I LE= ' I  NT . DAT ' ) 
aFT.N(UNn=20f  STATUS='NEU'  »  F  ILE=  '  BB  .  DAT  '  ) 

C  READ  parameter  VALUE'S  FROM  MIC  FILE 
C 

READ<5f20)N 
READ<S»40)A20 
READ<5»40)A21 
READ(5f40)A22 
READ(S»40)A31 
READ(5»40)A32 
READOf  40)  A42 
READ(5»40)GG 
RT  ADCSf  40)BBB 

1  L  AD  (  5  *  ^0  )  F!^ 

20  FORMAT (14) 

40  FOFvMAT  (FIO  .  4  ) 

4  1  FORMAT (I3F3XfI3»3X,F10.4,3XfF10.4) 

42F0RMAT( I3f3XfI3F3XfF10.4) 

43F0RMAT( 3<  2X  f 1 3 ) » 5 ( 3X r F8 . 4  )  ) 

44F0RMAT( I3»3X»F10.6f3XrF10,4) 

C 

0  N=NUMBER  OF  SAMPLE  POINTS  IN  ONE  DIMENSION 

0  A20=4S  dea  ASTIGMATISM 

0  A21=F0CAL  SHIFT 

0  A22=0  dea  ASTIGMATISM 

0  A31=X  COMA 

2  A32=^  20MA 

C  A42=SPHERICAL 

0  GG=UIDTH  OF  GAUSSIAN 

0  BBB=OFF-AXIS  DISTANCE  OF  ONE  POINT 

C  Fp=FOCAL  LENGTH  OF  '^IRST  _ENS 

PI=3. 1415927 
ai=3.YPI/4 , 

DL=3.^PI/2. 


•■nnnoo3>  non 


CIRCULAR  APERTURE  WITH  VARIABLE  TRANSMITTANCE 
AND  ABERRATIONS 


DO  220  I=1»N 
DO  200  J=1»N 

RAD=SQRT( ( I-(18.+0.5) )**2+( J-( lB.+0.5> )»*2) 

X=(l-<ia,+0,5) >/17.0 
Y=( J-(18.+0.5> )/17.0 
CALCULATE  CONTRIBUTIONS  OF  THIRD  ORDER- 
ABERRATIONS 

USING  ZERNIKE  MONOMIAL  REPRESENTATION 
ARG20=A20*2.»X»Y 

ARG21=A21»(-1 .+2.*Y*Y+2.*X*X)/2. 

AR-G22  =  A22*(  Y*Y-X»X) 

ARG31=AS1» ( -2.»X+3.*X*Y*Y+3.*X**3)/3. 

ARG32=A32» ( -2 . *Y+3 . *Y**3+3 . ♦X*X*Y ) /3 . 

ARG4  2A  =  A42»  <1.-6. »Y*Y-6 . *X*X  +  6 . *Y**4 ) 

ARG42  =  ARG42A  +  A42»  <12.*X*X»Y*Y  +  6.*X»*4) 

ARG42=ARG42/6. 

ARGl =(ARG20+ARG21+ARG22+ARG31+ARG32+ARG42)»2.#PI 
RG2=2.»PI*B&B!ltY/FF/.6328 

CALCULATE  CONTRIBUTIONS  OF  VARIOUS  APODISERS 
GAUSSIAN  APODISER 
TT=EXP ( -GG*RAD*RAD/17 . 0/17.0) 

SUPERGAUSSIAN  APODISER 

=  1  . 

T=EXP(-(GG!*:RAD*RAD/17./17.  )**F> 
r(RAD.LT. 16.S)G0  TO  60 
r(RAr'.GT.  17.0)G0  TO  100 
:=r-f.:AL'i  I  7 . 0 
C  TO  OC 

■f  C  A  I  1  F  J  )  =  CMPLX  (  TT*COS  (  ARGl  )  »2  .  » COS  (  ARG2  ;  .  TT*SIN  (  ARG  1  )  ) 

GO  TO  120 

80  A(I»J)=CMPLX(SC*TT»C0S(ARG1)*2.!»CQS(ARG2)  .  SC»TT  J*;SI  N  (  ARG 1  ) 
GO  TO  120 
100  A(IfJ/=<0.f0.) 

120  CONTINUE 
200  CONTINUE 
220  CONTINUE 
DO  230  I=lF4e 
DO  225  J=1f4Q 
Y  a  F  J ) =CABS  <  A ( I F  J  )  ) 

CUif.  ITE(24f42!  I  fJfMC  I  F^J; 

225C0NTINUE 

cc:  :  r^us 

c; 

C  PERFORM  FAST  FOURIER  TRANSFORM  USING  IMSL 


PROGRAM  EDGE 


THIS  PROGRAM  CALCULATES  THE  INTENSITY  IN  THE 
COHERENT  IMAGE  OP  AN  EDGE  THROUGH  AN 
ABERRATED  AND  APODISED  SYSTEM 


INTEGER  IrfK(25000) 

REAL  WK(2  500C) ,H(4096) , INT(2O40)  ,E(2C0  ) 
COMPLEX  G(2048) ,F(4096) 

0PEN(UNIT-22,DEVICE-*DSK',ACCESS-'SEQ0UT 

>  PILE* 'PH. DAT’ ) 

0PEN(UNIT-23,DEVICE-’DSK’,aCCESS-*SE00UT 

>  FILE-'MOD. DAT’ ) 

0  PEN (UN  IT-24 , DEVICE- ’DSK ACCESS- 'SEQOUT 
V  PILE-'E. DAT’ ) 

OPEN(UNIT-25,DEVICE-’DSK’ , ACCESS- ' SEOOUT 
1  PILE-’BB.DAT’ ) 

READ  PARAMETER  VALUES  FROM  MIC  FILE 

READ(5 , 2C)N 
READ(5 , 4n)A2C 
RSAD(5 , 4C )a2> 

READ(5 , 4C )a22 
READ(5.4C)a3V 
READ(5 , 4C )a32 
3SAD(5 , 4C )a42 
REaD(5,4C)GG 
2C  F0RMaT(I4) 
ic  ?  :p.hat(  Fi c.  4 ) 

FCRMAT(7^0,4.3X,?>C.4) 

?ORMaT(I3, JX, 13, 3X, FtC. 4  ) 

F0RMAT(3(2X. 13)  .5(3X, P8. 4)) 

FORMAT'.  13.3  X.  7 

FORMAT( >2X, FlC. 6, 3X. F>C. 6  ) 
F0RHAT(I4.3X,PVG.4,3X,F>C.4) 


A2C-45  deg  ASTICMATISM 
A2>" FOCAL  SHIFT 
A22-C  deg  astigmatism 
A3>-X  COMA 
A  3  2  -M  C  OM  A 
A -  S  P  H  E  .  0  A  „ 

GG-rfIDTH  OF  GAUSSIAN 
PI-3.  UV5927 


C  INPUT  FUNCTION  -  AN  EDGE 
C 

DO  6C  I->,VC24 
GU  )-CHPLX(>.  ,0.  ) 

6C  .  CONTINUE 

DO  70  I->025,2048 
G(I )-CMPLX(0. ,C. ) 

CONTINUE 

PERFORM  FIRST  FOURIER  TRANSFORM 


CALL  FFTCC(G,N,  liK.WK) 

DO  9C  I->,2C48 
rfRITE(L5.46)l.G(l) 

CONTINUE 

MOVE  ZEuO  SPATIAL  FREQ  TO  CENTER  OF  APERTURE 

G(255)-G(>2a) 

DO  1-CC  I-1>,>27 
G(>27*I)-G{I) 

G (I )-G (V92t^I  ) 

CONTINUE 
DO  1>C  I->,2C49 
WRITE(22,46  )I ,G(I  ) 

CONTINUE 

V-D  APERTURE  rfITH  VARIABLE  TRANSMITTANCE  AND 
A  B  E  ?  R  A  T  I  0  N  3 


DO  2CC  I-V.255 
Y-(I->23)/>27. 
ARG2>-A2V*(->.*2,*Y*T)/2. 

ARC22-A22*'r *Y  ) 
ARa32-A32'*(-2.*Y»3.*Y**3)/5. 
ARG42-A42*(li.  -6.*Y*Y*6.*Y»*4  )/6. 
ARa>-(ARG2>*ARG22*ARG32*ARG42)*2.*?I 
3  A  03  3  1  AN  AP0DI3EH 

A-EXP(-(GG*Y*Y  )  ) 

P(I)-CMPLX(a*C0S(aRC>).A*3IN(aRGV))*3( 

CCNTINl'E 

DO  22C  1-236, 4C9i 
F( I )-CMPLX(C. ,0. ) 

:  I’iT  INiJE 

DO  23C  I->,4096 
WRITE(24,46  )I ,F(I  ) 

CONTINUE 


3n 
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C  PERFORM  SECOND  FAST  FOURIER  TRANSFORM 
C 

M-4C96 

CALL  FFTCC  (F,M,IrfK,¥K) 

DO  24C  I-»,4096 
C  WRITE(25,46)I,P(I ) 

240  CONTINUE 

C 

C  CALCULATE  INTENSITY  :  USED  TO  NORMALIZE 
C 

DO  244  I-2C0C,2>5C 
H(I )-((CABS(F(l ))/>E*3)**2. )*.23S4 
X -(I -2C5C  )*>3. 85/96 
WRITE(24.4V)X,H(I) 

244  CONTINUE 

C 

C  GEOMETRIC  EDGE 
C 

DO  26C  I -2000,2050 
E(I  )-0. 

260  CONTINUE 

DO  270  I-2C5V,2>5C 
Ed  )•>. 

270  CONTINUE 

DO  230  1-2000, 2>50 
X-(l-2C50)*ii3. 85/96. 

•  P.  :TE(2  5,4>)X.E(I  ) 

230  CONTINUE 

CLOSE' 'JN:t-2  2  } 

C  LCSS(UNIT-23) 

C  LOSS(  J.NIT-24  ) 

CLOSE(jNIT-25  ) 


ooo  ooooooo 


PROGRAM  SLIT 


THIS  PROGRAM  CALCULATES  THE  INTENSITY  IN  THE 
COHERENT  IMAGE  OF  A  SLIT  THROUGH  AN  ABERRATED 
AND  APODISED  SYSTEM 


INTEGER  IrfK(25000) 

REAL  UK(250CC) ,H(4096) , INT(2C40 ) , E(700 ) , X (5CC ) 
REAL  V(700) 

COMPLEX  G(2C48),P(4C96) 

0  PEN (UN  IT-22,  DEVI CE-’DSK* , ACCESS- ' SEOOUT ‘  , 

>  FILE-'PH.DAT* ) 

0  PEN (UN  IT-2  3. DEVICE- *DSK ‘ , ACC  ESS  - ' S EQO UT '  , 

\  FILE-'MOD.  DAT*  ) 

0PEN(UNIT-24, DEVICE- ‘DSK* , ACCESS- * SEQOUT  *  , 

>  FILS-'B.DAT’ ) 

0  PEN (UN  I T-2 5. DEVICE- ’DSK* .ACCESS -’SEQOUT '  , 
y  FILE-'BB.DAT* ) 

READ  Parameter  values  from  mic  file 

READ(5,2C)N 
READ(5,40)a2C 
READ(5,4C)a21 
READ(5  .  4C  )a22 
REaD(5 ,4C)a3> 

READ(5 , 4C )a32 
READ(5 ,40)a42 
HEaD(5 , 4C  )gg 


2C  ?0RMaT(I4) 

4C  ?0RMAT(F>C,4) 


4  1 

fcrmat( ?>c 

.4, 

3X. 

FIC . 4  1 

4  2 

?GHMaT( I  3. 

3X, 

•  T 

♦  >  f 

3X,?>0. 4; 

43 

?ORMaT(3(2X, 13) , 

5(3X,F5. 4)) 

44 

F0RMAT(I5. 

3X. 

FA.C 

.4) 

4  “ 

?OHMAT( >2X 

.  ?  *' 

7  T  7  >  r  '1 

46 

r 

FORMAT(I4 , 

3X, 

FVC 

•  4,  4, 

4) 

n 

N -NUMBER  CF 

SA 

A# 

•  •  •  w 

^  POINTS  IN  ONE  DIM 

ENS  i: 

C  A2C-45  deg  astigmatism 

C  a2>-FOCaL  shift 


:  A22-C  deg  ASTI  IMATISV. 

C  A3>-X  COMA 

C  A32-Y  COMA 

C  A42-SPHERICAL 

:  gg-«:dth  of  Gaussian 

PI -3. >4>5927 


176 


C 

C  INPUT  FUNCTION  •  A  SLIT 
C 

DO  6C  1»».2CC 
G(I  )-CMPLX(*.  ,C.  ) 

60  CONTINUE 

DO  70  I-20>,2048 
G(I  )-CHPLX(0. ,C.) 

7C  CONTINUE 

C 

C  PERFORM  FIRST  FOURIER  TRANSFORM 
C 


C 

9C 

C 

C 

c 

c 


>cc 


CALL  FFTCC(C,N. IWK.WK) 

DO  90  I->,2C40 

H(I )-(Ca3S(G(I )))**2. 

WRITE(23,46  )I  ,K(I  ) 

CONTINUE 
GO  TO-  23> 

MOVE  ZERO  SPATIAL  FREQ  TO  CENTER  OP  APERTURE 

C(255)-G(V28) 

DO  VOO  I->,V27 
G(>27*I  )*G(I  ) 

C  (I  )-0(>92MI ) 

CONTINUE 


w- 

C, 


DO  >>C  I-V,2C48 

.  .  ^ 

C 

WRITE(22,46  )I ,C(I  ) 

CONTINUE 

C 

C  1  - 

D  APERTURE  rfITH  VARIABLE  TRANSMITTANCE  AND 

■.-.l-v 

C  A3 

SRRATIONS 

•  * 

DO  2CC  I-».255 

■**.‘“***. 

I •; I ->28)/v27. 

AH02*-A2>*(-‘»,»2.*Y*Y)/2. 

ARG22-A22*  , 

rr 

ARG32-a32*(-2.»Y»3.*Y**3)/3. 

AP.G42-A42*{>.  -6.»T*Y*6.*T**4  )/6. 

A?.0<  AHG2'*  ^=022  *  A  P.G32  ♦  AHGic  )  *2  .  ♦?  : 

:  CA 

■JSSIAN  APODISSH 

A-EXP(-(GG*Y*Y  )  ) 

?(  :  )-:MFLa  ;  A*  cos  (  ARG1  )  ,  A*3  IN  (  ARGO  )*G  ( I  ) 

P 

2:c 

CONTINUE 

DO  22C  I«256,2C48 

•7  '  •  '  «  '  V  B  •  »  '  '  "  ' 

•  * 

22C. 

CONTINUE 

DO  23C  I->.2C45 

‘-S' 

.  •  -  • 

C 

4FITE(24,46)I,F(i) 

2’C 

CONTINUE 

177 


C 

C  PERFORK  SECOND  FAST  FOURIER  TRANSFORM 
C 

N-4096 

23^  CALL  FFTCC  (F.N.WK.WK) 

DO  24C  I-%,4096 

Hd  )-((CABS{F(l  ))/>E»3)**2.  )/4.C933 
C  lrfRITE(25.46)l,H(l  ) 

240  CONTINUE 

C 

C  CALCULATE  INTENSITY  :  USED  TO  NORMALIZE 
C 

DO  244  I->,V2C 

H(I*4  99)-((CaBS(F(I ))/>E*3)**2. )/4.  >92C 
244  CONTINUE 

DO  246  I*>,499 

H(I)-((CABS(F(35  97*I ) ) /> E ♦ 3  ) **2 . ) /4 . > 92C 
246  CONTINUE 

DO  248  I»>,602 
V(I  )-  (l-3Cv)»?I/%6. 

WHITE(24.40V(I  ),H(I  ) 

248  CONTINUE 

C  GO  TO  333 

C 

C  GEOMETRIC  SLIT 
C 


26C 

DO  26C  I»%,607 

E(I  )-0. 

CONTINUE 

27C 

DO  27C  1-97, 

S(I  )->. 

CONTINUE 

DO  28C  I ->,602 

V(I  )-(l.3C>)*?I/li6. 
•  H :te(3  .  -1  >  'v . E'l 

C 

4RITS(22,4>)V(I ) ,E( 

2  3C 

CONTINUE 

333 

CL0SE(UNIT-22) 

>  1  CC 

CLOSE(UNIT-23  ; 
CL0SE(UNlT-24 ) 
CL0SE('JNIT-25) 

STOP 

END 

APPENDIX  2  THE  FABRICATION  OF  GAUSSIAN  AMPLITUDE  APODISERS 


The  Gaussian  apodisers  used  in  this  experiment  were  produced  using  a 
photographic  method  perfected  by  Weller^?.  The  method  was  origirsated  in  1964,  by  P. 
Jacquinot  and  B.  Roizen-Dossier^^  and  refined  by  Hee^o  (1974),  Kay*^  (1976),  and  Krisl82 
(1978).  The  essence  of  the  method  was  the  imaging  of  a  rotating  mask  onto  a 
photographic  plate  The  mask  is  mostly  black  with  a  generally  wedge-shaped  region  in 
white  (see  Fig  A2. 1)  The  tip  of  the  wedge  was  at  the  center  of  rotation  and  the  angular 
extent  of  the  wedge,  9,  varied  with  radial  distance,  r,  from  the  center.  Thus  the  exposure 
of  the  photographic  plate  varied  radially,  if  the  shape  of  the  wedge  was  chosen 
appropriately  the  exposure  would  be  Gaussian  By  working  in  the  linear  region  of  the 
film,  the  developed  plate  had  an  amplitude  Gaussian  transmittance. 

While  this  method  proved  satisfatory,  it  should  be  noted  that  other  methods  of 
achieving  a  Gaussian  apodiser  are  possible.  A  method  suggested  by  G.M.  Morris83 
involves  illuminating  a  film  plate  with  a  laser  beam  directly  If  the  laser  is  operating  m 
tne  owest  order  transverse  mode,  the  illumination  would  have  had  a  Gaussian  profile 
Another  method  uses  a  digital  film  plotter  such  as  the  Photowrite  P-1500  by  Optronics 
'r'e'navo''a\.  'nc  ’’’his  device  can  focus  a  light-emitting  diode  light  source  onto  a  ’2  5 
micron  diameter  spot  on  unexposed  film  The  light  can  be  modulated  (256  energy  levels) 
and  the  spot  moveo  ’o  any  location  on  j  ’O-mch  bv  tO-mch  area.  3y  orooerly  tailoring 
the  write  instructions  to  this  device,  almost  any  desired  transmittance,  including  a 


Gaussian,  can  be  achieved. 


The  Procedure  for  Making  the  Gaussian  Filter 


The  objective  was  to  make  a  transmitting  filter  which  had  an  amplitude 
transmittance  of  the  form 


(A2.1) 


where  r  was  the  radial  distance,  a  is  a  constant  which  controls  the  width  of  the  Gaussian, 
and  the  subscript  A  refers  to  the  amplitude  transmittance  of  the  filter.  In  making  the 
mask,  a  was  chosen  so  that  T/^ir)  =  0  05,  where  r  =  1 .0  inch.  When  a  was  evaluated  and 
substituted  into  (A2. 1)  the  result  was 


>(')= 


•f-)’ 

\0.S8/ 


(A2.2] 


The  intensity  transmittance  of  the  filter  was 


_rj 

'o.n 


(A2.3) 


where  the  subscriot  I  refers  to  the  intensity  transmittance  of  the  filter 

In  the  linear  region  of  the  film,  tne  intensity  transmittance  will  oe  proporuonai  to 
the  exposure.  Thus  the  opening  angle  6  which  controls  the  exposure  of  the  film  varied  as 


e|rj  =  Bjl-e  I  . 

whe-’e  B  was  chosen  such  that  0(1)  =  ru2.  Choosing  ‘his  value  of  9  insured  that  the  mas< 
variation  occurs  only  in  the  first  quadrant  The  other  quadrants  of  the  mask  were 


entirely  black.  Upon  evaluating  B,  (A2.4)  yielded 


0^^rj  =  1.58  1-e 


Using  the  Euler  relationships 

x=rsinQ  and  >=rcos0,  (A2.6) 

the  x-y  coordinates  of  the  outline  of  the  mask  were  determined. 

The  mask  profile  was  calculated  on  a  OECIO  computer  and  plotted  on  a  Tektronix 
4662  plotter.  The  profile  was  then  transferred  to  a  diffusely  .refleaing  white  paper; 
standard  computer  paper  worked  well.  The  mask  shape  in  the  first  quadrarrt  was  cut 
carefully  with  a  razor  blade  and  mounted  with  paper  cement  to  a  very  black  piece  of 
paper  The  resultant  mask  is  shown  in  Fig.  A2. 1 

The  mask  was  then  centered  on  a  turntable  which  had  been  centered  underneath 
a  camera.  The  setup  as  well  as  the  exposure  and  development  times  are  the  same  as 
those  used  by  Weller^a  The  process  which  initially  appeared  straightforward  was,  in 
fact,  time-consuming  and  produced  results  which  were  quite  variable.  Nevertheless, 
several  satisfactory  filters  were  produced 

Two  filters  (Nos.  6  and  52)  were  charactenzec  oy  rneasunng  the  Diffraction  pattern 
produced  by  them  The  optical  configuration  is  shown  m  Fig  5  1  Initially  the  apodiser 
was  not  in  the  optical  system.  The  lens  was  a  Jaegers  #1158  telescope  doublet  ana 
lens  L2  was  a  Jaegers  #958.  The  iris  haa  an  aoerture  diameter  of  0.4  inches  for  apodiser 
#  52  and  0.525  inches  for  apodiser  #  6  in  these  cases,  the  system  was  essentailly 
unaberrated.  The  aberrations  were  measured  with  a  ooint  diffraction  interferometer  as 
exolained  in  Chaoter  S  and  found  to  be  on  the  order  of  0  1  waves  for  third  order 
astigmatism,  coma,  and  spherical,  and  much  less  for  the  higher  order  aberrations 


The  results  when  these  filters  were  used  in  the  system  are  displayed  in  Figs.  A2.2  and 
A2  3  for  filters  #6  and  #52,  respectively  In  each  figure  the  photographs  are  of  the 
impulse  response  of  the  system  both  unapodisedfon  the  left)  and  apodisedfon  the  right). 
In  each  case  the  ringing  nature  r*  unapodised  impulse  response  was  suppressed  by 
the  use  of  the  apodiser.  In  th:-  :  er  half  of  each  figure,  theoretical  and  experimental 
data  are  displayed.  The  solid  lines  are  the  prediaed  values  of  a  central  slice  of  the 
irradiance  impulse  response  with  the  measured  aberrations.  The  experimental  data  are 
shown  as  dashed  curves  on  these  same  plots.  The  experimental  data  are  from  the  l-SCAN 
device  which  was  centered  on  the  impulse  response.  The  apodised  impulse  response  was 
measured  in  two  orthogonal  directions. 

The  agreement  between  theory  and  experiment  is  within  1%  for  these  data.  This 
agreement  indicates  that  the  filters  are  indeed  Gaussian  in  transmittance  and  are 
effective  in  removing  the  ringing  structure  from  the  irradiance  impulse  response  in  the 
presence  of  a  small  amount  of  aberration. 


The  characterization  of  filter  #6.  The  photographs  are  of  the 
jnapcoised(upper  ,eftj  ano  apooived  mpuise  responses  of  the  test  iystem. 
Below  each  photograph  are  plots  of  the  corresponding  thcoreticalfsolid 
curve)  and  experimental  plots  of  irradiance  along  a  line  through  the  center  of 
the  impulse  responses.  Experimental  data  are  from  two  orthogonal  slices. 
The  ordinate  of  each  piot  .s  tn  terms  of  the  ioganinm  of  relative  irradiancs  (I) 
and  the  abscissa  is  in  terms  of  the  normalized  distance  v. 
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APPENDIX  3  DATA  COLLECTION  SYSTEM 


The  data  collection  system  consisted  of  three  main  units,  as  shown  in  Fig.  A3.1. 
The  first  unit  was  the  1-SCAN  card  which  contained  a  256-element  linear  CCD  array  and 
some  supporting  electronics.  The  sec  :  s  the  data  acquisition  unit  which  converted 

the  analog  output  from  the  l-SCAN  into  an  8-bit  digital  format  and  stored  it  in  temporary 
memory.  The  last  unit  was  the  microcomputer,  which  controlled  the  data  collection 
process  and  transfered  the  data  from  the  temporary  location  to  an  internal  floppy  disk 
The  microcomputer  performed  some  simple  preprocessing  on  the  collected  data  as  well 
and  transferred  of  the  data  to  a  large  mainframe  computer  for  a  more  detailed  analysis. 


(•SCAN 

Detector 

Array 


Pig  A3  1  Block  diagram  of  the  data  collection  system 

The  i-SCAN  consisted  of  a  Fairchild  CCD1 1 1  256-element  line  scan  sensor  mounted 
on  a  orinted  circuit  card  that  contained  all  the  necessary  CCDH’  ooerating  electronics. 
The  data  were  available  on  one  of  the  output  pi  ns  of  the  card  edge  connector  m  the  form 
of  3  0.5-2.C  MHz  analog  ^Viaeo)  data  st'eam  Also  availaoie  was  a  comcensation  iCOMP' 
signal  which  was  identical  to  the  video  signal  except  for  the  lack  of  input  ooticai 
information  The  other  two  signals  used  were  a  0. 5-2.0  MHz  clock,  (MC2)  and  an. 
exaosure  sync  53)  'hese  four  sigrais  as  weil  as  oower  and  ground  .vere  connected  :o 
the  data  acquisition  unit. 


£ 
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The  purpose  of  the  data  acquisitior>  unit  was  to  condition  the  video  signal,  convert 
it  to  8-bit  digital  words  and  to  temporarily  store  it  in  a  memory  chip.  The  temporary 
storage  was  necessary  because  the  data  flowed  from  the  t-SCAN  at  a  rate  which  was 
much  faster  than  the  -  3uter  could  handle.  Additionally,  this  card  contained 

several  logical  circuits  to  coniroi  these  processes. 

The  schematic  for  the  data  acquisition  unit  is  shown  in  Fig.  A3  2.  The 
interconnections  to  the  l-SCAN  card  are  shown  m  the  upper  left-hand  corner.  Below  that 
the  power  supply  for  both  the  l-SCAN  and  the  data  acquisition  card  is  shown.  The  power 
supply  was  physically  located  within  the  data  acquisition  unit. 

The  video  and  compensation  signals  were  connected  through  47  ohm  resistors  R5 
and  R4  to  pins  6  and  5,  respectively,  of  the  Analog  Devices  ADLH0032CG  operational 
amplifier.  With  these  connections  the  operational  amplifier  acted  as  a  differential 
amplifier.  Feedback  resistors  R3.  R6,  and  R7  controlled  the  gain  of  the  op  amp.  The  gain 
was  set  to  provide  a  0.0  to  -5.0  volt  signal  on  the  input  (pin  19)  of  the  Analog  Devices 
MAH-0801-1  analog-to-digital  (A/D)  converter. 

The  A/D  converter,  when  triggered  by  the  encode  command  pulse,  (inout  on  pm 
3),  converted  the  analog  signal  into  an  eight-bit  digital  word.  When  the  data  'eady 
Signal  ipin  1  of  the  A/D)  went  low  a  valid  data  word  existed  on  the  outout  oins  5  'hrougn 
8  and  1 1  through  14 

A  two-stage  ouise  shaomg  circuit  ('J2A  and  'J2B)  generatea  the  encode  commanc 
pulse  in  response  to  the  applied  IVIG7  signal  The  MG2  signal  was  at  the  same  data  rate 
as  the  video  "he  first  stage.  U2A,  created  a  long  pulse  wmch  was  adjustaoie  m  length 
tnrougn  resistor  t '  The  second  state.  J23,  generated  the  encoae  command  ouise  wmcn 
was  approximately  100  nsec  long  The  encode  command  occurred  in  time  just  after  the 


start  of  the  data  pulse  on  pin  19  of  the  A/D  converter.  Resistor  R1  of  the  first  stage  was 
used  to  insure  this  coincidence. 

The  digital  data  generated  by  the  A/D  converter  was  connected  to  the  data  bus 
through  the  bus  transciever  U7(74LS24S).  This  chip  allowed  the  data  to  be  transferred  to 
the  memory  chip  U1 1{TC5517APR)  when  the  read  enable  signal  (from  the 
microcomputer)  was  low.  However,  when  data  was  being  transferred  from  the  memory 
chip  to  the  microcomputer,  the  read  enable  line  was  brought  high,  isolating  the  A/D 
converter  from  the  data  bus. 

The  memory  chip  required  several  time-coincident  signals  to  function  properly  in 
its  read  and  write  modes.  During  the  write  mode,  data  from  the  A/D  converter  was 
stored  into  the  chip.  When  valid  data  were  on  tne  data  bus.  the  read/write  signal  (R/W) 
on  pin  21  was  held  low  for  approximately  500  nsec.  This  negative  pulse  was  generated  by 
U6A  in  response  to  the  data-ready  signal. 

The  particular  address  within  the  memory  chip  where  each  data  word  was  stored 
was  controlled  by  the  address  counters  U8.  U9,  and  UlO.  These  binary  counters  (93L16) 
sequentially  counted  from  binary  zero  to  binary  2048  m  response  to  the  data-ready 
transitions.  The  first  address  iZe^oi  was  set  oy  the  .master  reset  (M.®)  signal  from  the 
microcomputer  before  data  acquisition  started  Then,  after  each  memory  write 
ooeration,  a  puise  from  U68  stepped  the  memory  counters  to  the  next  aodress  U63  was 
'riggered  bv  the  data-ready  transition 

The  MR  signal  not  only  reset  the  memory  aodress  counters  to  zero,  but  it  also  reset 
the  three-input  logical  switch  U4  When  U4  was  reset  the  output  on  pm  10  was  held  low 
This  pin  was  connected  to  one  of  the  inputs  of  the  AND  gate  U3A,  thus  when  it  was  low 
the  encode  command  pulse  could  not  reach  the  A/D  converter  and  data  were  not 


converted. 


The  data  collection  process  was  initiated  by  a  MR  followed  by  a  carriage  return 
(CR)  signal  from  the  microcomputer.  The  latter  signal  set  U4  so  that  the  next  exposure 
sync  (ES)  signal  from  the  l-SCAN  brought  the  output  (pin  10)  of  U4  high.  These  ES  pulses 
occurred  in  time  just  before  a  serial  stream  of  analog  data.  The  ES  pulses  were  about  1 
msec  apart  and  the  data  stream  was  about  SOO  nsec  long  (for  an  0.5  MHz  data  rate). 

The  data  collection  process  up  to  this  point  can  be  summarized  in  the  following 
series  of  steps. 

1.  A  MR  from  the  microcomputer  reset  the  three-input  switch  U4  and  zeroed  the 
memory  address  counter. 

2.  A  CR  set  U4. 

3.  The  next  ES  forced  pin  10  of  U4  high,  allowing  the  encode  commands  to  pass 
through  the  AND  gate  U3A. 

4.  Each  encode  command  initiated  an  A/D  conversion  which  was  stored  in  the 
memory  chip, 

5.  The  data-ready  transition  generated  two  pulses  The  first  signalled  the  memory 
chip  to  store  the  data  currently  on  the  data  bus  After  the  oata  word  was  stored, 
the  second  oulse  incremented  the  memory  adcress  for  the  next  data  word 

There  were  256  (2^)  detector  elements  on  the  CCDl'l  line  scan  detector  "his 
constituted  one  frame  of  data.  When  the  last  address  of  a  frame  was  reached  the 
memory  address  counters  generated  a  terminal  count  (TC)  signal  which  was  useo  to 
terminate  the  data  collection  process  for  that  frame  The  memory  chip  could  hold  eight 
frames  so  the  process  was  repeated  seven  times.  The  TC  aaed  as  a  memory  reset  for  U4 
but  did  not  reset  the  memory  address  counters  Instead,  a  software  clock  signal  from  the 
microcomputer  incremented  the  address  counters  by  one  address.  Then  a  CR  initiated 


another  frame  of  data.  The  microcomputer  monitored  the  TC  signal  and  after  the  eighth 
one  was  sensed,  the  data  collection  process  was  terminated. 

Next,  the  data  were  transferred  to  the  microcomputer.  The  first  signal  in  this 
process  was  a  MR  which  set  the  memory  address  to  binary  zero.  Next,  read  enable  was 
brought  high  isolating  the  A/D  converter  and  forcing  the  memory  chip,  through  the 
signal  on  pin  20,  into  the  data-read  mode.  At  this  time,  the  contents  of  the  first  memory 
location  were  on  the  data  bus.  When  the  microcomputer  had  stored  that  data,  it  issued  a 
software  clock  pulse.  This  pulse  incremented  the  memory  address  thus  placing  the  next 
word  of  data  on  the  data 'bus.  This  process  continued  until  all  2048  memory  locations 
were  read  and  stored  in  the  microcomputer. 

The  microcomputer  was  a  Healthkit  H89with  64K  (8-bit)  internal  memory  and  two 
100K  single-density,  single-sided  disk  drives.  Once  all  eight  frames  of  data  were  within 
the  microcomputer,  a  software  program  on  the  microcomputer  performed  some  simple 
preprocessing.  First,  the  eight  frames  of  data  were  averaged,  on  a  pixel  by  pixel  basis 
Then  the  data  was  displayed  in  a  x-y  graphic  format  on  the  microprocessor  CRT  screen. 
Finally,  the  average  frame  was  stored  in  a  file  on  disk  Additional  software  allowed  one 
frame  of  data  to  be  subtraaed  from  another  (SUB.  DAT)  and  the  data  to  be  transferred  to 
a  OEC10  comouter  (via  MODEM7)  The  data  collection  process  was  initiated  bv  entering 
DATA  FiLE.DAT'  into  the  microprocessor  The  data  was  stored  in  the  file  named 
=iLE  DA”  ”he  software  jsea  ;n  this  data  jccjisition  orocess  was  develooed  bv  Mr  .ee 


